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Abstract 

Dispersive and Strichartz estimates for solutions to general strictly hyper- 
bolic partial differential equations with constant coefficients are considered. 
The global time decay estimates of — L'^ norms of propagators is discussed, 
and it is shown how the time decay rates depend on the geometry of the prob- 
lem. The frequency space is separated in several zones each giving a certain 
decay rate. Geometric conditions on characteristics responsible for the partic- 
ular decay are investigated. Thus, a comprehensive analysis is carried out for 
strictly hyperbolic equations of high orders with lower order terms of a gen- 
eral form. Results are applied to time decay estimates for the Fokker-Planck 
equation and for semilinear hyperbolic equations. 
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1 Introduction 

These notes are devoted to the investigation of dispersive and Strichartz estimates 
for general hyperbolic equations with constant coefficients. The analysis that we 
carry out is also applicable to hyperbolic systems either by looking at characteristics 
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of the system directly, or first taking the determinant of the system (the dispersion 
relation). 

There are several important motivations for the analysis. First, while hyperbolic 
equations of the second order (such as the wave equation, dissipative wave equation, 
Klein-Gordon equation, etc.) are very well studied, relatively little is known about 
equations of higher orders. At the same time, equations or systems of high orders 
naturally arise in applications. For example, Grad systems of non-equilibrium gas 
dynamics, when linearised near an equilibrium point, are examples of large hyperbolic 
systems with constant coefficients (see e.g. |Rad03] . |Rad05] ). Here one has to deal 
with hyperbolic equations of orders 13, 20, etc., depending on the number of moments 
in the Grad system. Moreover, there are important families of systems of size going 
to infinity, or even of infinite hyperbolic systems. For example, the Hermite-Grad 
method for the analysis of the Fokker-Planck equation for the distribution function 
for particles for the Brownian motion produces an infinite hyperbolic system with 
constant coefficients. Indeed, making the decomposition in the space of velocities 
into the Hermite basis, and writing equations for the space-time coefficients produces 
a hyperbolic system for infinitely many coefficients (see e.g. |VR03j . [VR04j . |ZR04j . 
and Section [8.51) . The Galerkin approximation of this system leads to a family of 
systems with sizes increasing to infinity. Although explicit calculations are difficult in 
these situations, the time decay rate of the solution can still be calculated ( |Ruzh06] ). 

One of the main difficulties when dealing with large systems is that unlike in the 
case of the second order equations, in general characteristics can not be calculated 
explicitly. This raises a natural problem to look for properties of the equation that 
determine the decay rates for solutions. On one hand, it becomes clear that one has 
to look for geometric properties of characteristics that may be responsible for such 
decay rates. On the other hand, a subsequent problem arises to be able to reduce 
these properties from some properties of coefficients of the equation. 

One encounters several difficulties on this path. One difficulty lies in the absence 
of general formulae for characteristic roots. For large frequencies one can use per- 
turbation methods to deduce the necessary asymptotic properties of characteristics. 
However, this approach can not be used for small frequencies, where the situation 
becomes more subtle. For example, for small frequencies characteristics may become 
multiple, causing them to become irregular. This means that if we use the usual 
representation of solutions in terms of Fourier multipliers, phases become irregular, 
while amplitudes are irregular and blow up. Thus, we will need to carry out the de- 
tailed analysis of sets of possible multiplicities using the fact that they are solutions 
of parameter dependent polynomial equations. Another difficulty for small frequen- 
cies is that there exists a genuine interaction between time and frequencies. In the 
case of homogeneous symbols it can be shown (see e.g. Section II. 2p that time can 
be taken out of the estimates, after which low frequencies can be ignored since the 
corresponding operators are smoothing and their estimates are independent of time. 
In the case of the presence of lower order terms, the time can no longer be eliminated 
from the estimates, so even small frequencies become large for large times and may 
infiuence the resulting estimates. 
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The purpose of this work is to present a comprehensive analysis of such problems. 
Despite the difficulties described above, we will be able to determine what geomet- 
ric properties of characteristic roots are responsible for qualitatively different time 
decay rates for solutions. Moreover, we will calculate these rates and relate them 
to geometric properties of equations. This will lead to a comprehensive picture of 
decay rates and orders in dispersive estimates for hyperbolic equations with constant 
coefficients. Such estimates lead to Strichartz estimates, for which our analysis will 
be applied, with further implications for the corresponding semilinear problems. 

Thus, in this paper we consider a problem of determining dispersive and Strichartz 
estimates for general hyperbolic equations with lower order terms. Therefore, we 
consider the Cauchy problem for general m^^ order constant coefficient linear strictly 
hyperbolic equation with solution u = u(t,x): 

homogeneous principal part general lower order terms 

' " ^ ' " ^ 

rn m~ 1 

< Drw+$^^.(^.)A'""'^+5Z c^,rD:Diu = o, t>o, (1.1) 

i=l 1=0 \a\+r=l 

^ Dlu{0, x) = fi{x) G Co°°(M"), / = 0, . . . , m - 1, X G M" , 

where Pj{0^ the polynomial obtained from the operator Pj{Dx) by replacing each D^^ 
by C,k, is a constant coefficient homogeneous polynomial of order j, and the Ca,r are 
(complex) constants. Here, as usual, a = («i, . . . , a„), = D'^^ ■ ■ ■ -D"^, D^f, = \dxf, 
and Dt = ^dt- The full symbol of the operator in (11. ip will be denoted by 

m m—l 
j=l 1=0 \a\+r=l 

where = C,i^ ■ ■ ■ . We will always assume that the differential operator in (II. ip 
is hyperbolic, that is for each G M", the symbol of the principal part, 

m 

i=i 

has m real roots with respect to r. For simplicity, unless explicitly stated otherwise, 
we will also assume that the operator in (II. ip is strictly hyperbolic, that is at each 
^ G M"\{0}, these roots are pairwise distinct. We denote the roots of Lmij,^) with 
respect to r by <y5i(^) < ■■■ < ipm{i), and if L is strictly hyperbolic the above 
inequalities are strict for ^ 7^ 0. 

The condition of hyperbolicity arises naturally in the study of the Cauchy problem 
for linear partial differential operators and it can be shown that it is a necessary con- 
dition for well-posedness of the problem; this is discussed in |ES92] and |Hor83b] , 
for example. Strict hyperbolicity is sufficient for well-posedness of the Cauchy 
problem for such an operator with any lower order terms; if the operator is only 
hyperbolic (sometimes called weakly hyperbolic) the lower order terms must satisfy 
additional conditions for well-posedness, the so-called Levi conditions. For this 



4 



reason, we only consider strictly hyperbolic operators with lower order terms, since 
our main interest is to understand the influence of lower order terms on the decay 
properties of solutions. 

The roots of the associated full characteristic polynomial L{t,C,) with respect to 
r will be denoted by ri(,^), . . . , rm(0 and referred to as the characteristic roots of 
the full operator. Clearly, if L is a homogeneous operator then the characteristic 
roots Tk{^), k = 1, . . . ,m, coincide, possibly after reordering, with the roots (PkiO^ 
k = 1, . . . ,m, oi the operator L^. However, in general there is no natural ordering 
on the roots Tfc(^) as they may be complex- valued or may intersect. 

The analysis here will be based on the properties of characteristic roots Tfc(^). If 
the problem (II. ip is strictly hyperbolic, we can derive their asymptotic properties 
in a general situation, necessary for our analysis. However, if the problem is only 
hyperbolic, functions Tfc(^) may develop singularities for large ^. If this does not 
happen and we have the necessary information about them, we may drop the strict 
hyperbolicity assumption. This may be the case in some applications, for example in 
those arising in the analysis of the Fokker-Planck equation. 

We seek a priori estimates for the solution u{t,x) to the Cauchy problem (II. ip . 
of the type 



where l<p<2, ^ + ^ = 1, A'p = Np{a,r) is a constant depending on p,a and r, 

and K{t) is a function to be determined. Here W^'' ' is the Sobolev space over 
with Np — I (fractional) derivatives. 

We note that sometimes, for example in [TreSOj . in the definition of a hyperbolic 
operator the polynomial L(ir, ^) is used as it is better suited to taking the partial 
Fourier transform in x, corresponding as it does to L{dt, Dx); in this case, one requires 
the roots with respect to r to be purely imaginary (in the cases when we will require 
them to be real). However, the definition that we give above is perhaps more standard, 
and thus adopted here throughout. 

For a hyperbolic equation with real coefficients we note that the constants Ca^r 
satisfy ^'""'"'"'cQ.^r ^ ^] the equation is written in the form above since our results 
may be used to study hyperbolic systems, which can be reduced to an m^^ order 
equation with complex coefficients. 

Most results presented here will apply to operators which are pseudo-differential 
in X and to hyperbolic systems via their dispersion equation. Moreover, most of 
results in this paper are in general sharp. 

In this work, we place the priority on obtaining a comprehensive collection of 
estimates for hyperbolic equations with constant coefficients. The case of variable 
coefficients is also of great interest, but we leave some extensions of our analysis to 
this case outside the scope of this paper. Let us mention that already in the case of co- 
efficients depending on time, some unpleasant phenomena may happen. For example, 
already for the second order equations the oscillations in time dependent coefficients 
may change the time decay rates for solutions to the corresponding Cauchy problem. 





1=0 
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For example, equations with very fast oscillations, or with increasing coefficients, 
have been analysed in |RY99l [RYOOj . to mention only a few references. Results even 
for the wave equations with bounded coefficients may depend on the oscillations in 
coefficients (see e.g. |ReS05j ). At the same time, many results of this paper are stable 
under time perturbations of coefficients. For example, in the case of equations with 
homogeneous symbols with time-dependent coefficients with integrable derivative, a 
comprehensive analysis has been carried out in [MR07j . We will not deal with such 
questions in this paper. Let us also mention that while dispersive estimates are de- 
voted to W — L'^ estimates for solutions, — U' estimates are also of interest. A 
survey of W estimates for general non-degenerate Fourier integral operators and their 
dependence on the geometry can be found in |RuzhOO] in the case of real- valued phase 
functions, while operators with complex-valued phase functions have been analysed 
in |Ruzh01] . L^-estimates for solutions to some classes of hyperbolic systems with 
variable multiplicities appeared in [KR07j . 

Let us now explain the organisation of these notes. In the following parts of the 
introduction we will review results for second order equations and for equations with 
homogeneous symbols, as well as give several more motivations for the comprehensive 
analysis of this paper. In Section [2] we will present results for different types of 
behaviour of characteristic roots, and also of corresponding phase functions in cases 
where we can represent solutions in terms of Fourier multipliers. Thus, in Section 
12.11 we will present results without and with multiplicities, when roots are separated 
from the real axis, in which cases we can get exponential decay of solutions. In 
Section we present results for roots with non- degeneracies, in which case we have 
a variety of conclusions depending on geometric properties of roots. In Section 12.31 
we present results for complex roots that become real on some set. A version of this 
type of statements (although not in the microlocal form used here) partly appeared 
in |RS05j . and those are improved here. In Section YTM we summarise the microlocal 
results and formulate the main theorem on dispersive estimates for general hyperbolic 
equations with constant coefficients. Theorem 12.181 is the main theorem containing a 
table of results, and the rest of this section is devoted to the explanation and further 
remarks about this table. In Section 12.51 we will outline our approach, indicating 
the relations between frequency regions and statements. In Section 12.61 we present 
results for non- homogeneous equations, as well as formulate corresponding Strichartz 
estimates with further applications to semilinear equations. In general, we leave such 
developments outside the scope of this paper since they are quite well understood (see 
e.g. [KT98] ). once the time decay rates are determined (as we will do in Theorem 

EUD. 

The subsequent chapters contain the detailed analysis and proofs. In Section [3] we 
establish necessary properties of roots of hyperbolic polynomials, as well as carry out 
the perturbation analysis for large frequencies. In Section |4] we investigate estimates 
for oscillatory integrals under certain convexity assumptions on the level sets of the 
phase function. In Section [5] we analyse the corresponding oscillatory integrals with- 
out convexity assumption. Section E] is devoted to dispersive estimates for solutions 
to the general Cauchy problem, and here we prove various parts of Theorem I2.18[ 
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Section [7] deals with multiple characteristics. Here we present a procedure for the 
resolution of multiplicities in the representation of solutions, enabling us to obtain 
estimates in these cases as well. Section 17.41 is devoted to multiple roots on the real 
axis. Here, we investigate solutions for frequencies very close to multiplicities (in 
some shrinking neighborhoods) as well as for larger, but still bounded, frequencies. 
Here we present several different versions of results dependent on possibly different 
assumptions. Finally, Section [S] is devoted to examples of the presented analysis with 
further applications. Thus, in Section 18.11 we deal with second order equations and 
give examples of how our results can be applied to investigate the interplay between 
mass, dissipation, and frequencies. Further, in Section [8^ we discuss some conditions 
on coefficients of equations, and in Section [8731 we give examples of non-homogeneous 
roots in terms of hyperbolic triples and Hermite's theorem. In Section 18.41 we show 
briefly how the results can be applied for strictly hyperbolic systems. And finally, in 
Section 18.51 we give an application to the Fokker-Planck equations. 

The authors are grateful to Jens Wirth for remarks about the preliminary version 
of the manuscript and to Tokio Matsuyama and Mitsuru Sugimoto for discussions. 

We will denote various constants throughout the paper by the same letter C. 
Balls with radius R centred at ,^ G will be denoted by Br{^). We will use the 
notation {i) = ^/lTW^ (D) = v^I^A and \D\ = \ - A\^/^ . The Sobolev space Wj, 
is then defined as the space of measurable functions for which {D)^ f G L^(M"). 

We will also use the standard notation for the symbol class = S^q, as a 
space of smooth functions a = ^) G C°°(M" x M") satisfying symbolic estimates 
\d^d^a{x,0\ < Cafsil + |<^|)^"'"', for all x,^ G M", and all multi-indices a, (3. 

If function a = a(^) is independent of x, we will sometimes also write a G Siq{U) 
for an open set U C M", if a = a(0 e C°°{U) satisfies |(9f a(OI < -|- for 
all ^ G [/, and all multi- indices a. 



1.1 Background 

The study of — decay estimates, or Strichartz estimates, for linear evolution 
equations began in 1970 when Robert Strichartz published two papers, |Str70aj and 
|Str70bj . He proved that if m = u{t,x) satisfies the Cauchy problem (that is, the 
initial value problem) for the homogeneous linear wave equation 

dtu{t,x) - A^u{t,x) = 0, (t, x) G M" X (0, oo) , 
u{0,x) = (l){x), dtu{0,x) =^Ij{x), xgM", 

where the initial data (p and ip lie in suitable function spaces such as C^(M"), then 
the a priori estimate 

||(ni(t,-),V.n(t,-))||L. <C(l + t)-^(i-^)||(V.0,^)||^iv, (1.4) 

holds when n>2, ^ + ^ = 1, l<p<2 and Np > — ^). Using this estimate, 
Strichartz proved global existence and uniqueness of solutions to the Cauchy problem 
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for nonlinear wave equations with suitable ("small") initial data. This procedure 
of proving an a priori estimate for a linear equation and using it, together with 
local existence of a nonlinear equation, to prove global existence and uniqueness for 
a variety of nonlinear evolution equations is now standard; a systematic overview, 
with examples including the equations of elasticity, Schrodinger equations and heat 
equations, can be found in [Rac92j . or in many other more recent books. 

There are two main approaches used in order to prove (11.41) : firstly, one may write 
the solution to (II. 3p using the d'Alembert (n = 1), Poisson [n = 2) or Kirchhoff 
{n = 3) formulae, and their generalisation to large n, 



u{t, x) 



1 



UA (2j - 1) 



n-3 
I 2 



n-1 



'dS 



n-1 



it 



Bt{x) 



dBt(x) 

ip dS 

dBt{x) 

\y - x\ 



(odd n > 3) 



dy 



dy 



y-x\ 



even n) 



(here f- stands for the averaged integral; for the derivation of these formulae see, for 
example, |Ev98j ). as is done in [vW71j and |Rac92j . Alternatively, one may write the 
solution as a sum of Fourier integral operators: 



r — m + ^ 



This is done in |Str70aj . jBre75j and |Pec76j . for example. Using one of these repre- 
sentations for the solution and techniques from either the theory of Fourier integral 
operators ( |Pec76] ). Bessel functions ( |Str70a] ). or standard analysis ( |vW71j ). the 
estimate (11.41) may be obtained. 

Let us now compare the time decay rate for the wave equation with equations 
with lower order terms. An important example is the Klein-Gordon equation, where 
u = u{t,x) satisfies the initial value problem 



dtu{t, x) - A^u{t, x) + f?u{t, x) = 0, {t, x) G M" X (0, oo) 
u(0,x) = (t>{x), Ut{0,x) = ipi^x), xgM", 



:i.5) 



where E C^(M"), say, and /z 7^ is a constant (representing a mass term); then 

||(M(t,-),Mt(t,-),V.M(t,-))||L. <C(l + ^)-^(^i)||(V.0,7/.)||^iv,, (1.6) 



where p,q,Np are as before. Comparing (11.41) to (11.61) . we see that the estimate for 
the solution to the Klein-Gordon equation decays more rapidly. The estimate is 
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proved in [vW71] . [Pec 76 j and |Hor97j in different ways, each suggesting reasons for 
tliis improvement: in [vW71j . tlie function 



v{x,Xn+i,t) := e '^''"+^u(t,x) 



is defined; using f ll.51) . it is simple to sliow tliat v satisfies tlie wave equation in M"+^, 
and thus the Strichartz estimate (11 ■4p holds for v, yielding the desired estimate for u. 
This is elegant, but cannot easily be adapted to other situations due to the importance 
of the structures of the Klein-Gordon and wave equations for this proof. In |Pec76j 
and [Hor97j . a representation of the solution via Fourier integral operators is used 
and the stationary phase method then applied in order to obtain estimate (II. 6p . 

Another second order problem of interest is the Cauchy problem for the dissipative 
wave equation. 



d^u{t, x) - A^u{t, x) + ut{t, x) = , {t, x) G M" X (0, oo), 
u{0,x) = (f){x), ut{0,x) = 4j{x), xgM", 



:i.7) 



where tp, ^ 



say. In this case. 



\did:u{t, 



< c{i + ty 



H 

2 



W' 



with some Np = Np{n,a,r). This is proved in [Mat 77] with a view to showing well- 
posedness of related semilinear equations. Once again, this estimate (for the solution 
u{t, x) itself) is better than that for the solution to the wave equation; there is an 
even greater improvement for higher derivatives of the solution. As before, the proof 
of this may be done via a representation of the solution using the Fourier transform: 



2e-*/2 sinh(|v/l-4|e|2) 



u{t, x) 



3,-1 



^""^^^(^^^-^.e-cos(|vW-T)l0(O 



2e-*/2 sin(|v/4|e|2-l) 



+ 



Matsumura divides the phase space into the regions where the solution has different 
properties and then uses standard techniques from analysis. 

It is, therefore, motivating to ask why the addition of lower order terms improves 
the rate of decay of the solution to the equation; furthermore, in the first instance, we 
would like to understand why the improvement in the decay is the same for both the 
addition of a mass term and for the addition of a dissipative term. It will follow from 
the analysis of the paper that the quantities responsible for the decay rates for the 
Klein-Gordon and dissipative equations are of completely different nature. In the first 
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instance the characteristic roots are real and he on the real axis for all frequencies, 
while for the latter equation they are in the upper complex half-plane, intersect at a 
point, and one of them comes to the origin. From this point of view, the same decay 
rates in the dispersive estimate for these two equations is quite a coincidence. On 
the example of the dissipative equation we can see another difficulty for the analysis, 
namely the appearance of the multiple roots. This may lead to the loss of regularity 
in roots and blow-ups in the amplitudes of a representation, so we need to develop 
some techniques to deal with this type of situations. 

These questions are even more important for equations of higher orders. Let 
us mention briefly an example of a system that arises as the linearisation of the 
13-moment Grad system of non-equilibrium gas dynamics in two dimensions (other 
Grad systems are similar). The dispersion relation (the determinant) of this system 
is a polynomial of 9*^^ order that can be written as 



with polynomials Qj{uj,C,) deflned by 



Qi 



Qi 



103 



21 



13 

y 

67 



25 



-UJ 



-UJ 



912 
2625 



a/3 



25 V 



-UJ 



1094 , 1381 



9 



-UJ 



19 



75 
497 
~ ^ 
2908 



125 
3943 . 



UJ 



UJ 



-UJ 



where 



4 



225 
178 



■UJ^ + 



375 
13 2 

-UJ 



-UJ 



2032 
" 6905 
832 



432 
~ 675 
264 



3943 



125 
159 



125 



1 



143 
" 330 
48 



159 



aP 



32 
325 



aP 



6 

25 



9 



3 



-UJ 



UJ 



45 
1)^ 



-UJ 



^(0 



r(0 



a 



1^19' ^ 1^19' 



A natural question of flnding dispersive (and subsequent Strichartz) estimates for the 
Cauchy problem for operator P{Dt, D^) with symbol P(r, ^) becomes calculationally 
complicated. Clearly, in this situation it is hard to flnd the roots explicitly, and, 
therefore, we need some procedure of determining what are the general properties of 
the characteristics roots, and how to derive the time decay rate from these properties. 
Thus, in |Rad03] and [VR04] it is discussed when such polynomials are stable. In 
this case, the analysis of this paper will guarantee the decay rate, e.g. by applying 
Theorem 12.161 for frequencies near the origin. Theorem 12.21 for bounded frequencies 
near possible multiplicities (independent of the structure of such multiplicities), and 
Theorem 12.11 for large frequencies. In fact, once the behavior of the characteristic 
roots is understood. Theorem 12.181 will immediately show that the overall time decay 
rate here is the same as for the dissipative wave equation. 



10 



1.2 Homogeneous symbols 

The case where the operator in (11. 1|) has homogeneous symbol has been studied 
extensively: 

' L^{D,, Dt)u = 0, gM" X (0,oo), 

Dlu{0,x) = fi{x), / = 0,...,m- 1, X e M", 



where Lm is a homogeneous m^^ order constant coefficient strictly hyperbolic differ- 
ential operator; the symbol of Lm may be written in the form 

^m(r,0 = {r-<^i{0)...{T-^m{0), with < ■■• < <^™(0 0). 

In a series of papers, |Sug94] , |Sug96| and |Sug98| , Sugimoto showed how the geomet- 
ric properties of the characteristic roots ipi{^), . . . , (pmiO affect the — U estimate. 
To understand this, let us summarise the method of approach. 

Firstly, the solution can be written as the sum of Fourier multipliers: 



m— 1 

1=0 k=l 



c{t,x) = J2[Ei{t)fi]{x), where Ei{t) = J] g^^ie^^^^^^^afc^OS^, 



and ak,i{0 is homogeneous of order — /. Now, the problem of finding an — L'^ decay 
estimate for the solution is reduced to showing that operators of the form 

Mr{D) :=3^-v^(«)|er'^x(09^, 

where v^(0 ^ C"^(]R" \ {0}) is homogeneous of order 1 and x ^ C°°(]R") is equal to 1 
for large ^ and zero near the origin, are — L'^ bounded for suitably large r > I. In 
particular, this means that, for such r, we have 

l|i^^Ki)/IU. <c||/||^.-,. 

Then it may be assumed, without loss of generality, that t = 1. Indeed, it can be 
readily checked that for t > and / G C^(M"'), we have the equality 

mt)f]{x) = t^[E,{i)f{t-)]{t-'x). 

Using this identity and denoting ft{-) = fit-), we have 

\\m)f\\% = t'im)ft]{t-'-)\\% = t^' [ \mi)mt-'x)\'^dx 



Then, noting that a simple change of variables yields 
we have. 
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r-l 
P 



hence, 

m— 1 

\Ht,.)U.<cf-'<~I^Y.\\f^\\wr^- 

1=0 

It has long been known that the values of r for which Mr{D) is — L'^ bounded 
depend on the geometry of the level set 

S^ = {eGR"\{0}:^(O = l}. 

In |Lit73j . |Bre75j . it is shown that if the Gaussian curvature of is never zero 
then Mr{D) is — L'^ bounded when r > -^^(^ — This is extended in |Bre77j 

where it is proven that Mr{D) is Lp - bounded provided r > - -), where 

p = min^^o rank Hess (p{C)- 



Sugimoto extended this further in Sug94 , where he showed that if is convex 
then Mj.[D) is — L'^ bounded when r > {n — (p " 5)' here, 

7(S) := sup sup 7(S; (T,P) , S C M" a hypersurface , 

0-6S P 

where P is a plane containing the normal to S at cr and 7(S; a, P) denotes the order 
of the contact between the line T^- fl P, T^r is the tangent plane at a, and the curve 
S n P. See Section H73l for more on this maximal order of contact. 

In order to apply this result to the solution of (11. 91) . it is necessary to find a 
condition under which the level sets of the characteristic roots are convex. The 
following notion is the one that is sufficient: 

Definition 1.1. Let L = L{Dt,Dx) be a homogeneous tv}^ order constant coefficient 
partial differential operator. It is said to satisfy the convexity condition if the matrix 
of the second order derivatives, Hess corresponding to each of its characteristic 

roots . . . , V2m(0j scmi-dcfinite for ^ 7^ 0. 

It can be shown that if an operator L does satisfy this convexity condition, then 
the above results can be applied to the solution and thus an estimate of the form 
ffT^ holds with 

K{t) = (1 J with some 7 < m, (1-10) 

where 7 can be related to the convex indices of the level sets of characteristics. Indeed, 
under the convexity condition one can show that (pk can be made always positive or 
negative by adding an affine function, the corresponding level sets S^^ = G M" : 
4>k{0 = 1} convex for each k = 1, . . . ,m, and that 7(2,^^,) < 2[m/2]. So the decay 
in (11.101) is guaranteed with 7 = 2[m/2]. 

Finally, if this convexity condition does not hold the estimate fails; in Sug96| 



and |Sug98| it is shown that in general, Mr{D) is — L'^ bounded when r > (n 
L . ) (- — -), where 

70(E) := sup inf 7(E; cr, P) < 7(E). 
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For n = 2, 7o(S) = 7(S), so, the convexity condition may be lifted in that case. 
However, in |Sug96| , examples are given when n > 3, p = 1, 2 where this lower bound 
for r is the best possible and, thus, the convexity condition is necessary for the above 
estimate. It turns out that the case n>3, l<p<2is more interesting and is 
studied in greater depth in Sug98| , where microlocal geometric properties must be 



looked at in order to obtain an optimal result. 

Two remarks are worth making; firstly, the convexity condition result recovers 
the Strichartz decay estimate for the wave equation, since that clearly satisfies such 
a condition. Secondly, the convexity condition is an important restriction on the 
geometry of the characteristic roots that affects the — L'^ decay rate; hence, in the 
case of an m^^ order operator with lower order terms we must expect some geometrical 
conditions on the characteristic roots to affect the decay rate of solutions. 



2 Main results 

We will now turn to analysing the conditions under which we can obtain U' — L'^ decay 
estimates for the general im}^ order linear, constant coefficient, strictly hyperbolic 
Cauchy problem 

L{D,,D,)^DTu + Y,Pj{D^)DT^'u + Y, E c«,.D:A'« = 0, t > 0, 

3=1 1=0 \a\+r=l 

D[u{0, x) = fi{x) e C^{W), / = 0, . . . , m - 1, X e . 

(2.1) 

Results of this section will show how different behaviours of the characteristic roots 
Ti{^) , . . . , TmiO affect the rate of decay that can be obtained. As in the introduction, 
the symbol Pj{^) of Pj{Dx) is a homogeneous polynomial of order j, and the Ca^r 
are constants. The differential operator in the first line of fl2.ip will be denoted by 
L{Dt, Dx) and its symbol by L{t,E,). The principal part of L is denoted by Lm- Thus, 
Lm{T,^) is a homogeneous polynomial of order m. In the subsequent analysis, ideally, 
of course, we would like to have conditions on the lower order terms for different rates 
of decay; in Section [H] we shall give some results in this direction. For now, though, 
we concentrate on conditions on the characteristic roots. 

First of all, it is natural to impose the stability condition, namely that for all 
^ G M" we have 

Imrfc(O>0 for A: = l,...,m; (2.2) 

this is equivalent to requiring the characteristic polynomial of the operator to be 
stable at all points ^ G M", and thus cannot be expected to be lifted. In fact, certain 
microlocal decay estimates are possible even without this condition if the supports of 
the Fourier transforms of the Cauchy data are contained in the set where condition 
(12. 2p holds. However, this restriction is only technical so we may assume (12. 2p without 
great loss of generality since otherwise no time decay of solution can be expected. 

Also, it is sensible to divide the considerations of how characteristic roots behave 
into two parts: their behaviour for large values of |,^| and for bounded values of |,^|. 
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These two cases are then subdivided further; in particular the following are the key 
properties to consider: 

• multiplicities of roots (this only occurs in the case of bounded frequencies |^|); 

• whether roots lie on the real axis or are separated from it; 

• behaviour as |.^| ^ oo (only in the case of large |.^|); 

• how roots meet the real axis (if they do); 

• properties of the Hessian of the root, Hessrfc(,^); 

• a convexity- type condition, as in the case of homogeneous roots (Section II .21) . 

For some frequencies away from multiplicities we can actually establish indepen- 
dently interesting estimates for the corresponding oscillatory integrals that contribute 
to the solution. Around multiplicities we need to take extra care of the structure of 
solutions. This will be done by dividing the frequencies into zones each of which will 
give a certain decay rate. Combined together they will yield the total decay rate for 
solution to (12.1 1) . Several theorems below will deal with integrals of the form 

f e^(-«+-(?)*)a(Ox(e)c^e, (2.3) 

which appear in representations of solutions to Cauchy problem (12.11) as kernels of 
propagators, where a(^) is a suitable amplitude and x(0 is a cut-off to a correspond- 
ing zone, which may be bounded or unbounded. Solution to the Cauchy problem 
(12. ip can be written in the form 

m—l 
j=0 

where propagators Ej{t) are defined by 

« m 

E^{t)f{x)= / e"-«(^e^-'=(«*4(t,0)x(0/(0rfe, (2-4) 

with suitable amplitudes Aj{t,^). In the areas where roots are simple, phases and 
amplitudes are smooth, and we can analyse the sum (12. 4p termwise, reducing the 
analysis to integrals of the form (12. 3p . In the case of multiple characteristics we will 
group terms in (12. 4p in a special way to obtain suitable decay estimates. Below we 
will give results for decay rates dependent on the different qualitative behaviours of 
the characteristic roots. 
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2.1 Away from the real axis: exponential decay 

We begin by looking at the zone where roots are separated from the real axis. If the 
roots are smooth, we can analyse solution (12. 4p termwise: 

Theorem 2.1. Let t : U ^ C be a smooth function, f/ C M" open. Let a G Sil^{U), 
i.e. assume that a = a(0 e C°°{U) satisfies \d^a{0\ < C^{1 + for all 

C, E U and all multi-indices a. Let x ^ 'S'i,o(^") ■^'"'^^ ^^^^ X = outside U. 
Assume further that: 

(i) there exists 6 > such that Imr(,^) > d for all ^ eU] 

(ii) \r{i)\<C{l + \i\)foraineU. 
Then for all t > we have 

DID:( I e^(-^+^(«Wa(0x(0/(0rf^) , < Ce-'i/||^.,.H..-., (2.5) 

VJjj„ / L9(Kg) Wp 

where ^ + -^ = l,l<p<2,Np>n{^-^),r>Q,aa multi-index and f e Cg°(M"). 
If p = 1, we take Ni > n. 

Moreover, let us assume that equation L(r, ^) = has only simple roots Tfc(^) 
which satisfy condition (i) above, in the open set U C M", for all k = 1, . . . ,m. Then 
solution u to 112. 1\} satisfies 

■m— 1 

||D[D:x(^)^(t,-)||L.{Rs) <C^e-^*$^||//||^Wi+^-^, (2.6) 

/=o 

where 1 < p < 2, ^ + ^ = 1, and Np, r, a are as above. 

The proof of Theorem 12.11 will be given in Sections 16.41 and 16.101 Note also that 
if we omit assumption (ii) in Theorem 12. H estimate (12. 5p with r = still holds. In 
the case of (12. 6p . it can be shown (see Proposition 13.81) that characteristic roots of 
operator L{Dt,D^) in (12.11) satisfy (ii). 

We also note, that we may have different norms on the right hand side of (12.61) . 
For example, we will show in Section 16.41 that under conditions of Theorem 12.11 we 
also have the following estimate: 

m—l 

\\DlD:x{D)u{t, OIlL.dRs) < Ce-'' ^ 1 (2.7) 

/=o ^ 

where 1 < p < 2, i + i = 1, A^^ > f (i - i), and iV^ > f for p = 1. Estimate 
will follow from (16.81) and Proposition 16.51 by interpolation. In turn, interpolating 
between (12.61) and (12.71) . we can obtain similar — estimates for all intermediate 
p and q. 

To be able to derive time decay in the case of multiple roots, we will group terms 
in (12. 4p in the following way. Assume that roots Ti{^) , . . . , tl{C,) coincide on a set 
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contained in some A4, that is M. D = ' ' ' = '''^(0} • e > 0, we define 

M' := e M" : dist(^, M) < e} . Choose e > so that these roots ti{^),..., tl{^) 
do not intersect with any of the other roots tl+i{C,), ■ ■ ■ ,Tm(0 in A^^- If different 
numbers of roots intersect in different sets, we can apply the following theorem to 
such sets one by one. We note that by the strict hyperbolicity is bounded. Here 
we will estimate the sum 



r ^ 



(2. 



Theorem 2.2. Let the sum ( [^.^[ ) be the solution to the Cauchy problem l{2.1]) . As- 
sume that roots ti{$^) , . . . , tl{^) coincide in a set contained in M. and do not intersect 
other roots in the set Ai'^ . Let x ^ C^{A4'^). Assume that there exists 6 > such 
that Imrfc(,^) > 6 for all ^ G M"^ and k = 1, . . . , L. 
Then for all t > we have 



did: 



L 
k=l 



< C{l+t) 



L-1, 



~5t 



LP 



where - 
p 



1, 1 < P < 2. 



Thus, if characteristic roots are separated from the real axis on the support of 
some X ^ (^^(M"), we can separate the solution (12.41) into groups of multiple roots 



for which the L^ — L'^ norms still decay exponentially as stated in Theorem 12.21 We 
also note that since Ai"^ is bounded, assumption (ii) of Theorem 12.11 is automatically 
satisfied and, therefore, it is omitted in the formulation of Theorem 12.21 Theorem 
12.21 will be proved in Section [7l2l 



2.2 Roots with non-degeneracies 

The following case that we consider is the one of roots satisfying certain non-degeneracy 
conditions. These may be conditions on the Hessian, convexity conditions, or simply 
the information on the index of the corresponding level surfaces. In this section we 
will give the corresponding statements. We always assume the stability condition 
(12. 2p but no longer assume that roots are separated from the real axis. 

First we state the result for phases with the non-degenerate Hessian. The behavior 
depends on critical points with Vr(^°) = and the behavior of the Hessian at such 
points. As usual, we say that the critical point C,^ is non-degenerate if the Hessian 
Hess t{^^) is non-degenerate. 

Theorem 2.3. Let U G MJ^ be a bounded open set, and let t : U ^ C be smooth and 
such that Imr(^) > for all ^ G U. Assume that there are some constants Cq and 
M such that I det Hessr(0| > Co(l + 1^1)"^ for all ^eU. Let x e ^i%(M") be such 
that X = outside U and let a G Siq{U). 
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Assume that r has only one non-degenerate critical point in U , and that U is 
sufficiently small. Then there is a constant C > independent of the position of U 
such that for all t > we have 



<cii + ty 



-(---) 



Li(Wi) 



(2.9) 



with 1 < p < 2, i 



1, iVp = f( 



For example, the case of the Klein-Gordon equation corresponds to M = n + 2 in 
this theorem. If we work with a fixed bounded set f/, the ||/|| jvp norm on the right 

hand side of 02.91) can be replaced by 1 1/| |lp- However, since we may also want to have 
estimate fl2.9p uniform over such U (allowing it to move to infinity while remaining 
to be of the same size), we have the Sobolev norm in (12. 9p . From this point of view, 
we assume that a behaves as a symbol in f/ - the meaning is that if the symbolic 
constants here are uniform over the position of f/, then also the constant in (12. 9p is 
uniform over such a and U . 

The condition that critical points are isolated and therefore can be localised by 
different sets U may follow from certain properties of r and will be discussed in 
Section [6751 in particular see Lemma 16771 and remarks after it. If, in addition, we take 
the size of U uniform, say of volume bounded by one, then constant C in (12. 9p is 
also uniform over all such sets U . We may also assume that if is a critical point 
of r, then Imr(^°) = 0. Otherwise we would have Imr(^°) > and so Theorem 
12.11 would actually give the exponential decay rate. The proof of this theorem is 

If we apply 



based on the stationary phase method and will be given in Section 16.51 
different versions of the stationary phase method under different conditions, we can 
reach different conclusions here. For example, we also have: 

Theorem 2.4. Let U C he a hounded open and let t : U ^ C he smooth and 
such that Imr(,^) > for all ^ & U. Let x £ *S'5'g(]R") he such that x = outside U 
and let a G Siq{U). Assume that r has only one critical point in U , and that U 
is sufficiently small. 

Suppose that there are constants Co, M > independent of the size and position of 
U and of$,^, with the following conditions. Suppose i/iai rank Hess r(^°) = k, that this 
rank is attained on an kxk suhmatrix A(^°) and that \ det A(^°)| > Co(l + |^°|)~^^ 
Then for all t > we have 



<C(l+t) 



" 2 p 



with 1 <p<2, ^ 



) -z^- 



The proof of this theorem is similar to the proof of Theorem 12.31 once we restrict 
to the set of k variables (possibly after a suitable change) on which the rank of the 
Hessian is attained on A(^°). 
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This result can be improved dependent on further properties of For exam- 

ple, if rank/l(^'^) = n — 1 and this is attained on variables ^i, . . . , the analysis 
reduces to the behaviour of the oscillatory integral with respect to C,n- If the l-th 
derivative of the phase with respect to C,n is non-zero, we get an additional decay by 
This follows from the stationary phase method, see, for example Hormander 
|Hor83al Section 7.7], or from an appropriate use of van der Corput lemma. We will 
not formulate further statements here since they are quite straightforward. 

The next theorem is an estimate of oscillatory integrals with real-valued phases 
under convexity condition. It will be shown in Proposition 13.81 (see also Proposi- 
tion I6.16P that for large ^ characteristic roots of the Cauchy problem (12. ip satisfy 
assumptions of these theorems given below, if the homogeneous roots of the principal 
part satisfy them. The convexity condition is weaker than (but does not contain) the 
condition that the Hessian of r is positive definite and the result can be compared 
with Theorem 12.31 dependent on suitable properties of roots. 

Let us first give the necessary definitions. Given a smooth function r : — > M 
and A G M, set 

= S,(r) := e M" : r(0 = A} . 

In the case where t{C,) is homogeneous of order 1 and r G C°°(M"\0), we will also 
write St- := Ei(r) — for such r, we then have ^\{t) = AS,-. There should be no 
confusion in this notation since we always reserve letters 0, r for phases and A for the 
real number. 

Definition 2.5. A smooth function r : M" ^ M zs said to satisfy the convexity 
condition if surface is convex for each A G M. Note that the empty set and the 
point set are considered to be convex. 

If the Gaussian curvatures of never vanish, is automatically convex (the 
converse is not true). This curvature condition corresponds to the case k = n — 1 in 
Theorem 12. 4[ Another important notion is that of the maximal order of contact of a 
hypersurface, similar to the one in Section [L2} 

Definition 2.6. Let H he a hypersurface in M" {i.e. a manifold of dimension n — 1); 
let (7 G S, and denote the tangent plane at o by T„. Now let P be a 2-dimensional 
plane containing the normal to at a and denote the order of the contact between 
the line T^j (1 P and the curve T^DP by 7(S; cr, P). Then set 

7(S) := supsup7(S;cr,P) . 

o-es p 

Examples 2.7. 

(a) 7(S") = 2, as 7(S"; a,P) = 2 for all a G and all planes P containing cr and 
the origin. 

(b) If Lpi{^) is a characteristic root of an m*'^ order homogeneous strictly hyperbolic 
constant coefficient operator, then 7(S^J < m; see lSug96^ for a proof of this. 
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Now we can formulate the corresponding theorem. 



Theorem 2.8. Suppose r : M" M satisfies the convexity condition and let x ^ 
C°°(M"') ; furthermore, on suppx? we. assume: 

(i) for all multi-indices a there exists a constant Cq > such that 



|5^(0i<c«(i + iei) 



(ii) there exist constants M,C > such that for all \C,\ > M we have \t{^)\ > C\^\] 

(hi) there exists a constant C > such that \d^r{\uj)\ > C for all u G A > 0; 

m particular, |Vr(0| > C for all ^ G M" \ {0}; 

(iv) there exists a constant /2i > such that, for all X > 0, 

is,(r) = i{eeM":r(O = A}c5^,(0). 

Also, set 7 := sup^^^g 7(^a(t)) and assume this is finite. Let aj = %({) G S^q be a 
symbol of order —j of type (1, 0) on M". Then for all t > we have the estimate 



L9(Eg) 



< c{i + ty 



7 \P 1 ) 



(2.10) 



where 1 + 1 = 1, 1 < p < 2, and the Sobolev order satisfies Npj^t > ~ q) ~ 3 f'^''" 
0<t<l, andN,^,^t> (^-^) (J-^)-jMt> 1. 



Theorem 12.81 wiU be proved in Section 16.61 where estimate (12.101) wiU fohow by 
interpolation from the — estimate combined with — L°° cases given in (16. lip 
for small t, and in (I6.16P for large t. See those estimates also for the case of p = 1 
in estimate (I2.10p . The estimate for large times will follow from Theorem 14.81 which 
gives the L°°-estimate for the kernel of (12.101) . As another consequence of Theorem 
14.81 we will also have the following estimate: 



Corollary 2.9. Under conditions of Theorem \2.8\ with % = 1, assume that a G 
C^(M"). Then for all x G M" and t > we have the estimate 



a{Od^ <C{l+t) 



7 



(2.11) 



In Proposition 13.81 we show that properties (i)-(iv) of Theorem 12.81 are satisfied 
for characteristic roots of L{Dt, D^) in (12. ip . while in Lemma [6. Ill we will show that 
the index 7 is also finite, both for large frequencies. 

Now we turn to the case without convexity. As in the case of the homogeneous 
operators (see Introduction, Section II. 2p we introduce an analog of the order of 
contact also in the case where the convexity condition does not hold. 
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Definition 2.10. Let Ti be a hypersurface in M"; set 

7o(S) := sup inf 7(S; a, P) < 7(S), 



where 7(S; cr, P) is as in Definition \2.6[ 
Remark 2.11. 

(a) When n = 2, 70(E) = 7(E); 

(b) If p{C) is a polynomial of order m, E = G M" : p{C,) = 0} is compact, and 
Vp(^) on E, then 70(E) < 7(E) < m; this is useful when applying the result 
below to hyperbolic differential equations and is proved in jSugQff^ . 



Theorem 2.12. Suppose r : M" - 
furthermore, on supp x, we assume: 



is a smooth function. Let x G C"^ 



(i) for all multi-indices a there exist constants Ca > such that 

|9^M0l<Ca(i + iei)'-i"i; 

(ii) there exist constants M,C > such that for all \C,\ > M we have |t(^)| > C|^|; 

(ill) there exists a constant C > such that \dujT{Xuj)\ > C for all uo G S"""*^ and 
A > 0; 

(iv) there exists a constant Ri > such that, for all X > 0, 

Set 7o := sup;^>o 7o(Sa(t)) and assume it is finite. Let aj = aj{^) G be a symbol 
of order —j of type (1, 0) on M". Then for allt>0 we have the estimate 



<c{i + ty 



^(i_i) 

70 \ P 9 / 



>,j-t ) 



where ^ + ^ = 1, 1 < p < 2, and the Sobolev order satisfies Npj^t > n( 
0<t<l, andNp,,,t> (^n-^y^^-l)-j fort>l. 



1 _ l^ 

p q' 



j for 



The proof of Theorem 12. 121 wiU be given in Section I^TTl As in the convex case, as 
a consequence of estimates for the kernel on Theorem 15.31 we also have the following 
statement: 



Corollary 2.13. Under conditions of Theorem \2.1S\ with x = 1, assume that a G 
C^(M"'). Then for all x G M" and t > we have the estimate for the kernel: 



a{^)d^ < c{i + ty 



TO 
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Again, in Proposition 13.81 we show that properties (i)-(iv) of Theorem 12.121 are 
satisfied for characteristic roots of L{Dt,Drc) in fl2.1l) . while in Lemma [6.141 we will 
show that the index 70 is also finite, both for large frequencies. 

As a corollary and an example of these theorems, we get the following possibilities 
of decay for parts of solutions with roots on the axis. We can use a cut-off function 
X to microlocalise around points with different qualitative behaviour (hence we also 
do not have to worry about Sobolev orders). 

Corollary 2.14. Let Q C be an open set and let t : Q ^ W be a smooth real valued 
function. Let x ^ Let us make the following choices of K{t), depending on 

which of the following conditions are satisfied on suppx- 

(1) // det Hess r(0 ^ for all ^ e Q, we set K{t) = (1 + t)"^^^"?). 

(2) //rank Hess r(0 = n - 1 for all ^ E n, we set K{t) = (1 



n-l /I 1 - 



(3) If T satisfies the convexity condition with index 'y, we set K{t) = (l + t) ^ *-p 1'. 

(4) // T does not satisfy the convexity condition but has non-convex index 70, we set 
Kit) = (1 + 



Assume in each case that other assumptions of the corresponding Theorems \2.'d[^r[^ 
are satisfied. Let l<p<2, ^ + ^ = 1. Then for all t > we have 



<CK{t)\\f\\LnR-)- 



We note that no derivatives appear in the L^-norm of / because the support of 
X is bounded. In general, there are different ways to ensure the convexity condition 
for r. Thus, we can say that the principal part Lm of operator L{Dt, D^) in (12. ip 
satisfies the convexity condition if all Hessians <^'l{C)i ^ = 1, . . . ,m, are semi-definite 
for all ^ 7^ 0. In this case it was shown by Sugimoto in Sug94| that there exists a 



linear function a(^) such that = + a have convex level sets S(<^i), and we have 
7(E((^;)) < 2 [y] . For large frequencies, perturbation arguments imply that the same 
must be true for Y,\{ti), for sufficiently large A. If we now assume that Ti\{Ti) are 
also convex for small A, then ti will satisfy the convexity conditions. Alternatively, 
if they do not satisfy the convexity condition for small A, we can cut-off this regions 
and analyse the decay rates by other methods developed in this paper. 



2.3 Roots meeting the real axis 

In this section we will present the results for characteristic roots (or phase functions) 
in the upper complex plane near the real axis, that become real at some point or in 
some set. 

For M C W\ denote M" = G : dist(^, A^) < e} as before. The largest 
number G N such that meas(AI^) < Ce"^ for all sufficiently small e > 0, will be 
denoted by codim A^, and we will call it the codimension of M. 
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We will say that the root meets the real axis at with order Sk if Imrfc(^°) = 
and if there exists a constant Cq > such that 

Coie-eT'^- <lmrfc(0, 

for all ^ sufficiently near Here we may recall that in fl2.2l) we already assumed 
lmrfc(0 > for all ^. 

More generally, if the root meets the axis on the set Zk = E M" : Im Tk{^) = 0}, 
we will say that it meets the axis with order s if 

codist(e,^fe)' < lmrfc(0. 

We will localise around each connected component of Zk, e.g. around each point of 
Zk, if it is a union of isolated points. As usual, when we talk about multiple roots 
intersecting in a set Al, we adopt the terminology introduced in Section \n\ Since 
we are dealing with strictly hyperbolic equations, roots can meet each other only for 
bounded frequencies, so we may assume that set M. is bounded. 

Theorem 2.15. Assume that the characteristic roots ri(^), . . . , tl{^) intersect in the 
set Ai of codimension i. Assume also that they meet the real axis in M. with the 
finite orders < s, i.e. that 

Codisi{i,Mr <\mTk{0. 

for some cq > and all k = 1, . . . , L. Assume that ^2.4^ is the solution of the Cauchy 
problem h2. 1\) and we look at its part ^2.9\) . Let x ^ C^(7V1^) for sufficiently small 
£ > 0. Then for allt>0 we have 

L 
fe=l 

<C(l+t)-^(i-i)+^-^||/lUp, (2.12) 

where - + - = 1, 1 < p < 2. 

p ' q ' — ^ — 

We assume e > to be small enough to make sure that the type of behaviour 
assumed in the theorem is the only one that takes place in Ai^. In the complement of 
Ai^ we may use other theorems to analyse the decay rate. Moreover, we assume that 
set is C^. In fact, it is usually Lipschitz, so in order to avoid to go into depth about 
its structure and existence of almost everywhere different iable coordinate systems, we 
make the technical assumption. The proof of Theorem 12.151 will be given in Section 



Let us now give a special case of this theorem where simple roots meet the axis 
at a point, so that we have L = 1 and £ = n. The following statement is also global 
in frequency, so we have the result in Sobolev spaces. 
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Theorem 2.16. Consider the m^^ order strictly hyperbolic Cauchy problem 1(2. 1\] for 
operator L{Dt, D^) , with initial data fj G W^'''^^"'^'^^ forj = 0, . . . ,m — 1, where 
1 < p <2 and 2 < g < oo are such that ^ + ^ = 1,t>0 and a is a multi-index. We 
assume that the Sobolev index Np satisfies Np > — ■^) for 1 < p < 2 and Ni > n 
for p = 1. 

Assume that the characteristic roots Ti{^), . . . , Tm{C,) of L{t, C.) = satisfy Imrfc > 
for all k, and also the following conditions: 

(HI) for all k = 1, . . . ,m, we have 

lim inf Im (^) > ; 



(H2) for each G M" there is at most one index k for which lmTk{^^) = and there 
exists a constant c > such that 

ie-eT<cimr,(o, 

for ^ in some neighbourhood of . Assume also that there are finitely many 
points with lmTk{^^) = 0. 

Then the solution u = u{t, x) to Cauchy problem Ii2.1\) satisfies the following estimate 
for all t>0: 

m— 1 

WD^Dyit, OIU. < + ^||/,.||^^^,+,.,+._,. (2.13) 

j=0 



Theorem l2.16l is proved in Section [6. 11[ where we wiU also give microlocal versions 
of this resuh around points from hypothesis (H2). In the complement of such 
points, we have roots separated from the real axis, so we get the exponential decay 
from Theorems 12.11 and 12.21 Moreover, in the exponential decay zone we may have 
different versions of the estimate, for example we can use estimate (12. 7p there instead 
of (12.61) . As a special case, such estimate together with (12.151) below (used with 
s = Si = 2), we improve the indices in Sobolev spaces over for the dissipative 
wave equation in (11.71) and (II. 8p compared to |Mat77j . 

If conditions of Theorem 12.161 hold only with = 0, namely if Imrfc(^'^) = 
implies = 0, we will call the polynomial L(r, ^) strongly stable. Such polynomials 
will be discussed in more detail in applications in Section 18.51 Now we will give some 
improvements of (I2.13P under additional assumptions on the roots: 

Remark 2.17. The order of time decay in Theorem \2.16[ may be improved in the 
following cases, if we make additional assumptions. If, in addition, we assume that 
Imrfc(^°) = in (H2) implies that = 0, then we actually get the estimate 



Dinyit,-) ^ <c(i + t)-nl4J-^^||/^ 



m— 1 

Np + \c 
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where here and further in this remark Np is as in Theorem \2.1(k 

Now, assume further that for all in (H2) we also have the estimate 



\Tkm<ci\^-er 



(2.14) 



with some constant ci > 0, for all C, sufficiently close to 

If we have that Imrfe(^°) = in (H2) implies that we have ^2.14^ around such , 
then we actually get 



m—l 



DiDyit,-) ^ ^^<c{i + tri^)it^)~'^J2\\f: 



j=Q 



And finally, assume that for all such that Imrfc(^°) = m (H2), we also have 
= and Iji2.14\ l around such Then we actually get 



DiDy{t, 



L9(MS) 



< c{i + ty 



s \ p q I 



M_I£l 

s s 



m—l 



(2.15) 



Estimate fl2.15p with s = Si = 2 gives the decay estimate for the dissipative wave 
equation in (11.71) . The proof of this remark is given in Remark 16.191 

Moreover, there are other possibihties of multiple roots intersecting each other 
while lying entirely on the real axis. For example, this is the case for the wave equation 
or for more general equations with homogeneous symbols, when several roots meet 
at the origin. In this case roots always lie on the real axis, but they become irregular 
at the point of multiplicity, which is the origin for homogeneous roots. In the case 
when lower order terms are presents, characteristics roots are not homogeneous in 
general, so we can not eliminate time from the estimates as was done in Section II. 2[ 
It means that we have to look at the structure of such multiple points by making 
cut-offs around them and studying their structure in more detail. In particular, there 
is an interaction between low frequencies and large times, which does not take place 
for homogeneous symbols. The detailed discussion of this topic and corresponding 
decay rates will be determined in Section 17.41 



2.4 Application to the Cauchy problem 

Putting together theorems from previous sections we obtain the following conclusion 
about solutions to the Cauchy problem (12. ip . We will first formulate the following 
general result collecting statements of previous sections, and then will explain how 
this result can be used. 

Theorem 2.18. Suppose u = u{t,x) is the solution of the m}^ order linear., constant 
coefficient, strictly hyperbolic Cauchy problem 112. Denote the characteristic roots 
of the operator by ti{^) , . . . , Tm{0 > ^'^^ assume that Im rfc(,^) > for all k = 1, . . . ,n, 
and allie M". 

We introduce two functions, K^^^t) and K^^\t), which take values as follows: 
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I. Consider the behaviour of each characteristic root, rfc(^), in the region |^| > M, 
where M is a large enough real number. The following table gives values for the 
function K^\t) corresponding to possible properties of Tk{C}\ if Tk{C,) satisfies 
more than one, then take K^l^ (t) to be function that decays the slowest as t 
oo. 



Location of r^i^) 


Additional Property 




away from real axis 




e^"^*, some 6 > 


on real axis 


detHessrfc(0 
rank Hess T^iC.) = n — 1 
convexity condition 7 
no convexity condition, 70 


(l+t) g) 
(l+t) ~^P^^' 

(l+t) -y 1' 
(1 +t)~i 



Then take K^^\t) = ma^k=i ...,n Ki\t) . 

II. Consider the behaviour of the characteristic roots in the bounded region \^\ < 
M; again, takeK^^\t) to be the maximum {slowest decaying) function for which 
there are roots satisfying the conditions in the following table: 



Location of Root(s) 


Properties 


K'^'^^t) 


away from axis 


no multiplicities 
L roots coinciding 


e"*^*, some 6 > 
(1 + t)^e~'^* 


on axis, 
no multiplicities * 


detHessrfc(0 ^ 
convexity condition 7 
no convexity condition, 70 


(l+t) 1^ 

(l+t) ^ 

1(1 1\ 

(l+t) ^o^p 9' 


on axis, 
multiplicities* ,** 


L roots coincide 
on set of codimension £ 


(l+t)^-l-^ 


meeting axis 
with finite order s 


L roots coincide 
on set of codimension i 


(l+t) s'-p q' 



* These two cases of roots lying on the real axis require some additional regularity assump- 
tions; see corresponding microlocal statements for details. 

** This is the - rate in a shrinking region; see Proposition \ 7.y\ for details. For 
different types of L"^ estimates see Section \7.4l and then interpolate. 

Then, with K{t) = max {^K^^\t),K'^^\t)), the following estimate holds: 

m—l 

mOluit, < C.,rK{t) J] II// II ^^-P-^ , 

1=0 ^ 

where l<p<2, ^ + ^ = 1, and Np = Np{a, r) is a constant depending on p, a andr. 
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The scheme of the proof of this theorem and precise relations to microlocal the- 
orems of previous sections will be given in Section 12. 5[ However, let us now briefly 
explain how to understand this theorem. Since the decay rates do depend on the 
behaviour of characteristic roots in different regions and theorems from previous sec- 
tions determine the corresponding rates, in Theorem 12.181 we single out properties 
which determine the final decay rate. Since the same characteristic root, say t^, may 
exhibit different properties in different regions, we look at the corresponding rates 
K^^^{t), K^^^t) under each possible condition and then take the slowest one for the 
final answer. The value of the Sobolev index Np = Np{a,r) depends on the regions 
as well, and it can be found from microlocal statements of previous sections for each 
region. 

In conditions of Part I of the theorem, it can be shown by the perturbation 
arguments that only three cases are possible for large ^, namely, the characteristic 
root may be uniformly separated from the real axis, it may lie on the axis, or it may 
converge to the real axis at infinity. If, for example, the root lies on the axis and, in 
addition, it satisfies the convexity condition with index 7, we get the corresponding 
decay rate K^^\t) = (1 + t)~~^p~'^\ Indices 7 and 70 in the tables are defined 
as the maximum of the corresponding indices 7(Sa) and 7o(Sa), respectively, where 
12x = {C, ■ Tk{^) = X}, over all k and over all A, for which ^ lies in the corresponding 
region. At present, we do not have examples of characteristic roots tending to the 
real axis for large frequencies while remaining in the open upper half of the complex 
plane, so we do not give any estimates for this case in Theorem 12.181 However, in 
Section we will still discuss what happens in this case. 

The statement in Part II is more involved since we may have multiple roots 
intersecting on rather irregular sets. The number L of coinciding roots corresponds 
to the number of roots which actually contribute to the loss of regularity. For example, 
operator (d^ — A){df —2A) would have L = 2 for both pairs of roots ±1^1 and ±V2\^\, 
intersecting at the origin. Meeting the axis with finite order s means that we have 
the estimate 

dist(e,Zfe)^<c|Imrfc(OI (2.16) 

for all the intersecting roots, where Z^. = {C, : lmTk{C,) = 0}. In Part II of Theo- 
rem 12.181 the condition that L roots meet the axis with finite order s on a set of 
codimension i means that all these estimates hold and that there is a (C^) set Ai 
of codimension £ such that Zk C Ai for all corresponding k (see Theorem 12.151 for 
details). In Theorem 12.161 we discuss the special case of a single root meeting the 
axis at a point ^0 with order s, which means that Imrfc(^o) = and that we have the 
estimate \^ — (^o|^ < c|Imrfc(^)|. In fact, under certain conditions an improvement in 
this part of the estimates is possible, see Theorem 12. 161 and Remark 12.171 

In Part II of the theorem, condition ** is formulated in the region of the size 
decreasing with time: if we have L multiple roots which coincide on the real axis on 
a set A4 of codimension i, we have an estimate 

m— 1 

\u{t,x)\<C{l + t)'^-'-'J2\\fih^, (2.17) 



26 



if we cut off the Fourier transforms of the Cauchy data to the e-neighbourhood Ai'^ of 
Ai with e = 1/t. Here we may relax the definition of the intersection above and say 
that if L roots coincide in a set then they coincide on a set of codimension i if the 
measure of the e-neighborhood A4'^ of A4 satisfies < Ce^ for small e > 0; here 
= E M.'' : dist(^, A^) < e}. The estimate fl2.17p follows from the procedure 
described in Section 17.11 of the resolution of multiple roots, and details and proof of 
estimate fl2.17p are given in Section [731 especially in Proposition 17.91 

We can then combine this with the remaining cases outside of this neighborhood, 
where it is possible to establish decay by different arguments. In particular, this is 
the case of homogeneous equations with roots intersecting at the origin. However, 
one sometimes needs to introduce special norms to handle L^-estimates around the 
multiplicities. Details of this are given in the part of Section 17.4.21 in particular 
in Proposition 17. 5[ Finally, in the case of a simple root we may set L = 1, and i = n, 
if it meets the axis at a point. 

2.5 Schematic of method 

Let us briefly explain some ideas behind the reduction of Theorem 12.181 to the pro- 
ceeding theorems. The realisation of the steps below will be done in Sections M and 

13 

Step 1: Representation of the solution. 

Using the Fourier transform in x, this reduces the problem to studying time- dependent 
oscillatory integrals, at least for frequencies with no multiplicities. In the case near 
multiplicities we will introduce a special procedure to deal with them in Section [71 

Step 2: Division of the integral. 

We reduce the problem to several microlocal cases using suitable cut-off functions. 
The problem is divided into studying the behaviour of the characteristic roots in three 
regions of the phase space — large |.^|, bounded |.^| away from multiplicities of roots 
and bounded |.^| in a neighbourhood of multiplicities. 

Step 3: Interpolation reduces problem to finding — L°° and — estimates. 
Step 4: Large |,^|: 

• root separated from the real axis (Theorem 12. ip : 

• root lying on the real axis (Theorems I2.4f[2.12p . 

Step 5: Bounded |.^|, away from multiplicities: 

• root away from the real axis (Theorem 12. ip : 

• root meeting the real axis with finite order (Theorem 12. 16p : 

• root lying on the real axis (Theorems I2.4f[2.12l) . 
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Step 6: Bounded |^|, around multiplicities of roots: 

• all intersecting roots away from the real axis (Theorem I2.2p : 

• all intersecting roots lie on the real axis around the multiplicity (Section 

El; 

• all intersecting roots meet the real axis with finite order (Theorem 12.151) : 

• one or more of the roots meets the real axis with infinite order (similar 
to Theorems [23H2II2]). 

2.6 Strichartz estimates and nonlinear problems 

Let us denote by Kp^q{L{Dt, D^)) the time decay rate for the Cauchy problem (12.11) . 
so that function K(t) from Theorem 12. 181 satisfies K(t) ~ t~'^^''^^^^ for large t. Thus, 
for polynomial decay rates, we have 

«:,,,(L) = - hm — Ai. (2.18) 

i^oo In t 

We will also abbreviate the important case i^{L) = k,i ,^{L) since by interpolation 
we have = 1^2,2^ + '^i,oo(^ ~ 1 < P < 2. These indices k{L) and Kp^pi{L) 
of operator L{Dt, D^) will be responsible for the decay rate in the Strichartz esti- 
mates for solutions to (12.11) . and for the subsequent well-posedness properties of the 
corresponding semilinear equation which are discussed below. 

In order to present an application to nonlinear problems let us first consider the 
inhomogeneous equation 



L{Dt,D,)u = f, t>0, 
Dlu{0,x)=0, / = 0,...,m- 1, X G 



with L{Dt, Dx) as in (11.11) . By the Duhamel's formula the solution can be expressed 
as ^ 

uit) = [ E^_,{t-s)fis)ds, (2.20) 







where Em-i is given in (12. 4p . Let k = Kp^pi{L) be the time decay rate of operator L, 
determined by Theorem 12.181 and given in (I2.18p . Then Theorem 12.181 implies that 
we have estimate 

\\E„^-i{t)g\\w^,<C{l+ty''\\g\ 



\WS: 



Together with (12.201) this implies 



||w(t)lkMRS)<C / {t-sY^\\f{s)\\wsds<C\t\~^*\\fmwi. 
Jo 

By the Hardy-Littlewood-Sobolev theorem this is L^(M) — L'''(]R) bounded if 1 < 
g < 2 and 1 — k = ^ — ^. Therefore, this implies the following Strichartz estimate: 
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Theorem 2.19. Let Hp^pi he the time decay rate of the operator L{Dt,Dx) in the 
Cauchy problem 112. iy\) . Let 1 < p, q < 2 be such that 1/p+l/p' = 1/q+l/q' = 1 and 
1/q — l/q' = l — Kp^p'. Let s G M. Then there is a constant C such that the solution 
u to the Cauchy problem 112. 19\) satisfies 

\M\Li'(Rt,W^,iM.^)) - C'll/lli''(Kt,Wp^(R5))) 

for all data right hand side f = f{t,x). 

By the standard iteration method we obtain the weh-posedness resuh for the 
following semilinear equation 



L{Dt,D,)u = F{t,x,u), t>0, 
Dlu{0,x) = fi{x), / = 0,...,m-l, xG^" 



Theorem 2.20. Let Hpy be the time decay index of the operator L{Dt, D^) in the 
Cauchy problem h2.21\) . Let p,q be such that 1/p + 1/p' = 1/q + 1/q' = 1 and 
1/q - 1/q' = 1 - hip^pi. Let s G M. 

Assume that for any v E L'^ (Mj, W^i{W^)), the nonlinear term satisfies F{t, x, v) G 
L'^{Wt,Wp{M.^)). Moreover, assume that for every e > there exists a decomposition 
— oo = to < ti <■■■< tk = +00 such that the estimates 

\\F{t,x,u) - F{t,x,v)\\L'i{ij,w^{Rr^)) < e\\u - L<,'(7^.,m/^. (Rn)) 

hold for the intervals Ij = {tj, tj+i), j = 0, . . . , k — 1. 

Finally, assume that the solution of the corresponding homogeneous Cauchy prob- 
lem is in the space L^' {Rt, WpiR"^)). 

Then the semilinear Cauchy problem 112. 2 1\} has a unique solution in the space 



3 Properties of hyperbolic polynomials 

In order to study the solution u{t,x) to fll.ip . we must first know some properties 
of the characteristic roots ti{C,) , . . . , Tm{C,) ■ Naturally, we do not have explicit for- 
mulae for the roots, unlike in the cases of the dissipative wave equation and the 
Klein-Gordon equation (i.e. for second order equations), but we do know some 
properties for the roots of the principal symbol. For general hyperbolic operators, 
the roots ifi{C,), ■ ■ ■ , fmiO ^^e characteristic polynomial of the principal part are 
homogeneous functions of order 1 since the principal part is homogeneous. Further- 
more, for strictly hyperbolic polynomials these roots are distinct when ^ 7^ 0. Since 
these two properties are very useful when studying homogeneous (strictly) hyperbolic 
equations, it is useful to know whether the characteristic roots of the full equation, 
Ti{^), . . . , rm(0) have similar properties. Indeed, if we regard the full equation as a 
perturbation of the principal part by lower order terms, we can show that similar 
properties hold for large |.^|; these results are the focus of this section. In the out- 
line of the method in Section [531 we subdivided the phase space into large |^| and 
bounded |^|, and it is these properties that motivate this step. 
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3.1 General properties 

First, we give some properties of general polynomials which are useful to us. For 
constant coefficient polynomials, the following result holds: 

Lemma 3.1. Consider the polynomial over C with complex coefficients 

m 

Z"^ + CiZ"^-^ + ■■■ + Cm-lZ + Cm = \[{z - Zfc). 

fc=l 

// there exists M > such that \cj\ < for each j = 1, . . . ,m, then \zk\ < 2M for 
all k = 1, . . . ,m. 

Proof. Assume that \z\ > 2M. Then 

> (2M)"(1 - 2-^ 2-('"-i) - 2"™) > 0. 

That is, no zero of the polynomial lies outside of the ball about the origin of radius 
2M; hence \zk\ < 2M for each k = 1, . . . ,m. □ 

Remark 3.2. // we replace the hypothesis \cj\ < by \cj\ < M for each j = 
1, . . . , m, then by a similar argument we obtain that |2:fc| < max{2, 2M}. The quantity 
max{2, 2M} appears because we need M > 1 for the sum on the right hand side to 
be positive. 

For general polynomials with variable coefficients, we have continuous dependence 
of roots on coefficients (we give an independent proof of this result here for the sake 
of completeness and for referencing, but analogue of this result can be found in many 
monographs dealing with hyperbolic polynomials). 

Lemma 3.3. Consider the m^^ order polynomial with coefficients depending on & 
W 

p(r,0=r™ + ai(Or™-^ + --- + a„(0. 

// each of the coefficient functions aj{C,), j = 1, . . . ,m, is continuous in M" then each 
of the roots ri(^), . . . , rm(0 '"^^^^ respect to r of pij,^) = also continuous in R". 

Proof. Define p : C" — > C™ by p{zi, . . . , Zm) = (ci, . . . , c^) where the Cj satisfy 

m 

Z^ + Ciz"'-^ + --- + Cm = X[{z-Z,). 

i=i 

By the fundamental theorem of algebra p is invertible (but the inverse is not unique 
modulo permutation of roots), and, moreover, p is: 

(a) surjective by the Fundamental Theorem of Algebra; 
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(b) continuous since each of the Cj may be written as polynomials of the zj (by the 
Vieta formulae); 

(c) proper (that is, the preimage of each compact set is compact) by Remark I3.2( 



properties (b) and (c) imply that p is a closed mapping. 
Now, fix e M". For any given e > 0, consider the set 

m 

aSSm fc=l 

where a = {ai, . . . ,am) & Sm denotes the set of permutations of {1, ...,m} (see 
Fig. [1] for a diagram of this). Note that U is, by construction, symmetric, i.e. if 



2e 2e 



-2 , 












-0), 



























































Figure 1: U = UiUU2 

{zi, . . . , Zm) G U then (za^, . . . , Za^) G U for all (ai, . . . , am) G Sm- Let F denote the 
complement to U : 

P= n {C = {Ci,---Xm)eC^:\U-n{e)\>eBk = l,...,m}. 

We need to show that there exists 6 > such that (ri(,^), . . . , rm(0) ^ U whenever 
|^-{°| < 5; note: 

• P~^{pi.F)) = F hy construction — if p{w) = p{w') then both w and w' give rise to 
the same polynomial, and hence their entries are permutations of each other, and 
so either both or neither lie in F; 

• by the surjectivity of p, 

p{u) = p{F^) = p{[p-\p{F))r) = p{p-\p{Fr)) = p{Fr ■, 

• p{F) is closed since F a closed set and p is a closed mapping; 

therefore, p{U) is open. Thus, there exists an open ball in p{U) of radius S' (for some 
5' > 0) about aie) = («i(e°), • • • , amie)) = P(ri(e°), • • • , r„(e°)): 

Bs'{a{e)) = {(ci, . . . ,cj G C" : |c, - a,{e)\ < 5' Vj = 1, . . . , m} C p{U). 
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By the continuity of the aj{C,), there exists 6 > such that 

le - e°l < ^ - aj{e)\ < S' for an j = l, . . . , m ; 

hence, 

\^-^^\<5 ^ (ai(0,...,a^(0)e5y(a(e°))Cp(f/). 

Finally, since p(ri(^), . . . , rm(0) = {ciiiO y ■ ■ ■ ^ ^miO) f/ is symmetric (this is 
needed as different root orderings give the same coefficients), we find that we have 
(ri(^), 7:^(0) ^ U when |^ — < 5 as required; this completes the proof of the 
lemma. □ 

Now, let us turn to proving properties of the characteristic roots. 

Proposition 3.4. Let L = L{Dt, D^) he a linear m}^ order constant coefficient dif- 
ferential operator in Dt with coefficients that are pseudo- differential operators in x, 
with symbol 

m m 

where Pj{XO = ^^PjiO for all X >> 1, \^\ » 1, and aj G S^~^, for some e > 0. 

Then each of the characteristic roots of L, denoted ti(^), . . . ,Tm(^), is continuous 
in MJ^; furthermore, for each k = l,...,m, the characteristic root rfc(^) is smooth 
away from multiplicities, and analytic if the operator L{Dt, D^) is differential. 

If operator L{Dt,Dx) is strictly hyperbolic, then there exists a constant M such 
that, if \$,\ > M then the characteristic roots Ti(^), . . . , rm(0 of L are pairwise dis- 
tinct. 

Proof. The first part of Proposition is simple. Let us now investigate the structure 
of the characteristic determinant. We use the notation and results from Chapter 12 
of jGKZ94j concerning the discriminant Ap of the polynomial p{x) = PmX"^ + ■ ■ ■ + 
Pix + po, 

where the Xj (j = 1, . . . , m) are the roots of that is, the irreducible polynomial 
in the coefficients of the polynomial which vanishes when the polynomial has multiple 
roots. We note that Ap is a continuous function of the coefficients pq, . . . ,pm of p{x) 
and it is a homogeneous function of degree 2m — 2 in them; in addition, it satisfies 
the quasi-homogeneity property: 

A(po, Xpi, . . . , A>^) = A-(— 1) A(po, . . .,Prn). 

Furthermore, Ap = if and only if p(a;) has a double root. 
We write L{t,^) in the form 
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where 

m 

i=i 

is the principal part of L{t,^); note that the Pj{^) are homogeneous of degree j and 
the aj{^) are symbols of degree < j. By the homogeneity and quasi- homogeneity 
properties of A^,, we have, for A 7^ 0, 

Al{\0 = A(P,„(AO + a„(AO, . . . , Pi(AO + ai(AO, 1) 
= A(A™[P„,(0 + . . . , A[Pi(0 + 1) 

= A™(2™-2)A(p^(e) + . . . , A-('"-^)[Pi(0 + A-™) 

(using that A is homogenous of degree 2m — 2) 

= A"'('"-^)A(P^(0 + . . . , Pi(0 + 1) 

(by quasi-homogeneity) . 

Now, since L is strictly hyperbolic, the characteristic roots (pi{^), ■ ■ ■ ,^Pm{0 of Lm 
are pairwise distinct for ^ 7^ 0, so 

AlJO = A(P„.(0, • • • , Pi(0, 1) ^ for e ^ 0. 

Since the discriminant is continuous in each argument, there exists 6 > such that 
if I ^^^7^ I < 6 for all j = 1, . . . , m then 

|A(P^(0 + ^,...,Pi(0 + ^,l)MO, 

and hence the roots of the associated polynomial are pairwise distinct. So, fix ,^ G 
{(^ e : |,^| = 1} and let X —>■ 00. Since the aj(^) are polynomials of degree < j it 
follows that when |^| > M, the characteristic roots of L are pairwise distinct. □ 



3.2 Symbolic properties 

In this section we will establish a number of useful properties of characteristic roots 
which will be important for the subsequent analysis. In particular, we will show that 
asymptotically roots behave like symbols, and we will show the relation between roots 
of the full symbol of a strictly hyperbolic operator with homogeneous roots of the 
principal part. 

Proposition 3.5 (Symbolic properties of roots). Let L = L{Dt,Dx) be a hyperbolic 
operator of the following form 

m m 
i=l i=l \a\+m-j=K 

where P,(AO = A^P,(0 for A » I, \^\ » I, and c^j E Here < K < m - 1 
is the maximum order of the lower order terms of L. Let Ti {(,),..., Tm{0 denote its 
characteristic roots; then 
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I. for each k = 1, . . . , m, there exists a constant C > such that 

|rfc(OI <C(l + |e|) for all ^ER^. 

Furthermore, if we insist that L is strictly hyperbolic, and denote the roots of the 
principal part Lmij, C,) by ipi{C,) , . . . , ^m{C)i then we have the following properties as 
well: 

II. For each Tk{C,), k = 1, . . . ,m, there exists a corresponding root of the principal 
symbol LpkiO {possibly after reordering) such that 

\rkiO-MO\<Cil + \^\f^'-"' for all ^eR^. (3.1) 

In particular, for arbitrary lower terms, we have 

\rk{0-MO\<C for all ^eR"". (3.2) 

III. There exists M > such that, for each characteristic root of L and for each 
multi-index a, we can find constants C = Ck,a > such that 

|afrfc(0| <C7|er^"l forall\^\>M, (3.3) 

In particular, there exists a constant C > such that 

|Vrfc(0|<C forall\^\>M. (3.4) 

IV. There exists M > such that, for each Tk{C,) a corresponding root of the prin- 
cipal symbol (pk{0 be found {possibly after reordering) which satisfies, for 
each multi-index a and k = 1, . . . ,m, 

\d^TkiO - d^MO\ < C|er+'-™-l"l for all \^\ > M (3.5) 

In particular, since K < m — 1, we have 

\&lTk{i)-d''^^k{i)\<C\i\-\-\ forall\i\>M, (3.6) 

for each multi-index a and k = 1, . . . ,m. 

First, we need the following lemma about perturbation properties of general 
smooth functions. Clearly, we do not need to require that functions are smooth, 
but this will be the case in our application. 

Lemma 3.6. Let p,q : C ^ C be smooth functions and suppose zq is a simple 
zero of p{z) {i.e. p{zq) = 0, p'{zo) ^ 0). Consider, for each e > 0, the following 
^^perturbation" of p{z): 

Pe{z) := p{z) + eq{z) , 
and suppose z^ is a root of ps{z); then, for all sufficiently small e > 0, we have 

\z,-zo\<Ce4r\ ■ (3-7) 
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Proof. By Taylor's theorem, we have, near zo, 



Peiz) = Peizo) +p'^{zo){z - Zq) + 0{\z - Zol"^) 

= eq{zo) + {p'{zo) + eq'{xo)){z - zq) + 0{\z - 2:0^) 
Thus, setting z = Zs, we get 

= eq{zo) + {p'{zo) + eq'{zo)){ze - zo) + 0{\ze - zoH • 



(3. 



Now, consider the function of e, z{e) := Ze] this is clearly smooth since p and q are 
smooth and Zq is a simple zero of Indeed, p'^^iz^) ^ p'{zo) ^ for small e, hence 
Zs is a simple root of p^- Thus, near the origin. 



z{e) = z{0)+ez'{0) + Oie'^). 
Combining fl3.8p and fl3.9l) . we get 

= eq{zo) + (p'izo) + eq'izo))iez'iO) + 0{e^)) + 0{e^) , 

or, 

= q{zo)+p'izo)z'iO) + Oie), 
for small e. Therefore, by the triangle inequality, for each e > small enough. 



(3.9) 



^'(0)1 < 



Ce 



|p'(zo)| P'{zo) 



and, thus, 



z'(0)\ < C 



Finally, combining (13.101) with (13. 9p . we obtain (13.71) as required. 
Proof of Proposition \3.5[ 



(3.10) 
□ 



IPart II We may write L(r, ^) in the form 

L(r, 0=r'" + diiOr'''-' + ■■■ + a^_i(e)r + a^O, 
where |aj(^)| < C{C,)^. Hence for all k we have |Tfc(^)| < C{C,) by Lemma [3.11 

IPart lit In the proof of this part, let us write L{t,C,) in the form 

B. 

2=0 
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where rQ = 0, m ~ ri = K (the maximum order of the lower order terms), 1 < ri < 
■ ■ ■ < Tfi < m, 

m 

and Lm-rX'^,i) = ^ c„j(Ot^ ioi 1 < i < R; 

\a\+j=m-ri 

here, as usual, the Pj{$,) are homogeneous in of order j. 
Denote the roots of 

I 

1=0 

with respect to r by (,^), . . . , r^(^). Note that Lq^t,^) = Lm{T,^), i.e. Lq^t,^) is 
the principal symbol with no lower order terms. Since Iji{t, ^) are strictly hyperbolic, 
we will look at |^| > Mq, where all r{(^), . . . , r^(^) are distinct, for all /. 

We shall show that there exists M > Mq so that, possibly after reordering the 
roots, for all /c = 1, . . . , m, 

- riiOl < C\^n^^+' for all / = 0, . . . , /2 - 1 and |e| > M . (3.11) 

Assuming this, and noting that r°(^) = (pk{0 t^(0 = ''"fe(0 ^ach k = 1, . . . ,m 
(possibly after reordering), we obtain 

IniO - M0\ < (0 - riiOl < when \^\ > M; 

1=0 

this, together with the continuity of the Tfc(^) and ^PkiO — thus the boundedness 
of \TkiO - M0\ in Bm{0), gives ([31]). Then, ([321) follows by setting K = m - 1. 
Here we also used ri = m — K. 

So, with the aim of proving (13.111) . we first introduce some notation: set 

Lm_r^ : C X ^ C : Zm_r^(r,u;) = Lm_r,(r,cj) , z = 0, . . . , i?, 
: (Mo,oo) X C X §"-1 ^ C : L;(p, r, cu) = p"'"Li(pr, pcj), / = 0, . . . , i?; 

observe that Lm-n is just the restriction of Lm-rXTiC) to C x Denote by 

^i(cij), (t^), • • • , ^miyj) the roots of Lmij, uj) = Lo(p, r, cu) with respect to r, and by 
7f(p,t^),7t(p,^), • • -^T^ip,^) those of Lfc(p,r,u;). 
We denote t = -^^ Since, 

L„,{t,^)=L„,{t,^) = |e|-'"L^(r,0 = |e|-'"Lo(r,0 = Lo(|e|,r,|) 
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for ^ e M", r G C, and 

U+i (p, r,uj) = p~'"Lz+i (pr, puj) = p^™ (pr, pu;) 



i+i 



i=0 



i=0 



|a|+i=m-n+i 

=Li(p, r, uj) + p~'''+'Ll^_^^^^ (p, r, u) 



|a|+i=m-rj+i 
Ca,j{puj) 



\a\ 



for e S''-^ p>Mo,TeC,l = 0,...,R-l. Here 



CajipUj) 



\a\+j=m-ri+i 



P 



\a\ 



We also have ^ 

ler"^L^(r,0 = L/(|e|,|,r). 

As the left-hand side of this is zero when t = r^C,), k = 1, . . . , m, and the right-hand 
side is zero when r = t|(|^|, |||), = 1, . . . , m, we see that l^lr^d^l, ||) = r^(0 for 
each k = 1, . . . ,m (possibly after reordering). Hence, for all |.^| > Mq, k = 1, . . . ,m 
and / = 0, . . . , -R — 1, we have 

\ri-'\0-rU0\ = \rl:''m4)-rim4)m. 
Next, observe that applying Lemma [3.61 with e = p~^'+'^ to 

Li{p, T,uj)+p~^^+^Ll_^.^^^ip,T,u) 



yields, for all G and k = 1, . . . ,m, 

\Ti+\p,Uj)-TUp,u)\<Cp-^-'+^ 



drU{p,rl{p,u;),uj) 



provided we take p > M' for a sufficiently large constant M' > Mq. Therefore, for 
all |.^| > M', k = 1, . . . ,m and / = 0, . . . , i? — 1, we have 



rl-^'iO - rim < C\^\ 



-n+i+i 



-n+i 



drU 



I) 



(3,12) 



Thus, it suffices to show the following two inequalities when |^| > M for some M > 
M': 
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• there exists a constant Ci so that, for all 1 < i < R, 



7-0 ZiM ±\ 



E 



\a\+j=m-ri 



(3.13) 



and 



there exists a constant C2 > so that, for all < / < i? — 1, 



)\ = \^\-+'\dMrU0,0\>C2. 



(3.14) 



Then, combining fl332|) . fl3:T3D and fl3J4|l gives fl3TTD . 

The first estimate (13.131) follows immediately from lPart II since the t^(0 ^ire roots 
of strictly hyperbolic equations, and from the fact that Caj G 5*'"'. 

The second, (I3.14p . in the case Z = is clear: the homogeneity of Lm{T,^) and its 
roots give 

which is never zero due to the strict hyperbolicity of Lm and hence (using that the 
sphere S"'~^ is compact and L^^r,^) is continuous and thus achieves its minimum) 
is bounded below by some positive constant as required. 

For 1 < / < i? — 1, we know that r^(,^), /c = 1, . . . , m, are simple zeros of Li(r, ^) 
for 1^1 > Mo by the earlier choice of Mq. Observe, 

{drUMiO.i) {drLmWm.O , (^.^rn-n ) (t^ (0 , 



Now, 



|^|,n-l 
1^1^-1 



i=l 



as 1^1 — s> oo 



for z = 1, . . . , /, because drLm-niT^O is a symbol of order m — — 1. Also, using the 
Mean Value Theorem, 

(9.L,„)(r^(0,0 = {drL^){^u{0.0 + RL„)(r^(0,0 - {drL^){^u{0.0] 

= {drLm){MO,0 + (9^Lj(r^(0,0 , 

where r^(0 li^s on the line connecting (pkiO '^liO ^ach ^G]R",/c = l,...,m 
and / = 1, . . . , -R — 1, and 



|(^^^r.)(f^(0,0| 

|^|,n-l 



<cier'^oas lei^oo. 



Therefore, for a sufficiently large constant M > M', there exists a constant C2 > 
such that 



This completes the proof of (I3.13p and thus of lPart III 
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IPart Hit We take M > so that for |^| > M, the roots ri(0, • • • , r^(0 are 
distinct. 

To prove the statement, we do induction on \a\. 

First, assume |a| = 1. Since L{Tk{C,),0 — each = 1, . . . , m, we have, for 



each i 



n, 



The first term is a symbol of order m — 1 in (rfc(^), hence, by IPart T\ there exists 
a constant C such that, when |.^| > Mi for some suitably large constant Mi > M, 

The inequality (13.31) for |a| = 1 (i.e. (13.41) ) then follows immediately from: 

Lemma 3.7. There exists constants C > 0, M2 > M such that, for each k = 
1, . . . ,m, 

dL 



Proof. Note that 



> 



dr 

dim 



m—1 



when 1^1 > M2 



(3.15) 



where Lm{T,C,) is the principal symbol of L and (pi{^), . . . , (pmiO are the correspond- 
ing characteristic roots, ordered in the same way as in IPart III We look at each of 
the terms on the right-hand side in turn: 

• By strict hyperbolicity, ^§^{'^k{0^0 non-zero for 7^ 0. Thus, for all 7^ 0, 



dir. 



= i^r"i^(|'^(|))| > (3-16) 



Observe 
dL 



dr 



dr 



dLjy, 
~d7 



Now, 
dLr 



dr 



= m{rk{Cl 



dL„ 
~d^ 



m—l 



dL„ 
~d^ 



m—l 



r=0 \a\+l=r 



l-l 



J — 1 



39 



and 



So, by IPart Tl and [Part III (specifically inequality fl3.2p ) and the fact that the -Pj(^) 
are homogeneous in ^ of order j, we have, for some suitably large M2 > M, 

-Q^irk{0,0-^{MO,0 <C\m-' when|e|>M2. 
This, together with 

I Yl lCaAOrk{0'''\<C\^r'<C\^r-' when|e|>M2, r = 0,...,m-l, 

\a\+l=r 



which again follows straight from IPart II yields 
fir r)T 



(3.17) 



The result now follows by combining flS.lSp , fl3.17p and fl3.16p . The proof of Lemma 
13.71 is complete. □ 

For |a| = J > 1, assume inductively that, 

|95"rfc(0| <C|er^"l when|e|>M, |a| < J - 1 , 

for some fixed M > max(Mi, M2). 

Then, for |a| = J, we use d^[L{Tk{^),^)] = 0, i.e. 



Ca,/3i,...,/3^ 



/3l + ---+/5'-<a, 



By the inductive hypothesis and the fact that d^dlL{Tk{Cj,i) is a symbol of order 
m— j — we have, for all multi-indices [3^, . . . , /J*" 7^ or a satisfying [3^+- ■ < a, 



< Cfc,„|er"l"l when 1^1 > M. 



Thus, using Lemma 13.71 again, we have 



\drLiniO,0\ 

which completes the proof of the induction step. 



< Cfc,,|er-I"l when lel > M, 
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IPart IVt Once again, assume that the roots Tk{^), k = l,...,m, correspond to 
fkiO^ k = 1, ■ ■ ■ ,m, in the manner of IPart III 

The proof of this part for general multi-index a is quite technical, so we first give 
the proof in the case |a| = 1 to demonstrate the main ideas required, and then show 
how it can be extended when |a| > 1. 

From L{Tk{^),C,) = = L^i^PkiO^Oy have for each i = 1, . . . ,n, 

|(^.(0,0 + f(r.(0,o|^(0-o. 



dr 



0. 



Therefore, 
dL 



^^"■(^.(0,0-%^(r.(0,0 



+ 



diPk 



^(^.(0,0-^(^.(0,0 



(r,(0,0- (3.18) 



It suffices to show that the right-hand side is bounded absolutely by C|^|"*~^ when 
1^1 > Ml for some suitably large Mi > Mq; this is because an application of 
Lemma [3.71 then yields 



— (0 - ^(0 



< 



CIO 



m-2 



'f (r.(0,0| 



<C|^|-^ for 1^1 >M^ 



where M = max (Mi, Mg). 

Since — Lm)(T, ^) is a symbol of order < m — 2 in (r, ^), it is immediately clear 
that the final term of (13.181) is bounded by C|^|™~^; here we have also used IPart I[ 
Also, noting that |%v^/t(OI < by the homogeneity of v^fc(0, we have, by (I3.17p . 



(0 



dLr^ 



dr 



■(^.(0,0 



Finally, by the Mean Value Theorem, 



dir. 



■(^.(0,0 



dLr, 



(r.(0,0 



~d7 



< C 



(rfc(0,0 <C|0 



m-2 



drdC 



where f lies on the linear path between (pk{0 Tfc(^) — which means that (using 
IPart II once more) |f| < C\^\ for |,^| > M. Since drd^^Lm{T,^) is a symbol of order 
m — 2 in (r, ^), and \ipk{C,) — Tfc(OI ^ bv IPart Ilf this term is bounded by C\C,\'^~'^, 
completing the proof in the case |a| = 1. 

For |a| = J > 1, we assume inductively that 



\d?Tk{0-d?MO\ < C|0~'"' for 10 > M,\a\ < J -I 
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As in the proof of IPart ITTI we have 

r 

/3lH h/3'"<a, J = l 

similarly, 

r 

/3l + .-- + /3''<a, i = l 

Thus, 

r 

/9l + --- + /3''<a, J = l 

r 

f3T- + ...+f3r<a, j = l 

r 

/3^ + -- + ISr<a, j = l 

We claim the right-hand side is then bounded absolutely by Cq|.^|™~^~'"', which, 
together with Lemma [3 .7^ yields the desired estimate. 
To see this, let us look at each of the terms in turn: 

• Id^fkiOl — C*Q,|.C|^~'"' by the homogeneity of </?fc(0; using this with fl3.17p gives 
the desired bound. 

• Using the Mean Value Theorem as in the case |a| = 1, we get 

coupled with \d^MO\ < C'a^r"''^', this gives the correct bound. 

42 



• By the inductive hypothesis, 
together with 

which follows from IPart II and the homogeneity of Lm(r, this gives the correct 
estimate. 

• To show the final term is bounded absolutely by first note that 

is a symbol of order < m — \a\ + \f]^\ + ■ ■ ■ + \(3^\ — r — 1; applying IPart IIll to 
estimate the d^\kiO terms, we have the required result. 

This completes the proof of (13. 6p : (13.51) is proved in a similar way in the proof using 
the set-up of the proof of IPart iTl The proof of Proposition 13.51 is now complete. □ 
We will now establish further symbolic properties of characteristic roots. A re- 
finement of this proposition concerning real and imaginary parts of complex roots r 
is given in Proposition 16.161 

Proposition 3.8. Suppose that the characteristic roots (pk, k = l,...,m, of the 

principal part Lm{T,C,) of a strictly hyperbolic operator L{t,^) in Ii2.1\) are non-zero 
for all C, 7^ 0. Then the roots r(^) of the full symbols satisfy the following properties: 

(i) for all multi-indices a there exists a constants M, Ca > such that 

|95V(OI<c«ier-'"'; 

for all 1^1 > M. 

(ii) there exist constants M,C > such that for all \^\ > M we have |t(,^)| > C\^\] 

(iii) there exists a constant Cq > such that \d^T{Xuj)\ > Cq for all uo G S"~^, 
A > 0; m particular, | Vr(OI > Cq for all ^ e M" \ {0}; 

(iv) there exists a constant Ri > such that, for all A > 0, 

is,(r) = G : r(0 = A} C 5^,(0) . 



Proof. • Property (i): by Proposition l331IPart IIH 

<Ca|er-'"' forall|e|>M^ 

for all multi-indices a. 
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Properties (ii) and (iii) : these follow by using perturbation methods. By Propo- 
sition [HSl [ParOYl there exists a homogeneous function ip{^) of order 1 such that, 
for all 1^1 > M and k = l,...,n, 

|r(0 - ^(01 < Co and |%r(0 - d^MOl < Ck\^r' , 

for some constants Co,Ck > 0. Now, the homogeneity of (p{^) implies that ip{^) = 
and Cfc ■ V(p{ek) = (p{ek), where = ( 0, ■ ■ . , 0, 1 , 0, . . . , 0), so 

k 

M0\ > C'lel for all e e M" and |9^</?(Acu)| > C for all u G §""\ A > , 
for some constant C > 0. Thus, 

|r(OI > 1^(01 - HO - viOl > C'\0 -Co> C\0 for |e| > M, (3.19) 
for some constants M, C > 0, and 

\d^r{Xuj)\ > \dMXuj)\ - IdM^uj) - d^T{Xuj)\ >C'- CkX-' > C > 

for all uj G S"~^ and suitably large A; for small A > 0, di^r^Xu) is separated from 
by the convexity condition, so \di^T{Xuj)\ > C > for all G S"^^, A > 0, as 
required. 

Property (iv)] — there exists a constant -Ri > such that, for all A > 0, jT.x(t) C 
-Siii(O) — holds by Proposition 13. 5[ IPart III and the fact that jJ^xif) = for 
the characteristic root of the principal symbol (p corresponding to r. 

□ 



4 Oscillatory integrals with convexity 

As discussed in Section 11.21 in the case of homogeneous m}^ order strictly hyperbolic 
operators, geometric properties of the characteristic roots play the fundamental role 
in determining the — L'^ decay; in particular, if the characteristic roots satisfy 
the convexity condition of Definition 11.11 then the decay is, in general, more rapid 
than when they do not. We will show that a similar improvement can be obtained 
for operators with lower order terms when a suitable 'convexity condition' holds. 
In Section 14.31 we shall extend this notion of the convexity condition to functions 
r : — > M and prove a decay estimate for an oscillatory integral (related to the 
solution representation for a strictly hyperbolic operator) with phase function r. 

First, we give a general result for oscillatory integrals and show how the concept 
of functions of "convex type" allows its apphcation to derive the time decay. 
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4.1 Estimates for oscillatory integrals 

The following theorem is central in proving results involving convexity conditions. In 
some sense, it bridges the gap between the man der Corput Lemma and the method 
of stationary phase, in that the former is used when there is no convexity but gives 
a weaker result, while the latter can be used when a stronger condition than simply 
convexity holds and gives a better result. Here, we state and prove a result that has 
no reference to convexity; however, in the following section, we show how convexity 
(in some sense) enables this result to be used in applications. An earlier version of this 
result has appeared in [Ruzh07j . with applications to equations with time dependent 
homogeneous symbols in [MR07j . For completeness we also include a more detailed 
proof here. 

Theorem 4.1. Consider the oscillatory integral 

I{X,u)= [ e^'''^y'^^Aiy,u)giy)dy, (4.1) 



where G N, / : [0, oo) x A/" — C, A/" any set of parameters v, and 

(11) there exists a hounded open set U C such that g G C^{U); 

(12) $(?/,//) is a complex-valued function such that Im $(?/,!/) > for all y E U, 
u eAf; 

(13) for some fixed z G M'^, some 5 > 0, and some 7 G N, 7 > 2, the function 

F{p, uj, u) := ^{pLU + z, u) 

satisfies 

\d,F{p,u,iy)\ > Cp'-^ and \d^F{p,u,u)\ < anp'-"'\d,F{p,u,u)\ 

for all {llJ, u) G E>^^^ x A/", all integers 1 < m < [N/'-f] + 1 and all p > 0, for 
which pu + z E U; 

(14) for each multi-index a such that |a| < [^] + 1, there exists a constant Ca > 
such that \dyA{y, < Ca for ally eU, u e M. 

Then there exists a constant C = Cn^-j > such that 

N 

\I{X,u)\<C{l + \y~ for all \e [0,00), u eAf. (4.2) 

Constant C in (14.21) is independent of A and u. 

Remark 4.2. This theorem extends to the case where A{y, u) is replaced by A{y, u'), 
where v' may he independent of the variable v appearing in the phase function $(?/, z/); 
these parameters do not have to he related in any way, provided the estimates in 



hypotheses (12) and (14) hold uniformly in the appropriate parameters. We will simply 



unite hoth sets of parameters and call this union v again. 
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Proof. It is clear that 04.21) holds for < A < 1 since |/(A, z/)| is bounded for such A. 

Now, consider the case where A > 1. Set y = pu + z, where u G S^~^ (using the 
convention that S'^ = { — 1, 1}), p > and z G is some fixed point; then 

I{\,iy)= / e'^'^^''^+''''^A{pu + z,iy)g{puj + z)p^-Updiu. 

By the compactness of it suffices to prove (14.21) for the inner integral. 

Choose a function x ^ C*o°([0, oo)), < x(s) < 1 for all s, which is identically 1 
on < s < I and is zero when s > 1; then, writing F{p,LJ,h') = ^{pu + z,^), we 
split the inner integral into the sum of the two integrals 



/l(A, U,LJ, z) 



e^'^^^'^^'-'^Aipu; + z, u)g{puj + z)x{\^ p)p'^~' dp, 
^^mP,^,^)A[puJ + z, v)g{pu + z){l- x)(A^p)p^-^ dp . 



Let us first look at Ii = /i(A, z/, u;, 2;); since x('^^p) is zero for Xi p > 1, we have, 

^ i 

by the change of variables p = Xi p, 

poo POO 

\h\ < C / x(A^p)p^-^ dp = C (p)^-iA-^x(p)A" dp 
Jo Jo 

N N 

<CX^^ / {p}^-Up = CX-^ , 
Jo 

where we have used \e^^^^^'^''^^\ < 1 since lmF{p,uj, u) > for all p,uj, v by hypothe- 
sis 



(12) this is the desired estimate for |/i 



In order to estimate I2 = hiX, u^ui, z), let us first define the operator L 
{iXdpF(p,cu, i'))'^^ and observe that 

Denoting the adjoint of L by L*, we have, for each / G N U {0}, 

e^^^^P^^^'\n\A{puj + z, u)g{puj + z){l - x)(A^p)p^-^] dp . 



Now, 



{d,Fr 



dp' 



where the sum is over all integers Si, . . . , Sp, p, r > such that Si 



Sp + r-p = I. 



By Hypothesis (13) 



d'p'F ...dp'F 



{dpFY+P 



(P, t^] 



< CpP~''' 



Cp 
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Also, we claim that, for r < [— ] + 1 



dp' 



■[A{puj + z, u)g{puj + z){l- x)(A^p)p^-'] 



N-l-r 



x{\p) 



(4.3) 



where x(A,p) is a smooth function in p which is zero for \^ p < |. Assuming this is 
true, we see that, for large enough I — it suffices to take / = [^J + 1, i-e. N — l'j < — we 
have, 



/•oo 



<CmX 



X 7 



CmX 



N-l-y 



_P 



Cn,"iX 



together with the estimate for |Ji|, this yields the desired estimate (14.21) . Here we 
need / > which means an application of {L*)\ or estimates on d'^F for |a| < /. 
This gives a restriction on the number m of derivatives in (13). 
Finally, let us check (14.31) . It holds because: 

(i) |9;(p^^i)| < Cr,Np^^^"' for all r G N. 

i r 1 i 

(ii) For each r G N, 9^[(1 — x)(A^p)] = — At (9^x)(A^p); now, {dsX){X'' p) is sup- 
ported on the set |(A,p) G (0, 00) x (0, 00) : | < \~i p < l|, so, in particular, 

on its support At < p^ ; therefore, 

|9:[(1-x)(Atp)]| <Cp-'-(9:x)(ATp) for all rGN, 



and {dlx){\i p) is smooth in p and zero for ATp < i. 



(iii) By hypothesis (14), \^'^A{puJ + z,u)\ < Cr for each r < [ 



1 (this can be 
and 



I y=puj+z 



seen for r = 1 by noting that dpA{puj + z,^) = lu ■ \/yA{y, u) 
then for r > 2 by calculating the higher derivatives). Also, g is smooth in 
f/, so, \d''p\A{puj + z,u)g{puj + z)]\ < Cr for r < [^] + 1- Furthermore, by 
hypothesis (II) , there exists a constant po > so that g{pu! + 2;) = for p > po; 
thus, dp[A{puj + 2, u)g{puj + z)] is zero for p > po; hence, for r <[^] + 1, 

mA{puj + z, v)g{puj + z)]\< CrPlp-"" . 



This completes the proof of the claim, and thus the theorem. 



□ 



4.2 Functions of convex type 



Hypothesis (13) of Theorem 14.11 is sufficient for the result of the theorem to hold; 
however, it is often difficult to check. For this reason, we now introduce the concept 
of a function of convex type — a condition that is far simpler to verify — and show that 
for such functions, (13) automatically holds. 
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Definition 4.3. Let F = F{p,v) : [0, oo) x T ^ C be a function that is smooth in p 
for each fixed t; G T, where T is some parameter space. Write its M*'* order Taylor 
expansion in p about in the form 



M 



F{p, v) = Y^ (^ji^)p' + Rm{p, v) , (4.4) 

j=0 

where Rm{p,v) = Jq 9g^~^^F{s,v) ^'^~^] ds is the M*'^ remainder term. 

We say that F is a function of convex type 7 z/, for some 7 G N, 7 > 2, and for 
some 5 > we have 

(CTl) ao{v) = ai{v) = for all v E T {i.e. the Taylor expansion of F starts from 
order > 2); 

(CT2) there exists a constant C > such that J2]=2\^ji'^)\ — ^ /'^'^ v E T; 

(CT3) for each v E T, \dpF{p, v)\ is increasing in p for < p < 6; 

(CT4) for each k E N, dpF{p, v) is bounded uniformly inO<p<S,vET. 



Remark 4.4. Note that, if F is real-valued, then (CT3) implies that we have either 
dpF{p, v) > for all < p < 5, or dpF{p, v) < for all < p < 5— this is 
because dpF[0, u) = 0. This is the connection with convexity, hence the name of such 
functions. 

Such functions have the following useful property: 

Lemma 4.5. Let F{p,v) be a function of convex type 7. Then, for each sufficiently 
small < 5 < 1 there exist constants C, Cm > such that 

\dpFip,v)\>Cp''-' (4.5) 
and \d^F{p,v)\ < Cmp'^"'\dpF{p,v)\ (4.6) 

for all < p < 6, V E T and m G N. 

Remark 4.6. A version of this lemma appeared in \Sug94^ for analytic functions 
without dependence on v and is based on Lemmas 3, 4 o-^^d 5 of Randol \Ran69^ 
(which also appeared in Reals \Rea8^ . Lemmas 3.2, 3.3). Lemma \4.5\ extends it to 
functions that are only smooth and which depend on an additional parameter, which 
will be necessary of our analysis. A limited regularity version of this lemma appeared 
in Ruzh07l The proof of lemma given here is based on estimating the remainder 



rather than on using the Cauchy's integral formula for analytic functions. 



Proof. First, let us note that, for < p < 1 we have, by (CT2) 



<P.v) := 5^ J>,(t;) 1/7^-1 > Cp^~' . (4.7) 

J=2 
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Thus, in order to prove (14.51) . it suffices to show 

\dpF{p,v)\ > C7r{p,v) for all < p < 5, t; G T; 
For 1 < m < 7, we have, using (14. 4p . 

7— m , , , , 

d'^F{p,v) = ^ —^ak+m{v)p'' + Rm,^-m{p,v) 



k=0 



where Rma-m{p, v) = v) 



ip-s) 



(4.8) 



(4.9) 



th 



m,-y-miP.v)\ < C,,„,.P < C,,^7r(p, t;)p2"™ for < p < 5 . 



Taylor expansion of d^F{p,v). By (CT4) and (14. 7p . we see 
\R 

Hence, for < p < 5 

7-1 



ds is the remainder term of the (7— m) 
md (14. 7p . we see 

(4.10) 



\dpF{p,v)\ = |^(/i; + l)afc+i(t;)p'' + i?i,^_i(p,t;)| 

fc=0 

I ^ 7 
^ ^ ^l,7-l(p.^) > Xl^^j(^)^''"^ -C^^ip,v)p- 



i=2 



i=2 



Now, by (CT3) , \dpF{p,v)\ is increasing in p for each t; G T and, by (CTl) 
dpF{0,v) = 0; therefore, 

\dpF{p,v) \ = max\dpF{a,v)\ 

0<(T<p 



> max > jaj(v)a-' 

~ 0<a<p\'^ ■'^ ' 
~ ~ i=2 



max C^7r{(j, v)a 

0<cr<p 



C^'K{p,v)p. 



since 7r((T, t')(j = ^J=2 j clearly achieves its maximum on < o" < p at 

a = p. Noting that 



0<CT<1 



max I > Zjff- 



and 



are norms on and, hence, are equivalent, we immediately get 

7 

\dpF{p,v)\ >CY,3\a^{v)\pi-'-C,^{p,v)p 
i=2 

> (C - C^(5)7r(p, v) = C^,5vr(p, w) , 
which completes the proof of (14. 8p . 
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To prove (14.61) . we consider the cases 1 < m < 7 and m > 'j separately. 



For m > 7, we have, by (CT4) 



KF{p,v)\ <Cm< C^.^p^+i-'" for < p < 5, 
since 7 + 1 — m < 0, and, thus, p^^^~^ > ^7+i-"i ^ g. t^q have 

\d;'F{p,v)\<C^,sp'~"'\dpF{p,v)\ for0<p<5, m>7. (4.11) 
For 1 < m < 7, we have the representation (14.91) . It is clear that 



fc=0 



which, together with (I4.10p and (14.81) . yields 

\d';F{p,v)\ < Cm,5p'-ndpF{p,v)\ for < p < (5, 1 < m < 7. 
This, together with (14. lip , completes the proof of (14. 6 p and, thus, the lemma. □ 



This lemma means we have the following alternative version of Theorem 14. 1[ 



Corollary 4.7. Hypothesis (13) of Theorem \4.1\ may be replaced by: 



(13') for some fixed z G M^, the function F{p,uj,v) := $(pu; + z,v) is a function 
of convex type 7, for some 7 G N, m the sense of Definition \4-3\ with (cj, v) G 
S^-i xM=T. 

4.3 Convexity condition for real-valued phase functions 

Using the results of the previous two sections, we can now prove a series of results 
for which a so-called convexity condition holds; here we recall Definitions 12.51 and 12.61 
from Section [2] and prove the basic result for real-valued functions. We recall that a 
smooth function r : M" ^ R is said to satisfy the convexity condition if Tj\ is convex 
for each A G M (and the empty set is considered to be convex). The maximal order of 
contact of a hypersurface S is defined as follows. Let ex G E, and denote the tangent 
plane at a by T^. Let P be a plane containing the normal to S at a and denote the 
order of the contact between the line T^DP and the curve SnP by 7(E; a, P). Then 
we set 

7(E) := supsup7(E;a,P) . 

o-gS P 

In the proof of Theorem 12.81 we will need a Besov space version of the estimate 
for the kernel. For this, let us introduce some useful notation for a family of cut-off 
functions Qr G C^(M'^), R G [0, 00): these functions will correspond to the cut-offs to 
annuli in the frequency space and we need to trace the dependence on the parameter 
R. Suppose g G C^(]R") is such that, for some constants co,Ci > 0, it is supported 
in the set 

: Co < < ci} , 
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and let go G C^(M" \ {0}) be another (arbitrary) compactly supported function. 
Then, for i? > 0, set 

[^o(0 if < i? < 1. 

Now we can prove the main convexity theorem: 

Theorem 4.8. Suppose r : M" — > M satisfies the convexity condition. Set 7 := 
s^Pa>o 7(^a(t)) and assume this is finite. Let a{^) be a symbol of order — n of 
type (1,0) on M" ; furthermore, on suppa, we assume: 

(i) for all multi-indices a there exists a constant > such that 

|9^M0l<c«(i + iei)'-'"'; 

(ii) there exist constants M,C > such that for all |^| > M we have |t(^)| > C|^|; 

(iii) there exists a constant Co > such that \d^T{\u)\ > Cq for all u G 
A > 0; m particular, | Vr(0| > Cq for all ^ G M" \ {0}; 

(iv) there exists a constant Ri > such that, for all A > 0, 

^S,(r) = G : r(0 = A} C 5^,(0) . 

Then, the following estimate holds for all R > 0, x G M", t > 1: 

e^(^<+^^^)^^a{09R{0 < Cr"^ , (4.13) 



where gji{^) is as given in ( [^. j^[ ) and C > is independent of R. 

Remark 4.9. For an integral of this type with some specific compactly supported 
function, x £ C^(M'^) say, in place of gn, we can just use the result for R = 0. In 
this way we obtain Corollaru \2.9[ . 

Proof. We may assume throughout, without loss of generality, that either r(^) > 
for all C, G M" or t{^) < for all ^ G M". Indeed, hypothesis (ii) and the continuity 
of r ensure that either t{C,) is positive for all |^| > M or negative for all |^| > M. In 
the case where t{^) is positive for all |.^| > M, set 

r+({) := r(0 + min(0, inf r(0) > for all ^ G M". 

Now, t(^)— r+(^) is a constant (in particular, it is independent of ^) and |e*['^(^)~'^+(^)l*| = 
1, so it suffices to show 



n-l 

< ct-— 
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In the case where r(^) is negative for |^| > M, set r(^) := — r(^) and by similar 
reasoning to above, it is sufficient to show 



< cr— 



where -t+(0 < for all ^ G M". 

We begin by dividing the integral into two parts: near to the wave-front set, i.e. 
points where Vg[x ■ ^ + T(^)t] = 0, and away from such points. To this end, we 
introduce a cut-off function k G C^(M"), < hi{y) < 1, which is identically 1 in 
the ball of radius r > (which will be fixed below) centred at the origin, Br{0), and 
identically outside the ball of radius 2r, i?2r(0). Then we estimate the following 
two integrals separately: 

JR" 

For l2{t,x) we have the following result: 

Lemma 4.10. Suppose a(^) is a symbol of order j G M. Then, for each I G N with 
I > n + j, we have, for all t > 0, 

\l2{t,x)\<Cr,lt-\ (4.14) 

where the constants Cr,i > are independent of R. 

Proof. In the support of (1 — K){t^^x + Vr(^)), we have |x + tVr(^)| > rt > 0, so we 
can write 

(x + tVr(0) _ ^ J .(x-i+rm] - pii^<+rm ■ 

z|x + tVr(OP ^ 
therefore, denoting the adjoint to P = -i^^^y^l^jp ■ by P*, we get 

hit,x)= I e'^-^^^^^^'\P*)'[a{09Rmi-^^){t-'x + Vr{0)]di 

JR" 

for each / G N. We claim that for each / there exists some constant C^,/ > indepen- 
dent of R so that, when t > 1, we have 

{P*y HOgnim - ^) {t-'x + Vr(0)] < a,.r'(l + 1^1)^-' ; (4.15) 
assuming this, we obtain, 

\l2it,x)\<Cr,lt~' [ -— 2— — rfe. 

JR" (1 + Kl) ^ 



Noting that dC, converges for / — j > n yields the desired estimate (14.141) . 
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It remains to prove fl4.15l) . Let / = f{C,;x,t) be a function that is zero for 
\x + tVr(^)| < rt and is continuously differentiable with respect to ^; then, 



i|a; + tVr(OrJ 2|x + tVr(0r i|a; + tVr(0|2 

2t(x + tVr(O) ■ [VMO ■ + ^Vr(O)] 



/• (4.16) 



z|x + tVr(0|^ 

Hence, using |x + tVr(OI > (hypothesis on /) and |9"r(0| < C(l + 
(hypothesis (i)), we have 

\P*f\<Crr'[{l + \^\)-'\f\ + \VJ\]. (4.17) 

Now, for all multi- indices a and for all ^ G M", we get 

• |5"a(0l < C'a(l + |el)^'"'°' for all ^ G as a G 5(o(^"); 

• \d^ [{I - K)(t-^x + Vr(0)] I < C„(l + |^|)-l°l, for all ^ G M'^— here we have used 
hypothesis (i) once more. Also, it is zero for each a when |x + tVr(^)| < rt by the 
definition of k. 

Furthermore, |9"^/j(0l = \d^9o{0\ < Ca(l + K|)"'"' for < < 1, since C^(M" \ 
{0}) C 5i%(M"). For i? > 1, we have: 

d'-gniO = d'-m/R)] = R-^''\d^g){i/R) and g G S%{W^) 

=^ |5"^7i?,(OI < c,i?^i"i(i + l^/i?l)^i"i < c,(i + l^l)-i°i . 

Therefore, 

|9"^r(0I < C'a(l + ICI)^'"' for all ieW and multi-indices a , (4.18) 

where the Cq, > are independent of R. 
Hence, by f l4.17p . we obtain 

\P*HOgR{m - ^){t^'x + Vr(0)]| < Crr\i + \^\y^' . 

To prove fl4.15p for / > 2 we do induction on /. Note that 

i(p*)'/i < + \^\y'\iPT'f\ + iva(^*)''V}i] • 

The first term satisfies the desired estimate by the inductive hypothesis. For the 
second term, repeated application of the properties of a(^), g{^) and (1 — K)(t~^x + 
Vr(^)) noted above to inductively estimate derivatives of (P*)' /, 1 — 2 

yields the desired estimate. This completes the proof of the lemma. □ 
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This lemma, with j = — n, means that it suffices to prove (14.131) for Ii(t,x), 
where \t~^x + 'Vt{^)\ < 2r. 

Let {'^£{^)}^^^ be a partition of unity in M" where ^ C°°(M") is supported 

in a narrow (the breadth will be fixed below) open cone Ke, i = 1, . . . ,L; let us 
assume that Ki contains the point e„ = (0, . . . , 0, 1) (if necessary, relabel the cones 
to ensure this) and also that each Ki, £ = 1,...,L, can be mapped onto Ki by 
rotation. Then, it suffices to estimate 

I[{t,x)= [ e^(^-«+-(«Wa(0^?i?(0*i(0«^(i-'a: + Vr(0)t/e, (4.19) 

since the properties of r(^), a(^), gn^^) and K^f^x + Vr(^)) used throughout are 
invariant under rotation. 



By hypothesis (iii) , the level sets Sa = {.^ G M" : t{^) = A} are all non-degenerate 
(or empty). Furthermore, the Implicit Function Theorem allows us to parameterise 
the intersection of the surface S';^ = jT^x and the cone Ki. 

K^nj:', = {iy,h,{y))--y^u} ; 

here U C IR"'^^ is a bounded open set for which p{U) = fl Ki where p{y) = 
{y, ^/T^^Jy\^), and hx : U ^ is a smooth function for each A > 0; in particular, 
each hx is concave due to r(^) satisfying the convexity condition, i.e. E'^j^ is convex 
for each A G M. Then, in the case that r(^) > for all ^ G M", the cone Ki is 
parameterised by 

K, = {iXy,Xhxiy)):\>0,yeU} , 
and when r(0 < for all ^ G M", 

K, = {{\y,\hx{y)):X<0,yeU} . 

Now, let n : Ki n T.'^ ^ be the Gauss map. 

By the definition of K,{t~^x + Vr(,^)), we have 

\r'x-i-Vr{^x))\<2r 

for each C,x & KiH T,'x that is also in the support of the integrand of (I4.19p . Hence, 
provided r > is taken sufficiently small, the convexity of T,'^ ensures that the 
points t~^x/\t^^x\ and — n(,^A) are close enough so that there exists z{X) G U (for 
each ^x ^ Ki nT,'^) satisfying 



n{z{X),hxiz{X))) = -t-^x/\r^x\ = -x/\x\ G §' 



n-l 
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Also, {—Vyh\{y), 1) is normal to S';^ at (y, hx{y)), so, writing x = {x', we have 
X _ i-\/yhxiz{X)),l) _^ Xn _ 1 



|x| |(-V,/.a(^(A)),1)| \x\ |(-V,/.a(^(A)),1)| 

and — 



x' _ -Vyhx{z{X)) _ XnVyhx{z{X)) 



\x\ |(-V,/ia(^(A)),1)| |x| 

therefore, —x' = Xn^ yh\{z{X)). We claim that x„, is away from provided the 
breadth of the cone Ki is chosen to be sufficiently narrow, so 

— = -V,/ia(^(A)). (4.20) 



To prove this claim, first recall that Yi'^ C Br-^{Q) for all A > (hypothesis (iv) ) 



and note that d^^T{^) is absolutely continuous on Bji-^{0) (it is continuous in 



taking Co > as in hypothesis (iii) , we get that 



there exists 5 > so that \r]^ — ?7^| < 6, where r]^,!]'^ G B^i^^O), 
implies \d^„r{r]^) - d^„T{r]^)\ < Co/4:. 

Then, fix the breadth of Ki so that the maximal shortest distance from a point ^ G 
Ki n (IJa>o ^a) {f^^n : /i > 0} is less than this 6, i.e. 



sup 



iinf|e-/ie„|:eeJ^in(|JS'A)l<5. 

l''^ A>0 J 



Now, observe that for any eW^, > 0, we have 

I tl > \di,Af^en)\ - \d^Me) - %r(/ie„)| - \^ + d^Me)\ • 

Choose ^° e Ki n T^'x fl supp[K(t"^x + Vr(^))] and ;U > so that - /ie„| < 6 and, 
hence, 

l%^(a-%^(/^e„)|<Co/4; 



also, by hypothesis (ni) , |9^,^r(/ie„)| > Cq, so 

\t~^Xn\ > 3Co/4-2r. 

Taking r sufficiently small, less than Co/8 say, (ensuring r > satisfies the earlier 
condition also) we get 

\xn\ >ct>0 (4.22) 

proving the claim. 

Before estimating 04.191] . we introduce some useful notation: by the definition of 
9r{0, <n>, when R > 1 

^ e supp^fij =^ Rco < 1^1 < Rci, 
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also, if < -R < 1, then there exist constants Cq,Ci > so that Cq < |^| < Ci for 
C, G supp(7/j. Thus, by hypotheses (i) and (ii) , there exist constants c'q,c[ > such 
that 

Rc'q < |t(^)| < Rc[ if i? > 1 and ^ G supp (yf/j, 
c'o < < c'l if < < 1 and ^ G supp Qr. 

Let G G be identically one on the set {s G M : Cq < s < c'^} and identically 

zero in a neighbourhood of the origin; writing TZ = max{R, 1), this then satisfies 



Also, for simplicity, write 



(4.23) 



this is a type (1,0) symbol of order — n supported in the cone Ki, and the 
constants in the symbolic estimates are all independent of R as each gniC,), -R > 0, is 
a symbol of order with constants independent of R (see f l4.18l) ). 

We now turn to estimating fl4.19l) . Using the change of variables ^ t-^ {Xy, Xhx{y)) 
and equality fl4.20p . it becomes 



I[{t,x) 



JU 



J[\x'-y+\x„hx{y)+T{\y,Xhxiy))t]Q^(^\y^ A/lA(y)) 



gniXy, Xhxiy))^iiXy, Xhxiy))K{t-'x + Vr(Ay, Xh^iy))) 



JU 



G{X/n)K{r^x + VT{Xy,Xhx{y))) 



diX,y) 



d{X,y) 



dydX, 



dy dX 



(4.24) 



where we have used T{Xy, Xhx{y)) = X (definition of Sa) in the last line. Here, note 
that 



d^ 



d{X,y) 



XI y 

XVyhxiy) dxlXhxiy)] 



X^~'{dx[Xhx{y)]-yVyhx{y)) 



where I is the identity matrix. Differentiating T{Xy, Xhx{y)) = X with respect to A 
in the first case and with respect to y in the second, gives 

y ■ V^.r(Ay, Xhxiy)) + dx[Xhx{y)]d^M^y, Xhx{y)) = 1 , 
AVer(Ay, Xhxiy)) + XVyhx{y)d^„T{Xy, Xhx{y)) = . 

Substituting the second of these equalities into the first yields 

{dx[Xhx{y)] - y ■ Wyhx{y))d^„T{Xy, Xhxiy)) = 1 . 

We claim that 



\d^„T{Xy,Xhx{y))\>C>0. 
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(4.25) 



To see this, first note that 

\di,A^y,^hx{y))\ > \d^„T{\fien)\ - \d^„r{Xfien) - d^„r{Xy, Xhx{y))\ 

where /i > is chosen as above so that |/ie„ — (y, hx{y))\ < 6; now, \d^^T{\fien)\ > Cq 
by hypothesis (iii) Also, by the Mean Value Theorem, there exists C, lying on the 
segment between {Xy, Xhx{y)) and A/ie„ such that 

|%r(A/ie„) - d^M^y,Xhx{y))\ < C\V^d^MO\^S < C\^\-'XS < C5; 



choosing 6 > small enough (also ensuring it satisfies condition (14.211) above) com- 
pletes the proof of the claim. Hence, 









d{X,y) 




dinr{Xy,Xhx{y)) 



< cx 



n-1 



(4.26) 



Also, note that this Jacobian is bounded below away from zero because |<9^„t(^)| < C 
for all ^ e M" (hypothesis (i)), which means that the transformation above is valid 
in Ki. 

Next, using the change of variables A = Ax„ = Xxnt in (14.241) . writing h{X,y) = 
hxiy) and setting x := t~^x (so Xn = t~^Xn), we obtain 



Ju 



G 



i)<5'+Vr(X,,J.fc(J.,j,)))_il^, 



t ^x^^dydX. 



Therefore, using |^^| < CA"~^|x„|~("~^)t~(""i) (by (g^S])) and recalhng that 
1^(^7)1 < 1; we have. 



n — 1 , , 71 — 1 



\i[{t,x)\ < cr—\xn 



dX , (4.27) 



where. 



I{X,J--z( J- 



u 



' Xnty^ Xnt^ixnfy 



With Theorem 14.11 in mind, let us rewrite this in the form of (14. ip : 

/(A, fi;z)= [ e'^^^y'^-'^^aoifiy, fih,{y))b{y) dy , 

with arbitrary A > 0, /i > and z G R"^^, where 
• ^{y, = h^{y) - h^{z) -{y- z)- Vyh^,{z)] 
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A 

t\Xn\ 



dy. 



with support contained in U. 



We shall show that the following conditions (numbered as in Theorem 14 . 1 1 and Corol- 
lary ST]) are satisfied by /(A,/i; z): 

(11) there exists a bounded set U C M""^ such that b G C^{U); 

(12) Im $(y, fi]z)>0 for ally eU, fi> 0; 

(13') F{p, u, fi; z) = ^{pu + z, fi] z), u G p > 0, is a function of convex type 7 

(see Definition I4.3p : 

(14) there exist constants Ca such that \dy[ao{fiy, fihf^{y))]\ < for all y G f/, p > 
and |a| < [^] + 1. 



Assuming for now that these hold, Theorem 14.11 (or, more precisely. Corollary 14.71 
states that, for all A > 0, /i > 0, 



|J(A,/i;2)| < C(l + A)-V < CX 
This, together with fl4.27p . gives 



n-1 
7 



, , ^ , n—1, , , n— 1 

\l[{t,x)\ < Ct-—\Xn\- — 



~ n-1 



A ^ d\\ 



then, setting v 


_ A 


we have 


\m,x) 




n-1 


1 < cr 






= cr 


n-1 

-< \Xn\ 



u'^G{u) du < cr- 



for alH > 1 



Here we have used that G is identically zero in a neighbourhood of the origin and that 
it is compactly supported and also (14.221) (|x„| >C > 0); also, note the constant here 
is independent of R. Since this inequality holds for J((t,x), it also holds for Ii{t,x); 
thus, together with Lemma [4.101 this proves the desired estimate (14.131) . provided we 



show that the four properties (II) (14) above hold. 



Now, clearly (II) holds automatically and (12) is true since hn{y) is real- valued, 
so Im $(y, fi; z) = for all y E U, fi > 0. 



For (13') and (14) , we need an auxiliary result about the boundedness of the 
derivatives of h\{y): 

Lemma 4.11. All derivatives of hx{y) with respect to y are bounded uniformly in y. 
That is, for each multi-index a there exists a constant G^ > such that 

\d;;hx{y)\<Ga forallyeU,X>0. 
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Proof. By definition, T{\y, \hx{y)) = A. So, 

{V^>r){Xy,Xh,{y)) + {d^^T){\y,Xh,{y))Vyh,{y) = \-'Vy[r{\y, Xh,{y))] = 
or, equivalently, 

V./^.(y) = -^M^^^- (4.28 



{d^„T)iXy,Xh,iy)) ■ 

Hypothesis [(i)]( I Si'r (01 < for all ^ G M") and (g^SD (|%r(A?/, A/ia(?/))| > 

C > 0) then ensure that {Vyhxiy)] < C for ally eU, X> 0. 

For higher derivatives, note that \{y,h\{y))\ < Ri by hypothesis 



IV 



so, using 

hypothesis (i) once more, for all multi-indices a, there exists a constant Ca > such 
that 

\{d^T)iXy,Xh,{y))\<C^X'-\-\. 

Then, differentiating fl4.28p . this ensures, by an inductive argument, that the desired 
result for higher derivatives of hx{y) holds, proving the Lemma. □ 



Returning to the proof of (14), note that, 

|a^"ao(OI <Ca(l + lel)~'"' for alUeM" 



since, a{^) is a symbol of order 



n-l 



n (see fl4.23p for its definition). Together with 



Lemma l4.11[ this ensures that dy[ao{fiy, fJ^hnly)) is uniformly bounded for all y & U, 
/i > and |a| < [^^] + 1 as required. 

Finally, we show (13') observe that for \p\ < 6', some suitably small S' > 0, 

F{p, u, /i; z) = hf,{puj + z) - h^{z) - puj ■ Vyh^{z) 

7+1 



k=2 \a\=k 

So, F{p,u!, p; z) is a function of convex type 7 if (using the numbering of Defini- 
tion USD 

(CT2) Etl\EH=k MWizP'^l > C > for all G S"-^, p>0,ze 
(CT3) \dpF{p, uj, p; z)\ is increasing in p for < p < 5, for each uj G p > 0; 



(CT4) for each k E N, dpF{p, u, p; z) is bounded uniformly in < p < 5', G § 
p > 0. 
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Condition (CT4) , follows straight from Lemma 14.111 The concavity of hfj,{y) 
means that 

dpF{p,u,p; z) = dp[h^{piu + z)] = cj^Hess hf_,{pLU + z)lu < 

for all < p < 5' and for each cu G p > 0, z G M""^; coupled with the fact that 

dpF{0,uj, p; z) = 0, this ensures Condition (CT3) holds. 
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Lastly, recall that, by definition, 7 > 7(Sa) for all A > 0, which is the maximal 
order of contact between and its tangent plane; furthermore, 7 is assumed to be 
finite; thus, for some /c<7 + l<oo, we have 

Now, dp[hf,{z + = J2\a\=k ^dyh^{z)uj°', so for some A; < 7 + 1, we have 



\a\=k 



> C>0 



for all u G §""2. Thus, condition (CT2) holds 



This completes the proof of conditions (II) (14), and, hence. Theorem 



□ 



5 Oscillatory integrals without convexity 



Theorem 14.81 requires the phase function to satisfy the convexity condition of Defini- 
tion I2.5j however, we will also investigate solutions to hyperbolic equations for which 
the characteristic roots do not necessarily satisfy such a condition. In this section we 
state and prove a theorem for this case. First, we give the key results that replaces 
Theorem 14.11 in the proof, the well-known van der Corput Lemma. We recall the 
standard van der Corput Lemma as given in, for example, |Sog93 Lemma 1.1.2], or 
in [5te93l Proposition 2, Ch VIII]: 

Lemma 5.1. Let $ G C°°(]R) be real-valued, a G C^(]R) and m > 2 be an integer 
such that $(J)(0) = /or < j < m - 1 and $('")(0) 0; then 



a(x) dx 



<C(1 + A) 



-l/m. 



for all A > 0, 



provided the support of a is sufficiently small. The constant on the right-hand side is 
independent of X and $. 

If m = 1, then the same result holds provided ^'{x) is monotonic on the support 
of a. 



5.1 Real- valued phase function 

In the case when the convexity condition holds the estimate of Theorem 14. 81 is given in 
terms of the constant 7; as in the case of the homogeneous operators (see Introduction, 
Section 11.21) we introduce an analog to this in the case where the convexity condition 
does not hold. Let S be a hypersurface in M"; we set 

7o(S) := sup inf 7(S; a, P) < 7(S) 

crGE P 

where 7(S; o", P) is as in Definition 12.61 

An important result for calculating this value is the following: 
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Lemma 5.2 ( |Sug96| ). Suppose S = {{y,h{y)) -.yeU}, he C°^{U), U C M""^ is 
an open set, and let 

F{p) = h{r] + puj) — h{ri) — pVh{r]) ■ to 
where rj e U , u e S"^^. Taking a = (77, h{ri)) G S, G S"^^ and 

P = {a + s{uj, Vh{ri) ■ u) + t(-V/i(r/), 1) G M" : s, t G M} , 

then 

7(S; (7, P) = min {A; G N : F^^\0) ^ O} =: -f{h; r], u) . 

Therefore, 

7(S) = supsup7(/i;?7,t^), 
7o(S) = supinf 7(/i;?7,u;) . 

Now we are in a position to state and prove the result for oscillatory integrals with 
a real-valued phase function that does not satisfy the earlier introduced convexity 
condition. This is a parameter dependent version of Corollary 12.131 

Theorem 5.3. Let a{E) he a symbol of order -^^ — n of type (1,0) on M". Let t : 
R" ^ W be smooth on suppa, set 70 := sup;^>o 7o(Sa(t)) and assume it is finite; 
furthermore, on suppa, we also assume the following conditions: 

(i) for all multi-indices a there exists a constant Ca > such that 

|9^M0l<Ca(i + iei)'-i"i; 

(ii) there exist constants M,C > such that for all \C,\ > M we have |r(^)| > C|^|; 

(iii) there exists a constant Cq > such that \di^T{Xuj)\ > Cq for all uj G S"~^, A > 0; 

(iv) there exists a constant Ri > such that, for all \ > 0, 

Then, the following estimate holds for all R> 0, x G M", t > 1: 



< Ct~~o 



where gji{^) is as given in ^4-12 ) and C > is independent of R. 
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Proof. We follow the proof of Theorem 14.81 as far as possible, and shall show how 
the absence of the convexity condition affects the estimate. Thus, as in the proof of 
Theorem 14.81 we may first assume, without loss of generality, that either t{^) > 
for all ^ G M*^ or t{C,) < for all ^ G M". We will always work on the support of a, 
so by writing ^ G M*^ we will mean G supp a. 
Divide the integral into two parts: 

h{t,x):= [ e'^^<^^'^^^'^a{09R{04t-'x + VriO)d^, 

JR" 

where k G C^(M"), < K{y) < 1, which is identically 1 in the ball of radius r > 
centred at the origin, Br{0), and identically outside the ball of radius 2r, i?2r(0). 
By Lemma 14.101 (which does not require the phase function to satisfy the convexity 
condition), we have 

\l2{t,x)\ < Crt-^^^'" for all t > 1. 

To estimate \Ii(t,x)\ we introduce, as before, a partition of unity {'^e{0}i=i 
restrict attention to 

where is supported in a sufficiently narrow cone, Ki, that contains e„ = 

(0, . . . , 0, 1). Parameterise this cone in the same way as above: with U C M""^, 



{{Xy, \h,iy)) ■.\>0,yeU} if r(0 > for all ^ G 
{{\y, \h,{y)) ■.\<0,yeU} if r(0 < for all ^ G 



Here the Implicit Function Theorem ensures the existence of a smooth function hx : 
U ^ M. for each A > 0, but there is one major difference: the functions hx are not 
necessarily concave, in contrast to the earlier proof. Using the change of variables 
^ I— {Xy, Xhx{y)) — note that 



< C < 



d{X,y) 



< CX 



n-l 



by the same argument as in the proof of Theorem 14. 8[ providing the width of Ki is 
taken to be sufficiently small — gives 

Jo Ju 

gniXy, Xhx{y))'^iiXy, Xhx{y))K{t-^x + Vr(Ay, Xhx{y))) dydX . 

Once again, let G G C^iM.) so that gni^i) = ^^(0^(^(0/^) (where 7^ = max(i?, 1)) 
and a{C,) = a({)5'_R({)^&i(0) which is a symbol of order ~ n supported in Ki and 
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with all the constants in the symbolic estimates independent of R. So, recalling that 
r(Ay, Xhxiy)) = A and writing h{\,y) = hxiy), we get 



Ju 



Ju 



G{X/n)K{r'x + Vr{Xy,Xhx{y))) 



M&yA^'y)+^-\( J-y, Xh.^^, 



diX,y) 



dy dX 



G 



TZXnt 



K X + Vr rA-y, jA-h{ rA-,y 



d{X,y) 

where x = tx and A = Ax„ = Xxnt. Thus, using \K,{f])\ < 1, we have 



x„ t dydX, 



i(x 



TZXnt 



dX 



(5.1) 



where 

/I 



At this point, we diverge from the proof of the earlier theorem since we cannot apply 
Theorem I4.lt instead, note that, for some b G C^(M"~^) with support contained in 
U, we have 



< 



a( J^y, J^hiJ^, y 

1 Xnt^^ Xnt \ Xnt^ ^ 



\x„\t 



Ky) dyi 



dy'. 



We wish to apply the van der Corput Lemma, Lemma 15. 1[ to the inner integral. Set 
/i; 2) := z ■ y + h^{y), which is real-valued, and consider the integral 



aoiy,iJ')Ky) dyi 



where ao(y, /i) := /x" ^^^"'°^a{ixy, fih^{y)). Recall that 

J:, = {{y,h,{y)):yeU}, 



,0)) =: m. 



so by Lemma [5. 2|, 

min |A; G N : d^^^^y, /x; z)\^^^^ ^ o} = -f{h^; 0, (1, 0, 

Fixing the size of U so that \dy^^^{y, fi; z)\ > e > for all ?/ G ?7 ensures that 
the hypotheses of Lemma 15.11 are satisfied. Thus, since the support of b is compact 
in W-\ is contained in U, and Oq is smooth, we obtain 



JX^{y,fj.;z) 



ao(y,/i)&(?/) dyi 



< CA"^/" 
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Carry out a suitable change of coordinates so that m = inf^7(/i^; 0,uj) (this is pos- 
sible due to the rotational invariance of all properties used); then, since m < 70 by 
definition, we have 



for all A such that 



A 



TZx nt 



G supp G (this is to ensure A is away from the origin) 



Combining this with (15.11) then gives the required estimate: 



\i[{t,x)\ < 



A 



TZXnt 



pec, ^ 

Jo 



□ 



6 Decay of solutions to the Cauchy problem 

Recall that we begin with the Cauchy problem with solution u = u{t,x) 



m—1 



Dru + Y,PjiD.)Dr'u + J2 E c.,rD:Dln = 0, t>0, 

j=l 1=0 \a\+r=l ("•-'-) 

Dlu{0,x) = fi{x) eG^{W), / = 0,...,m-l, xeM", 

where PjiO^ the polynomial obtained from the operator Pj^D^) by replacing each 
derivative D^^ = \dxf. by ^fc, is a constant coefficient homogeneous polynomial of 
order j, and the Ca^r are constants. In this section we will prove different parts of 
Theorem 12.181 



6.1 Representation of the solution 

Applying the partial Fourier transform with respect to x yields an ordinary differential 
equation for u = u{t,^) := e~^^'^u{t, x) dx: 

m m—1 

DTu + Y,PAODT~'u + Y, E c«,.rA'^2=0, (6.2a) 

3=1 1=0 \a\+r=l 

Dlu{0,O = M), / = 0,...,m-l, (6.2b) 

where (t,0 e [0, 00) x W and P,(0 are symbols of PjiD^). Let Ej = Ej{t,^), 
j = 0, . . . , m — 1, be the solutions to fl6.2al) with initial data 

Then the solution u of (16.11) can be written in the form 

m—1 

u{t,x) = Y,i^''Ej3^fj)it,x), (6.3) 

j=0 
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where 3^ and 5"^^ represent the partial Fourier transform with respect to x and its 
inverse, respectively. 

Now, as fl6.2al) . fl6.2bl) is the Cauchy problem for a linear ordinary differential 
equation, we can write, denoting the characteristic roots of (16. ip by ri(^), . . . , rm(0! 

■m 

k=l 

where Aj{t,C,) are polynomials in t whose coefficients depend on ^. Moreover, for 
each k = 1, . . . ,m and j = . . . , m — 1, the A^(t,^) are independent of t at points 
of the (open) set e M" : rfc(^) 7^ T/(^)V/ 7^ k}; when this is the case, we write 
Aj{t,^) = Aj{^). In particular, there exists M > such that if |^| > M, the roots 
are pairwise distinct. For Aj{^), we have the following properties: 

Lemma 6.1. Suppose ^ G Sk := G : Tfc(^) 7^ ti{^)WI 7^ k}; then we have the 
following formula: 

m—j—l 

4(0 = — ^^"r"^"'-^^'" — , (6.4) 
n (^Ko-r.(o) 

where means sum over the range indicated excluding k. Furthermore, we have, 
for each j = 0, . . . ,m — 1 and k = 1, . . . ,m, 

(i) Aj{C,) is smooth in Sk] 

(n) A';{0 = 0(|er^') as \^\ ^ 00. 
Proof. The representation (16. 4p follows from Cramer's rule (and is done explicitly 
in |Kli67] ): Aj(^) = where V := {tI~^{0)u=i Vandermonde ma- 

trix and Vj' is the matrix obtained by taking V and replacing the k^^ column by 
( P ..^01 . . . 0)T. 

j 

Smoothness of Aj{^) then follows by Proposition [3]1] and the asymptotic behaviour 
is a consequence of lPart 11 of Proposition 13.51 since (16. 4p holds for all |.^| > M. □ 

6.2 Division of the integral 

We choose M > so that all roots rfc(^), k = 1, . . . ,n, are distinct for |,^| > M. Let 
X = xiO e C;f=(M''), < x(0 < 1, be a cut-off function that is identically 1 for 
|(^| < M and identically zero for > 2M. Then (16.30 can be rewritten as: 

m— 1 m—1 

u{t,x) = Y,3'~\E,x3'f))it,x) + ^3^-1(^,(1 - x)3'f,){t,x) . (6.5) 
i=o j=o 
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Large |^|: The second term of fl6.5p is the most straightforward to study: by the 
choice of M, we have 

m 
k=l 

therefore, since each summand is smooth in M"', we can write 

m—1 



^> ,=0 k=l 



j=0 

m—1 m 



(2-) 

Each of these integrals may be studied separately. Note that, unlike in the cases of 
the wave equation, Brenner jBre75] . and the general m**^ order homogeneous strictly 
hyperbolic equations, Sugimoto Sug94| , we may not assume that t = 1. The — L'^ 
estimates obtained under different conditions on the phase function for operators of 
this type are given in Section 16.31 below. 



Bounded |^|: We turn our attention to the terms of the first sum in fl6.5l) . the case 
of bounded frequencies, 

1 r 

■J~\E,x^f){t,x) = -— e-«(5^e'-(«)*4(t,0)x(e)/(0^e- (6.6) 

Unlike in the case above, here the characteristic roots ri(^), . . . , rm(0 cire not neces- 
sarily distinct at all points in the support of the integrand (which is contained in the 
ball of radius 2M about the origin); in particular, this means that the 

may 

genuinely depend on t and we have no simple formula valid for them in the whole 
region. 

For this reason, we begin by systematically separating neighbourhoods of points 
where roots meet — referred to henceforth as multiplicities — from the rest of the re- 
gion, and then considering the two cases separately. In Section 16.91 we find V — L'^ 
estimates in the region away from multiplicities under various conditions; in Section [7] 
we show how these differ in the neighbourhoods of singularities. First, we need to 
understand in what type of sets the roots Tk{^) can intersect: 

Lemma 6.2. The complement of the set of multiplicities of a linear strictly hyperbolic 
constant coefficient partial differential operator L{Dt, D^), 

5 := G M" : r,(0 7^ r,(0 for all j ^ k} , 

is dense in M". 
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Proof. First note 

5 = {eeM":Ai(O^0}, 

where is the discriminant of L(r, ^) (see the proof of Proposition 13.41 for definition 
and some properties). Now, by Sylvester's Formula (see [GKZ94] ). Aj, is a polynomial 
in the coefficients of L(r, which are themselves polynomials in ^. Hence, A^, is a 
polynomial in ^; as it is not identically zero (for large |^|, the characteristic roots are 
distinct, and hence it is non-zero at such points), it cannot be zero on an open set, 
and hence its complement is dense in M". □ 

Corollary 6.3. Let L{Dt, D^) he a linear strictly hyperbolic constant coefficient par- 
tial differential operator with characteristic roots ri(^), . . . , Tm{Cj- Suppose, for k ^ I, 
that Aiki C M" is the set of all ^ such that Tfc(^) = r/(^). For e > 0, define 

Mil ■■= {eeM":dist(e,Alfc/)<£}; 

denote the largest z/ G N such that meas(A^|;) < Ce'^ for all sufficiently small e > 
by codimA^fci. Then codimAiki > 1- 

Proof. Follows straight from Lemma 16.21 the fact that Ai ki has non-empty interior 
(it is an algebraic set) ensures that its ^-neighbourhood is bounded by Ce in at least 
one dimension for all small e > 0. □ 

We can note that if L{Dt, D^) is not differential, but pseudo-differential in D^, 
the rest of the analysis goes through in a similar way, but we may need to assume 
that codimA^fci > 1. 

With this in mind, we shall subdivide the integral (16. 6p : suppose L roots meet 
in a set M with codimA^ = £; without loss of generality, by relabelling, assume the 
coinciding roots are ri(^), . . . , ri(^). By continuity, there exists an e > such that 
they do not intersect other roots tl+i, . . . , in Ai^. Furthermore, we may assume 
that dAi"^ G C^: for each e > there exists a set with boundary such that 
C Se and meas{Se\M^) ^ as e ^ 0. Then: 

1. Let XM,e ^ C°°(M") be a smooth function identically 1 on Ai^ and identically zero 
outside Ai^^; now consider the subdivision of (16. 6p : 

/ e'^<Ej{t,Omd^= I e'^'-^Ejit^OXMAOmd^ 

+ [ e^<E,{t,i){l-XMM)mdi; 

JB2m{0) 

for the second integral, simply repeat the above procedure around any root mul- 
tiplicities in 52Af (0) \ Ai"". 
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2. For the first integral, the case where the integrand is supported on spht off 
the coinciding roots from the others: 

B2m(0) 

L 



■«( ^ e-'=(«)*4(t,e))xA4,.(0/(0rfe (6.7) 

3. For the first integral, we use techniques discussed in Section [7] below to estimate 
it. 

4. For the second there are two possibilities: firstly, two or more of the characteristic 
roots tl+i(0) ■ ■ ■ ! ''"m(0 coincide in i?2A/(0) — in this case, repeat the procedure 
above for this integral. Alternatively, these roots are all distinct in i?2A/(0)\A1^ — 
in this case, it suffices to study each integral separately as the AKt,^) are inde- 
pendent of t, and thus the expression (16.41) is valid and we can write 



» m 

m „ 

,1 JBonriO) 



fc=L+l 



estimates for integrals of the type on the right-hand side are found in Section [6T9] — 
note that in this case we may use that the region is bounded to ensure that all 
continuous functions are also bounded. 

Continue this procedure until all multiplicities are accounted for in this way. 

Finally, let us recall the following result that can be found in |BL76t Theorem 6.4.5]: 
Theorem 6.4. Suppose T is a linear map such that it maps 

where Sq 7^ Si, 1 < Po)Pi < 00; then T also maps: 



where < 6 < 1 and 



11- 1 1 - 

— = + — , — = + — , se = {l-0)so + 0si. 

Pe Po Pi qe qo qi 

That is, 1 1 T/ 1 1^99 < C||/||^^6i and C is independent of f E . 
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In particular, this means that if we have estimates 

||T/|U^ < Ct'^-'ll/ll^ivo , IIT/IU. < Ct'^^ 11/11^^.. , 

then 

r/iu, <c(i+tr''ii/ii^^., 

where i + i = 1, 1 < p < 2, iV^ = 7Vo(J - J) + fiVi and = rfo(J - \) + fc^i- As 
usual, this reduces our task to finding — L°° and — estimates in each case. 



6.3 Estimates for large frequencies 

Via the division of the integral above, it suffices to find — L'^ estimates for integrals 
of the form 

where aj(^) = as |^| — >■ oo is smooth (or is zero in a neighbourhood of 

0), and r(^) is a complex-valued, smooth function which is 0(|^|) as |.^| — > oo and 
Imr(^) > for all ^ G M". Note that t{^) does not have to be homogeneous. 

By further judicious use of cut-off functions, we can split the considerations into 
the following cases of Theorem I2.18t 

1. t{^) is separated from the real axis, i.e. there exists 6 > such that Imr(.^) > 6 
for all 1^1 > M (Theorem EI]); 

2. t{C,) lies on the real axis (this case is contained in Theorems I2.3H2.121 since r is 
real valued); 

Let us look at each of these in turn. We will not consider the case of t{^) tending 
asymptotically to the real axis as |.^| — ^ oo since it is not part of Theorem 12.181 and 
since we do not have at present any examples of such behaviour. 



6.4 Phase separated from the real axis: Theorem 12.11 

In this section, we consider the case where characteristic root r(^) is separated from 
the real axis for large |^|; let us define 5 > to be a constant such that Imr(^) > 6 
for all 1^1 > M. Again, x is a cut-off to the region (which may be unbounded) where 
these properties hold. 

We claim that, for all t > 0, we have 



did: 



did: 



,i(x-5+r(0i) 



< Ce 



~5t\ 



|^^iVi + |a|+r-j , 



,i{x-i+T{C)i) 



L2 



< Ce 



-St 



w. 



o:\+r-j , 



where A^i > n, r > 0, a multi-index. Indeed, these follow immediately from: 
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Proposition 6.5. Let t : U ^ C be a smooth function, U C M" open, and aj 
djiO ^ ^loi^)- Assume: 

(i) there exists 6 > such that Imr(^) > 6 for all ^ & U; 

(ii) \r{0\<C{l + \^\)forall^eU. 
Then, 



,i(x-C+r(5)t) 



u 



< Ce 



-St 



Ni + \a\+r-j 



and 



u 



< Ce 



^5t 



for all t > 0, Ni > n, multi-indices a, r G M and f G C^{U). 
Note that in the case of r = 0, condition (ii) may be omitted. 

Proof. By the hypotheses on r(^) and aj(^), we can estimate 



U JU 



-ImT(5)i| 



5t 



U 



Aa\+r-j\ 



'U 



Ni + \a\+r-j 



1/(011 



< Ce~ ll^iVi + lcl+r-i 



-St 



This proves the first inequahty. For the second, note that Plancherel's theorem imphes 



L2(RS) 



a,(0rr(07(0| 



then. 



e^-«)*a,(orr(07(orrfe 



u 



< / e-2^--«)*ia,(orieri°i|r(ori/(or^e 



< Ce 



2(\a\+r-j) 



\fm'd^<Ce 



-2St\\ ru2 



U 



This completes the proof of the proposition. 



□ 
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We note that there may be different version of the L°°-estimate for the integral in 
Proposition I6.5[ For example, applying Cauchy-Schwartz inequality to the estimate 



< Ce 



-St 



(0 



\a\+r-j 



u 



u 



1/2 



2N[+2\a\+2r-2j 



u 



1/2 



established in the proof, we get 

r 

from which we obtain the estimate 

3*(-«+-(«)*)a,(0rr(07(0 < Ce^^i/ll (6.8) 

W2 

witlfl A'^^ > ^. Interpolating with the L^-estimate from Proposition 16.51 yields esti- 



mate (12. 7p in Section [2TT1 

From Proposition 16. 5[ by the interpolation Theorem 16.41 we get 



did: 



3^(-€+-(«)*)a,(0x(0/(0rf^ 



< Ce 



-5t 



LI 



Np + \a\+r-j , 



where i + ^ = 1, 1 < p < 2, Ap > n(i - i) , r > 0, a a multi-index and / G C^(M"). 
Thus, in this case we have exponential decay of the solution. This proves the first part 
of Theorem 12. 1[ The second part of the statement of Theorem 12. II is a straightforward 
consequence. 



6.5 Non-degenerate phase: Theorems 12.31 and 12.4 

In this section, we will prove Theorems 12.31 and 12.41 and discuss the behavior of critical 
points of the phase. In fact, we will prove Theorem 12.31 since the proof of Theorem 
12.41 can be given in the same way after restricting to a subset of variables on which 
the non-degenerate matrix A{C,^) is attained (possibly after a coordinate change). We 
will not write a further cut-off function x to a set U as in Theorems 12.31 and 12.41 to 
ensure that the results that we obtain are uniform over the positions of such sets U. 
However, we will keep in mind that we are only interested in the local in frequency 
region here, so all the integrals are convergent. So, we first consider the case where 
we have 

f e^(-«+-(€))*a(0/(Orfe, 

and det Hess t{^) ^ for all ^ G supp a. Here we denote x = t~^x. To estimate this, 
we first consider the oscillatory integral 

[ e*(^-«+^(«))*a(Of/e, 
^Here N{ does not have to be an integer. 
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where a = a(0 e S^^^, some G M, lmr(0 > for all ^ G M", and, for some ^° G M", 
X + V^r(^°) = and detHessr(^°) 7^ 0; we refer to as a (non-degenerate) critical 
point and we microlocalise around it. Let us assume that is the only such critical 
point — if there are more than one, we use suitable cut-off functions to localise around 
each separately (we assume the set of critical points has no accumulation points). 
Indeed, let ^ G C(f (M") be supported in a neighbourhood V of so that there are 
no other critical points in V. Then consider separately 



aiO^iOd^ and 



The second integral, which we may assume contains no critical points in its support 
(otherwise introduce further cut-off functions around those), can be shown to decay 
faster than any power of t: note that away from the critical points, we can use the 
equality 



X 



Vr(0 



so, integrating by parts repeatedly shows that for any G N sufficiently large. 



3'("-«+"(«))*a(0(l-?9)(0 



Let us return to the case when there is a critical point. We may assume that 
Imr(^°) = since otherwise Imr(^°) > in view of (12.21) . and then Theorem 12.11 
would actually give the exponential decay rate. We now claim that 



aiomd^ < cr"/VetHess(e°)r^/xe°)x(e° 



<Cr"/Vet Hess(e° 



^1/2/ 



\e\)-r 



(6.9) 



This is a consequence of the following theorem, see e.g. [Hor83a[ Theorem 7.7.12, 
p. 228]: 



Theorem 6.6. Suppose $ = $(x, y) G C°°(]R" 
function in a neighbourhood of the origin (0, 0) G 

• Im<l> > 0; 



X 



is a complex-valued smooth 
■^^ such that: 



• lm$(0,0) = 0; 

• *'.(0,0) = 0; 

. det $^,(0,0)^0. 

Also, suppose u G C^{K) where K is a small neighbourhood of (0,0). Then 



o:'-'^(^'y)u{x,y)dx- 



N-l 
j=0 



72 



for some choice of operators L^j, where the notation G^{y) {where G{x,y) is the 
function) means the function of y only which is in the same residue class modulo the 
ideal generated by d^/dxj, j = 1, . . . ,n. 

The proof of this result uses the method of stationary phase; similar results (with 
slightly differing conditions and conclusions) can be found in |Sog93 (1.1.20), p. 49], 
|Ste93l Ch. VIII, 2.3, Proposition 6, p. 344], |Dui96l Proposition 1.2.4, p. 14] and 
|Tre80l p. 432, Ch. VIII, (2.15)-(2.16)], for example. 

So, we have fl6.9p as a simple consequence of this theorem; now, in order to show 
that 



(6.10) 



we must choose /i G M suitably. In the sequel we may assume that M is even; if M 
is odd, the result follows by a standard interpolation argument taking the geometric 
mean. 

Assume that |det IIessr(,^)| > C{1 + |C|)~^^ for some M e M; then taking jj, = 
M/2, we have this estimate. This extends the case of Klein-Gordon equation (which 
is done in |Hor97j pp.146-155) where detHessr(0 = (1 + l^|)~""^ so M = n + 2. 

Let us now apply this result to our situation. We have 

where we may now think of as G S^q to ensure uniformity, and aj(0 = C^d'Cl""') 
as 1^1 — s> oo; we assume |det Hessr(^)| > C(l + 1^1)"^''^. Now, for each z/ G N, we have 



%(0 = (i + iei 



2\-u/ 



where aj^^) = Cq,(1 + |^n~''^°'aj(0 is of order —j — 2u+ \a\. Moreover, aj^a'& is 
of order —j — 2v+\a\ uniformly over -d (satisfying the necessary uniform symbolic 
estimates). Taking v = M/2— j and using that |a| < i/, we can ensure that the worst 
order of any of these symbols is —M/2. Then, 

f ^^i^<+rima^^^)^(^^f{^)d^= J2 [ e^^~^<^^^^^^\^{0mD^f{0di 



Then 



IE 

\a\<u 



\a\<i/ 

E 

\a\<iy 



a,AOmd^*Dy{x] 



a,AOmd^*Dy]{x 



\a\<iy 



ajAOmd^ 
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where we used estimate (16.101) . Thus, we have an — L°° estimate in this case. To 
find an — estimate is simpler: by the Plancherel's theorem, we have 



<c\\{o-'mh.<c\\f\\ 



w: 



Using the interpolation Theorem 16.41 and noting that all integrals are bounded for 
small t, we obtain Theorem 12.31 

Behaviour of Critical Points: Above, we assumed that S,^ was the only critical 
point of the phase function; this is not such an unreasonable assumption as the 
following observation shows: 

Lemma 6.7. // the matrix of second order derivatives Hessr(,^) is positive definite 
for all ^, then the integral 

has only one critical point. 

Proof. Suppose G M" are two such critical points. So x+V ^t{^£}) = x+'V^t{^'^), 
or d^.T{^^) = d^.T^^"^) for each j = 1, . . . ,n. Thus, by the fundamental theorem of 
calculus, for all j = 1, ... , n, we have 



Jo 



ds 



But this means that (^^ - ^'^) Hess t{^^ + s{^'^ - ^^))(^^ - ^'^) = for all s since the 
Hessian is positive definite; and since it is never zero, we have that — = 0, which 
shows that there is at most one critical point. □ 

An example of such an operator is the Klein-Gordon equation. 

Remark 6.8. In general, another consequence o/Hessr(^) being positive definite is 
that the level sets Sx = {C, E M."' : r(,^) = A}, A e M are all strictly convex; indeed, 
if we take a smooth curve C,{s) G ^a, s > 0, where C,{0) = and, by assumption, 
^(s) 7^ 0, then Vr(,^(s)) ■ C,{s) = (differentiate t{C,{s)) = X), and (differentiating 
again) 

i{sf ■ Hess r(e(s)) ■ i{s) + Vr(e(s)) ■ ^{s) = 0. 

Then, since Hessr(.^) is positive definite, the first term in this sum is positive, hence 
the second is negative — which means that the angle between Vr(,^(s)), that is, the 
normal to the level set, and ^{s) is strictly greater than 7i/2, so the level set is strictly 
convex. In particular, this shows that imposing the condition Hessr(.^) positive defi- 
nite is stronger than imposing the convexity condition of Definition \2.5\ and making 
it clear why we get a faster rate of decay in this case (see the next section for that 
case). 
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Remark 6.9. //rank Hess r(^) = n — 1, then a similar argument can be used to prove 
the corresponding part of Theorem \2.18\ i.e. that there is decay of order — This 
is a consequence of an extension to Theorem \6.6\ — see Hormander \Hdr83a\ Section 
7.7] . 



6.6 Phase satisfies the convexity condition: Theorem 12.81 

The case of real roots and real-valued phase functions subdivides into the following 
subcases, each of which yields a different decay rate: 

(i) detHessr(,^) 7^ 0; in this case we use the method of stationary phase in the 
same way as in Section [6751 with same result; 



(ii) det Hessr(^) = and t{^) satisfies the convexity condition of Definition I2.5j in 
this case we use Theorem I4.8t 

(iii) the general case when det Hessr(^) = (i. e. t{^) does not satisfy the convexity 
condition); in this case, we use Theorem I5.3[ 

We assume throughout that r(0 > for all ^ G M" or r(0 < for all ^ G M". This 
is valid because for the characteristic roots lying on the real axis, there exists a linear 
function r(^) such that T^iC,) '■= TkiO ~ ^(0 either everywhere non- negative or 
everywhere non-positive, and, if Tfc(^) satisfies the convexity condition, so does rfe(^). 
A proof for this in the case of homogeneous symbols is given in Sug94| and we recall 
this result here for completeness: 

Proposition 6.10. Let ^PkiO; k = 1, . . . ,m, be the characteristic roots of a strictly 
hyperbolic operator with homogeneous symbol of order m, ordered as ^Pi{^) > V^2(0 > 
■ ■ ■ > (fmiO fo'fi 7^ 0- Suppose that all the Hessians (p'l.{0 ^"^^ s emi- definite for 7^ 0. 
Then there exists a polynomial a{^) of order one such that fm/2{0 > '^(0 > fm/2+1 
{if m is even) or a{^) = ip(^rn+i)/2{0 {if is odd). Moreover, the hypersurfaces 
Sfc = G M"; ifk = ±1} with (pk{i) = ^k{0 — ^(0 (A; 7^ (m + l)/2) are convex and 
7(Sfc) < 2[m/2]. 

The generalisation of this proposition to the case of non-homogeneous symbols 
follows using the perturbation results in Section [31 

Assume that r(^) satisfies the convexity condition of Definition 12. 51 Set 7 = 
7(r) := sup;^>o 7(^a(t)), where, as before, 

S,(r) = G M" : r(0 = A} . 

and 

7(SA(r)) := sup sup 7(SA(r); a, P) 

aeS;,(r) P 

where the second supremum is over planes P containing the normal to ^\{t) at a 
and j{T,\{t); a, P) denotes the order of the contact between the line T^ fl P — T„ is 
the tangent plane at a — and the curve ^x{t) fl P. 

We have the following results which ensures that this is finite: 
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Lemma 6.11. Suppose t : M" — M a characteristic root of a linear m^^ order 
constant coefficient strictly hyperbolic partial differential operator. Then, there exists 
a homogeneous function of order 1, </?(0) characteristic root of the principal symbol, 
such that 

7(SA(r)) 7(Si(v9)) as A ^ oo . 
// we assume that ^{J1\{t)) < oo for all A > 0, then we have 7(r) < oo. 

Proof. This is true because: 

(a) by Proposition I3.5[ IPart III Sa(^) is near to ^^(v?) for large A in a suitable 
metric; 

(b) by the homogeneity of </:>, if |A — A'| is sufficiently small, then Tj\{ip) is near to 
Ea'(</?) for large A in the same metric; 

(c) Proposition [231 [ParOYl ensures that T^{t) is near to T^i^^) (because derivatives 
of r tend to those of ip) for large A; 

(d) so, with Sa(t) and T„{t) near to (in a suitable sense) the corresponding data 
of ifi for large A, it is clear that the 'j(T,x{t); a, P) is near to 7(Sa(</3); a, P), and 
hence 7(SA(r)) is near to 7(Sa(v5)); 

(e) finally, 7(Si(v?)) = 7(Sa(v5)) by homogeneity. 

□ 

In order to prove Theorem 12.81 we shall show that if G S^q is a symbol of order 
— j, then we have the estimate 



a.iOfiO < C(l + , (6.11) 



for all t > 0, where i + i = l,l<p<2, and/G C^(M"). The Sobolev order 
Np,j (which does not have to be an integer here) is worse for small times, being 
Npj > — ^) — j . It can be actually improved for large times, which will be done 
in estimate fl6.16p . 

Besov Space Reduction: We begin by following Brenner [Bre75j and also Sugi- 
moto |Sug94| in using the theory of Besov spaces and Paley decomposition to reduce 
this to showing, for all t > 0, the estimate 

:r-i(e-«)S(0*KO/(0) < C{l + t)-"-^i'^-^)\\f\\ ; (6.12) 

LI ''^p 

here {$i(0}z^o ^ Hardy-Littlewood partition: let $ G C(f (M") be such that 
supp<l> = |e G M" : i < 10 < 2| , $(0 > for ^ < 10 < 2, 

oo 

and J2 "^(S^^O = 1 for ^ ^ 0, 

fc=— oo 



76 



and set 

oo 

1=1 

Now, recall the definition of a Besov space, as given in, for example, Bergh and 
Lofstrom |BL76j : 

Definition 6.12. For suitable p, g, s G M define the Besov norm by 



1=1 



the Besov space B^^ is the space of functions in 5'(]R") for which this norm is finite. 

This result is the main one we shall need: 
Theorem 6.13 ( |BL76j . Theorem 6.4.4). The following inclusions hold: 

Bl^ C W; C Bl^ and B^, C W'^ C B^^ 

for all s eR, 1 < j> < 2, 2 < g < oo. 

There are some weaker versions of these embeddings for p = 1. Using this theorem, 
we have 



e^(-«+-«Wa,(0/(Orfe ^^^^^^ = (27r)-||9^-i(e-«)*a,(0/(e))(t,^)|| 



Li 



<C||3^-i(e-(«)*a,-(0/(0)Lo 

q,2 

CO 

1=0 

oo l+l 



|2 nV2 



2 ^1/2 



1=0 r=l-l 

in the final line we have used that Xl!-=l-i *^'r(0 = 1 supp by the structure of 
the partition of unity. Now, assuming that fl6.12p holds, this can be further estimated: 

i+i 



;=o r=l-l 



LI 



°° ■ 2n1/2 



/=0 r=l-l ^ 
oo l+l 

1=0 r=l-l 

oo 



1/2 
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Finally, using Theorem 16.131 once again, we get 



1=0 



2 N 2 

N 



1 



1=0 |a|<iVp,j 

^ 1/2 

H<Np,j 1=0 

c E P"/IIbo,<^II/II<..- 

\a\<Np,, 



2 \ 2 

LP 



Combining these estimates shows that (16.121) implies (16.111) as desired. So, it suffices 
to prove (16.121) : moreover, as shown above, this requires us to show two estimates 
and then interpolate — Theorem 16.41 yields: 



(6.13) 
(6.14) 



where Ni > n. 



— estimate: Since r(^) is real-valued and aj{^) = 0{\^\ ^) as |^| — * oo, by 
Plancherel's theorem we get 



w- 



Note that C is independent of / because is uniformly bounded in M". This 

proves the required estimate (16.141) . 

L} — L°° estimate: First, suppose < t < 1; then 



J\^\>M 
<^ll/ll<i-, 



(6.15) 



where Ni > n. 

For t > 1, we show 



n-l 

<Ct -I 11./ || ^_n^ 



(6.16) 
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Together (I6.15P and fl6.16l) will imply fl6.13p . We claim now that it suffices to prove 
that there exists a constant C > which is independent of / such that, for all t > 1, 



(6.17) 



Indeed, 



(2vr)"9^7ije-«)*a,(0<^>KO]*/(^) 



,j(x-5+r(0i) 



and, by the definition of the symbol of (D^), we have 
/ n \™- — - 



also, 

\\g * hWiac < \\g\\L'^ II^IIli , 

for all g e L°°(M"), /i e /.^(M"). Combining all these shows that fleTTD implies (16161) . 

In order to prove fl6.17p . we can use Theorem 14.81 as r : M" — > M is assumed to 
satisfy the convexity condition; let us check that each hypothesis holds. In addition 
to properties ensured by Proposition 13.81 we have: 

• Property (i) suffices for the hypothesis (i) of Theorem 14.81 to hold since aj (^) is 
supported away from the origin. 



-n+j 



is a symbol of order — n since a G S* and because it is zero 



in a neighbourhood of the origin. 
• the partition of unity {"^'KOIi^i is in the form of gR{C,) as required by Theorem 
Also, 7 < oo by Lemma [6 . 1 1 1 ab ove . Therefore, for t > 1, we get 



< cr 



Hence, we have fl6.16p . which, together with fl6.15p . proves fl6.13p : this completes the 
proof of Theorem 12.81 on real axis with convexity condition 7. 
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6.7 Results without convexity: Theorem 12.121 



The general case depends upon Theorem I5.3[ just as the case where the convexity 
condition holds depends upon Theorem 14.81 Here we assume that r is real valued. 
We introduce 70 = 70 (r) := sup;,>o 7o(SA(r)), where, 

7o(SA(r)) := sup inf 7(SA(r); a, P) 

(all notation as before). For this quantity we have the analogous result to LemmaEHH 
which can be proved in the same way: 

Lemma 6.14. If r : M" — > R a characteristic root of a linear m*^ order constant 
coefficient strictly hyperbolic partial differential operator, then, there exists a homo- 
geneous function of order 1, v?(0, a characteristic root of the principal symbol, such 
that 

7o(SA(r)) ^ 7o(Si((/))) as A ^ 00 . 
If we assume that 7o(SA(r)) < 00 for all A > 0, then we have 70 (t) < 00. 
We shall show 



^i(x-5 + T(5)t) 



1n \P II 



Li 



Np, > n(i - i) - J and 



for all t > 0, where ^ + ^ = I, I < p < 2, f e Q 
Nij > n — j. Similarly to fl6.16p . the Sobolev order Npj can be improved for large 
times. 

As in the case of Section 16. 6[ this can be reduced, via a Besov space reduction 
the interpolation Theorem 16. 4[ to showing 



< 0(1 + ty 

\3^-\e-^^^\{om)fmit,x)i,<c\\f\\ 



Ni-j 



\w:r-> ' 



where the partition of unity above and A^^i > n. 

The estimate follows by the Plancherel's theorem in the same way as before. 

For the — L°° estimate, the case < t < 1 is as in (16.151) : for t > 1 it suffices 
to show (see the earlier argument). 



This follows by Theorem 15.31 the hypotheses of this hold by the same arguments as 
above (see Proposition 13. 8p — the convexity condition is not required for the pertur- 
bation methods employed — and Lemma [6.141 This completes the proof of I2.12[ 
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6.8 Asymptotic properties of complex phase functions 

Here we consider what happens when the phase function t{^) is complex valued and 
look at its behaviour for large frequencies. In particular, this is related to the case 

lmr(0 ^ as 1^1 ^ oo. 

Unlike in the case of the phase function t{C,) lying on the real axis, here we do 
not consider a case where the phase function satisfies a "convexity condition". The 
reason for this is twofold: firstly, there is no straightforward analog of the convexity 
condition for real-valued phase functions as the presence of the non-zero imaginary 
part causes problems; secondly, there are no common examples of this situation, and 
hence it does not seem worthwhile developing a complicated theory for this situation. 

If detHessr(^) ^ 0, the analysis can be done in exactly the same way as that in 
Section 16. 5[ since Theorem 16.61 holds for integrals with complex phase functions. 

In general, we can derive certain properties of real and imaginary parts of r(^) us- 
ing perturbation arguments of Section [31 For example, for the index 70 = 7o(Rer) = 
supA>o7o(SA(Rer)) we can note the following: 

Lemma 6.15. If t : W"' ^ C is a characteristic root of a linear m*^ order constant 
coefficient strictly hyperbolic partial differential operator such that Imr(^) — * as 
1^1 00, then, there exists a homogeneous function of order 1, (f{$,), a characteristic 
root of the principal symbol, such that 

7o(SA(Rer)) 7o(Si(v9)) as \ ^ 00 . 

In particular, 7o(Rer) < 00. 

Proof. The hypothesis that the imaginary part goes to zero as |^| — > 00 implies that 
|r(^) — Rer(^)| — > as |^| — > 00. With this additional observation, the proof of 
Lemma 16.111 can then be used once more. □ 

In addition to Proposition l3.8[ we will now prove the following refined perturbation 
properties: 

Proposition 6.16. Suppose t : M."- —>■ C is a characteristic root of the strictly 
hyperbolic Cauchy problem U.l\) . Assume that it is a smooth function satisfying 
Imr(^) > 0. Assume also that the roots (pkiO' ^ — ^, ■ ■ ■ ,1^, of the principal part 
Lm are non-zero for all C, 7^ 0. Then we have the following properties: 

(i) for all multi-indices a there exist constants M, Ca, C'^ > such that 

|ae"Rer(0| <C,(l + |el)'-'"' 

and 

|a5Mmr(0l<c;(i + iei)-i'^i; 

for all 1^1 > M. 
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(ii) there exist constants M,C > such that for all |^| > M we have |Rer(^)| > 

ciei; 

(iii) there exists a constant Cq > such that |(9^ Rer(Au;)| > Cq for all ui G S"""*^ 
and sufficiently large A > 0; 

(iv) there exists a constant Ri > such that, for all sufficiently large A > 0, 

i{eeM":Rer(O = A}ci?«,(0). 

Proof. (i) The statements follow by Proposition 13.51 IPart Illl implies that for all |^| > 
and multi-indices a, 

\d^ReT{0\<\d^T{0\<m'~^^\, 

which suffices for the first part of (i). Furthermore, IPart IVl tells us that for all 
|,^| > and multi-indices a, 

\d^[ReT{0 - viO] + ^^^lmT{0\ = \d^r{0 - 9^X01 < ^1^1"'"' , 

where ^p{C,) is a characteristic root of the principal part (and is thus real- valued by 
definition of hyperbolicity) ; this implies that, for all |^| > iV and multi- indices a, 

|a|'[Rer(0-V^(0]l <C|er'"' and lmr(0| < . (6.18) 

The second of these gives us the second part of (i). 

(ii) We note that there exist constants C, C, C", M > such that, for all |.^| > M , 

|Rer(OI > |r(OI - |lmr(0| > C'\^\ - C" > C\i\ . 

Here we have used (13.191) . which did not require r to be real- valued (nor to satisfy 
the convexity condition), simply to be a characteristic root of a linear constant 
coefficient strictly hyperbolic partial differential equation, and the second part 

of dSISD. 

(iii) This follows in a similar way: using fl6.18|) . we have, for A > M, some M > 0, that 

|(9^Rer(Acj)| > \d^T{Xuj)\ - \d^lmT{Xu)\ > C" - C'A"^ > C . 

(iv) This follows from |Rer(0 - < C for all ^ e which holds in all M" by 
IPart III of Proposition 13.51 

□ 
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6.9 Estimates for bounded frequencies away from multiplic- 
ities 

In the following sections we find — L'^ estimates for integrals of the kind 

n 

where C is open and bounded, / G C^{W), a G C^{n), r E and 
Im r(0 > for all ^ E n. 

As in the case of large \^\, we can further split this into three main cases by using 
suitable cut-off functions: 

1. r(^) is separated from the real axis for all ^ ^ ^ (Theorem 12.11) : 

2. t{^) meets the real axis with order s < oo at a point E Q (Theorem 12.161) : 

3. r(^) lies on the real axis for all ^ G f2. 
We look at each in turn. 



6.10 Phase separated from the real axis: Theorem 12JJ again 

Similarly to the case for large \^\, we show that when the phase function t{^) is 
separated from the real axis (here, for ^ G fi, is a bounded set). 



did: 



< Ce 



-5t 



LI 



LP 



(6.19) 



where ^ + ^ = 1, l<p<2, r>0, aa multi-index, / G C^(]R"), 5 > is a constant 
such that Imr(e) > S for all ^ ^ ^ and C = Cfi^r,a,p > 0. So, in this case we also 
have exponential decay of the solution. 

By interpolating (Theorem 16. 4p . it suffices to show for such t{^) 



did: 



did: 



a(0/(0 



< Ce 



L2 



< Ce 



-St 



L2 



for t > 0, where r > and a is a multi-index. 

These are proved in a similar way to Proposition 16. 5[ but noting that the bound- 
edness of ^2 and the continuity in Vt of T{^Ya{^) ensure there exists a constant 
Cn^r,c, = C >Q such that |r(0r|a(0llel'''' < C for all ^eVL. Then, for alH > and 
r, a as above, we can estimate 



did: 



< C 



■i--(«)*ia(oiie||"||r(on/(oi^e 



< c 



[ e-^--(«)*|/(0|rfe<Ce' 



St 



L°°{n) 



< Ce 



-St 
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and 



n 

<Ce-''\\f\\r2<n.<Ce-'' 



lL2(n) 



We have now completed the proof of Theorem 12.11 

6.11 Roots meeting the real axis: Theorem 12.161 

In the case of bounded |^|, we must also consider the situation where the phase 
function r(^) meets the real axis. Suppose G f2 is such a point, i.e. Imr(^°) = 0, 
while in each punctured ball around B'^{S,^) C ^2, e > 0, we have Imr(^) > 0. 
Then, is a root of Imr(,^) of some finite order s: indeed, if were a zero of Imr(^) 
of infinite order, then, by the analyticity of Imr(^) at (which follows straight from 
the analyticity of t{^) at it would be identically zero in a neighbourhood of 
contradicting the assumption. 

Furthermore, we claim that s > 2, s is even, and that there exist constants 
Co, Ci > such that, for all ^ sufficiently close to we have 

coie-eT<|imr(oi<ciie-eT- 

Indeed, the Taylor expansion of Imr(^) around 



1=1 



is vahd for ^ e Be{i^) C Vt for some small e > 0. Now, if ^ e Be{^^), then + 2^° G 
Bir{^^) also. However, 



Imr(-e + = - 5^ % Imr(e°)te - + 0{\^- ^^'^ 

i=l 



thus, for e > chosen small enough, this means that either Imr(^) < or ImT(— ^ + 
2^°) < 0. In view of the hypothesis that Imr(^) > for all ^ G fi, we must have 
d^. Imr(^'^) = for each i = 1, . . . ,n. In conclusion, Imr(^) = 0(|^ — fo^' 
^ G which means that the zero is of order s > 2, and a similar argument 

shows that s must be even; also, this means that there exist Cq, Ci > so that the 
above inequality holds for ^ G B^iC,^), proving the claim. 

Now, we need the following result, which will be useful in the sequel. Moreover, we 
will give its further extension in Proposition 17.31 to deal with the setting of Theorem 
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Proposition 6.17. Let : f/ — * M, f/ C M" open, be a continuous function and 
suppose G t/ is such that 4>{^^) = and such that 0(^) > in a punctured open 
neighbourhood ofC,^, V\{^^}. Furthermore, assume that, for some s > 0, there exists 
a constant Cq > such that, for all ^ ^ V, 

0(0 > cole -eT- 

Then, for any function a{C,) that is bounded and compactly supported in U , and for 
allt>0, f e C^{W), and r eR, we have 

e-*(C)t|^ - eTl«(OII/(OI < C{1 + t)-^-^^y-^\\fU^ , (6.20) 

V 

and 

l|e-^«)*ie - er«(0/(OL.(v-) < ^(i + tr^ifh^ ■ (6.21) 

The constant C depends on U,V,Co and ||a||L°° , but not on the position of C,o- 

First, we establish a straightforward result that is useful in proving each of these 
estimates: 

Lemma 6.18. For each p,M > and a,c > there exists C = Cp^(j^M,c > such 
that, for allt>0, we have 

I xPe-^''''' dx < C(l + t)-(^+i)/- and sup x^e"'^^''* < C(l + t)""/'^ . 

Jo 0<x<M 

Proof. For < t < 1, each is clearly bounded: the first by and the second 

by M^. For t > 1, set y = xt^^"; with this substitution, the first becomes 

Mt^/" POO 

Jo 



while the second becomes 



sup yn-p'" e-^y" < t-p''' sup y^e-^^^ ■ 



These estimates imply those of Lemma [6 . 1 81 since both the integral and the supremum 
in the right hand sides are bounded. □ 



Proof of Proposition \6Al\ As for the proof of (l6.2Up . since a{^) is bounded in U by 
assumption, we have 



e-'^(«)*ie-eTia(Oii/(OMe<c / e-^«)*ie-e°ri/(OMe: 

V Jv 
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where V = V (1 supp a; this, in turn, can be estimated in the following manner using 
the hypothesis on 0(^) and Holder's inequality: 

v Jv 



'(V) 



Then, transforming to polar coordinates and using the Hausdorff- Young inequality, 
we find that, for some R> (chosen so that V C Bjj{$^^), which is possible since a(^) 
is compactly supported), we have 



Jv J§"-i Jo 



Finally, by the first part of Lemma 16.181 we find 

Jv Jo 

<C(l + t)-("+^)^||/|Ui. 

This completes the proof of the first part. 

Now let us look at the second part. By the second part of Lemma 16.181 we get 

Jv 

Now, it follows that 



f ^-co|«-,o|.,|j^^^|2^^ < sup|e--l^-«"l^*|||/||i.(^,) < C\\f\\ 
Jv v 



2 

L2 



Together these give the required estimate fl6.2ip . □ 
So, using this proposition, we have, for all t > 0, and sufficiently small e > 0, 



did: / e^(^-«+^«)*)a(0/(Oc?e 

JB.ie) 



L°°(RS) 



< / e-i--«)*|a(oiir(oneii/(OMe<c(i+t)-"/^ 



and, using the Plancherel's theorem, we get 
did: [ e^(^-«+^(«)*)a(0/(Ot^e 

JbJ£0) 



c||e-(^)V(o^ra(o/(OlL.(s^(50)) < cii/iu. 
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here we have used that |^|'"'|t(^)|'' < C on B^{^^) for r G N, a a multi-index. 
Thus, by Theorem 16 ■4[ for alH > 0, we get 



did: 



,i{x-C+r(Ot) 



LP(R5) 



< cil + ty 



-(---) 

s \ p q / 



Li 



(6.22) 



where l<p<2, ^ + ^ = 1. This completes the proof of Theorem 12.161 for roots 
meeting the axis with finite order and no multiplicities. 

Remark 6.19. If = 0, then Proposition 1 6'. i 7| further tells us that 

n ( 1 l\ \a\ 

< c{i + ty~yp^i)~~ 



did: 



L9(Mn) 



LP 



If, in addition, we have |r(^)| < ci|^ — ^''|'^\ for ^ near , then we also get 



did: 



n ( 1 l\ rsi 

< C(i + 



LP 



If both assumptions hold, we get the improvement from both cases, which is the esti- 

n { 1 l\ ]a\ rsi 

mate by C{1 + t) "^p i^ " 

From this, we obtain the statement of Theorem 12.161 in the frequency region 
-Be(^°)- Since there are only finitely many such points by hypothesis (H2) of Theorem 
12.161 hypothesis (HI) guarantees that on the complement of their neighborhoods we 
have ImTfc > 0. There we can apply Theorems 12.11 and 12.21 to get the exponential 
decay. In may happen that the roots are multiple, but Theorem 12.21 provides the 
required estimate in such cases as well. The Sobolev orders in Theorem 12.161 come 
from large frequencies as given in Theorem 12. 1[ This completes the proof of Theorem 
\21B and of Remark I^TTl 



6.12 Phase function lies on the real axis 

As in the case of large |^|, we can subdivide into several subcases: 

(i) detHessr(0 ^ 0; 

(ii) detHessr(^) = and r(^) satisfies the convexity condition; 

(iii) the general case when detHessr(,^) = 0. 

For the first case, the approach used in Section [631 can be used here also, since there 
we do not use that |^| is large other than to ensure that t{C,) was smooth; here, we are 
away from multiplicities, so that still holds. Therefore, the conclusion is the same, 
giving Theorem 12.31 

The other two cases are considered in the next section alongside the case where 
there are multiplicities since it is important precisely how the integral is spht up for 
such cases. 
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7 Estimates for bounded frequencies around mul- 
tiplicities 



Finally, let us turn to finding estimates for the first term of (\6.7\i . which we may write 
in the form 

L 
k=l 

where the characteristic roots ri(^), . . . , ri(^) coincide in a set 7W C of codimen- 
sion £ (in the sense of Section [2?T|) . C is a bounded open set and x ^ 

As before, we must consider the cases where the image of the phase function(s) 
either lie on the real axis, are separated from the real axis or meet the real axis. One 
additional thing to note in this case is that in principle the order of contact at points 
of multiplicity may be infinite as the roots are not necessarily analytic at such points; 
we have no examples of such a situation occurring, so it is not worth studying too 
deeply unless such an example can be found — for now, we can use the same technique 
as if the point (s) were points where the roots lie entirely on the real axis, and the 
results in these two situations are given together in Theorem I2.18[ We study this 
very briefly nevertheless to ensure the completeness of the obtained results. 

Unlike in the case away from multiplicities of characteristic roots, we have no 
explicit representation for the coefficients Aj{t,C,) (as we have in Lemma 16.11 away 
from the multiplicities), which in turn means we cannot split this into L separate 
integrals. To overcome this, we first show, in Section [7?T| that a useful representation 
for the above integral does exist that allows us to use techniques from earlier. Using 
this alternative representation, it is a simple matter to find estimates in the case 
where the image of the set Ai mapped by the characteristic roots is separated from 
the real axis (this is Theorem 12.21) and when it arises on the real clXlS cLS db result of 
all the roots meeting the axis with finite order, and these are done in Sections 17.21 
and [731 respectively. 

The situations where the roots meet on the real axis and at least one has a zero 
of infinite order there (either because it fully lies on the axis, or because it meets the 
axis with infinite order) is slightly more complicated; this is discussed in Section 17. 4[ 

7.1 Resolution of multiple roots 

In this section, we find estimates for 

L 

5^e^-'=«)*4(t,0, 

k=l 

corresponding to (12.81) . where Ti{C,) , . . . , tl{^) coincide in a set A4 of codimension i. 
For simplicity, first consider the simplest case of two roots intersecting at a single 
point, so that we have L = 2 and Ai = {C,^}', the general case works in a similar way, 
and we shall show how it differs below. So, assume 

Tiie) = r2{e) and r,(e°) ^ ni^) for A; = 3, . . . , m ; 



/ 

Jn 
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by continuity, there exists a ball of radius e > about C,^, B^l^'^), in which the only 
root which coincides with ri(^) is T2(^). Then: 

Lemma 7.1. For allt>0 and ^ G Bi;{^^), we have 

2 

|^e^"^(«)MJ(t,0 < C(l + t)e"™°(^'""^(«)'^'""^(«))*, (7.1) 

k=l 

where the minimum is taken over ^ G i?e(^*^). 
Proof. First, note that in the set 

5 := G : ri(0 7^ r^iO VA; = 2, . . . ,m andr2(0 7^ r^O V/ = 3, . . . ,m} 

the formula (16. 4p is valid for ^j(0 ^iid ^j(0- Now, recall that the sum Ej{t,C,) = 
e*'^''^^^*y4^(t, ,^) is the solution to the Cauchy problem fl6.2a|] . fl6.2cp . and thus is 
continuous; therefore, for all 77 G M*^ such that Ti{ri) 7^ Tfc(?7) and r2(?7) 7^ Tfc(?7) for 
= 3, . . . , m (but allow Ti{ri) = T2{ji)), we have 

2 

y e-.Wt^fc(i^^) = lin^ (e'-i(f)MJ(O + e'"2(«)*v4^(e)) , 
fc=i 

provided ^ varies in the set S (thus, ensuring e^'^^^^'^^ A}j{^) + e*'^2*^^^M|(,^) is well- 
defined). Hence, to obtain fl7.ip for all ^ G B^^i^^^), it suffices to show 

|gm(?)t^i(-^^ _^ e'''^^^^^*A^({)| < C(l + t)e~™"(^™'^i(^)'^™^2K))< 

for alH > 0, e e 5:(a = \ 

Now, note the following trivial equality: 

pi{yi-y2) _ p-i{yi-y2) pKyi-y2) 1 p-i(yi-y2) 

= K.e'y^ + — K.e'y^ 

2 2 

p-i{yi-y2) _ pi(yi-y2) p-i{yi-y2) 1 pi{yi~y2) 
+ ^ ^ K2e'y' + p K2e^y^ 

= sinh(|/i - y2)[Kie'y^ - Kse^^^] + cosh(?/i - ?/2)[i^ie*^2 + i^se^^^] . 
Using this, we have, for all ^ G -B^(^''), t > 0, 

= sinh[(ri(0 - r2(0)t](e^^^(«M)(0 - e^^^(«)*A2(0) 

+ cosh[(ri(0 - T2mKe'^'^^'"A]iO + e'^^^^^'A^O) ■ (7-2) 

We estimate each of these terms: 
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"sinh" term : The first term is simple to estimate: since 

sinh[(ri(0-r2(0)t] _ , ... ^ .... _^ . 

or, equivalently, as ^ — * through S, and /l^(^)(ri(^) — r2(^)) is continuous in 
Bei^°) for A; = 1, 2, it follows that, for all ^ G t > 0, we have 

|sinh[(ri(0 - r2(0)t](^](0e^^^(^)* - Aj(Oe^^^«)*) | 

< Ct[|e*"'2(^)*| + |e*^^(^)*|] < C'te"""''^^"'^^^^)'^™""^^^))* . (7.3) 

"cosh" term : Estimating the second term is slightly more complicated. First, 
recall the explicit representation (16.41) for the ^4^(0 &t points away from multi- 
plicities of Tk{C,) 



4(0 



m—j—1 
l<si<---<Sm-j-i<rn q=l 



n (^Ko-r.(o) 

So, we can write 

cosh[(n(0 - T,{OmA]iOe'^'^^^' + A2(0e-^(«)*) 

cosh[(ri(0 - T^m] e-^(€)*i^y,(0 - e-^(€)*F,yi(0 

nr=3(^^'(o - nmrkio - r^io) nio - t,{o 

where 

(m—i \ m 

n^^.(^) n (^^(o-r.(o)- 
l<si<---<Sm_i<m (?=! / k=l,k^p,(T 

Now, (cosh[(n(0 - r2(0)t])/(nr=3(^fc(0 - n(0)(^fc(0 - r2(0)) is contin- 
uous in S*, hence it is bounded there, and, thus, absolutely converges to a 
constant, C > say, as i ^ through S. This leaves the [e*^2(?)ii?i|2^(^) _ 

e-^«)*i^^i(0]/(ri(0-r2(0)term. 
For this, write 

m— 1 

where the Q^^'iiO polynomials in the Tfc(^) for k ^ p.a (which depend on i); 
also, note QllliC) = Q^iiO- Then, we have 



■ (7.4) 



n(e)-r2(0 

_ EZo K:?+i(0(^2(0'^e-«)*-n(0''e-^?)*)] 



^(0-^2(0 
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Let us show that this is continuous in B^{^^) and is bounded absolutely by 
(-;^g-mm{Ai,A2}. f^j. ^ g^j^j fQj, T , s ^ N , t > , have 

yi - 2/2 

y^^yl^e^y^^ - e^y^') ^ yle'y'\y\-yl) ^ (^^''yliyl-yl) 

yi - 2/2 2/1 - 2/2 2/1-2/2 

Furthermore, for all yi,y2 & C, t ^ [0, oo), s G N, 



2/1 - 2/2 



< Cnte-™'^^^"''^^'^'^^^)* and 



2/1 -2/2 



2/1 - 2/2 



for some constants Co,Cs. Using these with yi = ri(^), y2 = T2{^), r = k, 
and s chosen appropriately for Q^'J+ilO) the continuity and upper bound follow 
immediately. Thus, for all ^ G -B^(^°), t >0, 

|cosh[(ri(0 - r2(0)t](4(0e^^^«)* + A^^(Oe-^(«)*)| 

< (^^g-mm(Imri(0,ImT2(0)i (7 5) 

Combining (17. 2p . (17.31) and (17.51) we have (17.11) . which completes the proof of the 
lemma. □ 

Now we show that a similar result holds in the general case: suppose the charac- 
teristic roots ri(^), . . . ,tl{^), 2 < L < m, coincide in a set 7W, and that ti(^) 7^ Tfc(^) 
for all ^ G 7W when = L + 1, . . . , m. By continuity, we may take e > so that the 
set A4'^ = G M" : dist(^,7W) < e} contains no points 77 at which ri(?7), . . . , Ti(?7) = 
Tkirj) for k = L + 1, . . . ,m. With this notation, we have: 

Lemma 7.2. For all t > and ^ G A4^, we have the estimate 

L 
fe=l 

where the minimum is taken over ^ G A^"^. 

Note that this estimate does not depend on the codimension of M.. 

Proof. First note that, just as in the previous proof, for all 77 G M" such that 
ri(?7) . . . , Tl{vi) 7^ Tkijf) when k = L+1, . . . , m (but allowing any or all of ri(?7), . . . , Ti{ri) 
to be equal), 

J2 e*"^(^)*4(t, r/) = lim (e*"i(«)*A](0 + ■ ■ ■ + e*"^(«)*A^^(0) , 
fc=i 

provided ^ to varies the set S := IJ^i "S";! where 
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to ensure that each term of the sum on the right-hand side is well-defined. Note 
that Lemma 16.21 ensures every point in Ai is the limit of a sequence of points in S 
in the case of differential operators. Thus, we must simply show, for all t > 0, 
^ G {M'^y = M'^ \ M, that we have the estimate 

+ . . . + e'^^^^^'Af{^)\ < C(l + t)^-ie"*°^''^'==i--^^^"^(«) . 

Now, we claim that we can write ELi 0> for C G (M^)' and t > 0, 

as a sum of terms involving products of ^^~2^^ ^^^^ cosh terms of differences of 
coinciding roots; to clarify, (17.21) is this kind of representation for L = 2, while for 
L = 3, we want sums of terms such as 

sinh[«i(ri(0 - r2{0)t] cosh[«2(ri(0 - T^{0)t] sinh[a3(r2(0 " ^3(0)^] , 

where the are appropriately chosen constants; incidentally, a comparison to the 
L = 2 case suggests that the term above is multiplied by 

in the full representation. 

To show this, we do induction on L; Lemma [7TT] gives us the case L = 2 (note that 
the proof holds with C,^ and -Be(^'') replaced throughout by Ai and A4^, respectively). 
Assume there is such a representation for L = K < m — 1. Observe, 

K+l K K+l 



k=l k=l k=l,k^K 



. K+l 



K 

k=2 

by the induction hypothesis, there is a representation for each of these terms by 

2 

sinh[afc,/(rfc(0 - n{^))t] and cosh[p k,i{^k{0 - niO)^ terms. 



means of products of 



where 1 < k,l < K + 1 and the ak^u(3k,i are some non-zero constants. Next, note 
that we can write (ti(,^) — T2{C)) (o^'j indeed, the difference of any pair of roots from 
ri(^), . . . , tk+i{C}) as a linear combination of the ^^^^^"f differences Tk{0 — such 
that l<k<l<K + l] that is 

sinh[ai,2(ri(0 - T2{i))t\ = sinh [ ^ «fe,/(T-fc(0 - n{0)i 

l<k<l<K+l 

for some non-zero constants ;; similarly, there is such a representation for cosh[/3i^2(Ti(0' 
T2(0)^]- Lastly, repeated application of the double angle formulae 

sinh (a ± 6) = sinh a cosh b ± cosh a sinh b , 
cosh (a ±b) = cosh a cosh b ± sinh a sinh b , 
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yields products of '^^^"'"^^ terms, which completes the induction step. 

Now, as in the previous proof, each of these terms must be estimated. The key 
fact to observe is that 

L 

4(0 n (^Ko-r/c(e)) 

l=l,l^k 

is continuous in for all /c = 1, . . . , L. Then, using the same arguments as for each 
of the terms in the earlier proof, and observing that the exponent of t is determined 
by the products involving either 

(a) (sinh[afc,i(rfc(0 - Ti{^)t)]) / {Tk{Cj - n{0) terms, or 

(b) (e^"^(«)* - e*"'(«)*)/(rfe(0 - r^O) terms (see ([73])), 

the estimate (17.61) is immediately obtained. □ 

7.2 Phase separated from the real axis: Theorem 12.21 

We now turn back to finding — L'^ estimates for 

L 



-'^ k=i 



when Ti{C,) , ■ ■ ■ , tl{0 coincide in a set of codimension £; choose e > so that 
these roots do not intersect with any of the roots tl+i{^), . . . , rm(0 The set 

fl is bounded, and we may take x ^ C^(7W^). 

In this section (under assumptions of Theorem 12.21) . we assume that there exists 
S > such that Imr,fc(^) > 6 for all ^ G Ai"^ — so, min^ lmTk{^) > 6 > 0. For this, we 
use the same approach as in Section [6.10] but using Lemma UTZl to estimate the sum. 
Firstly, the — L°° estimate: 



•^^ k=i 

L 

( ^ e^-^(«)*4(t, Ork{ir)Cx{im) dx 



k=l 



L°°(RS) 



<maxsup|r,(Or / Ij^ e^^''(«)*4(t, \^\^''W)\dx 
< C{1 + t)^-^e-^i/|Uo.(^.) < C(l + t)^^'e^''\\fU. . 
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Similarly, the L? — estimate: 



k=l 
L 



L2(IRJ) 



k=l 



L2{n) 



< C{1 + t)^-^e-^i/|U.(n) < C{1 + tr~'e-'mfh2 . 



L~l-St\ 



Then, Theorem 16.41 yields 



•^^ k=i 



<C{l + tf-'e-''\\fh. 



where - + - = 1, l<j9<2. Once again, we have exponential decay. This, together 
with fl6.19p gives the statement when there are multiplicities away from the axis and 
completes the proof of Theorem 12. 2[ 



7.3 Phase meeting the real axis: Theorem 12.151 

We next look at the case where the characteristic roots ti{^) , . . . , tl{^) that coincide 
in the set A4 of codimension i meet the real axis in Ai with finite orders. If 
there are more points in A4 at which the above roots meet the axis with finite order 
(or even with infinite order/lying on the axis), they may be considered separately in 
the same way (or using the method below where necessary), while away from such 
points, the roots are separated from the axis, and the previous arguments and results 
of Section 12.11 may be used. 

Since the characteristic roots are not necessarily analytic (or even different iable) 
in Ai, we must look at each branch of the roots as they approach the real axis; 
set Sk to be the maximal order of the contact with the real axis for Tfc(,^), that is, the 
maximal value for which there exist constant cq > such that 

codist(e,^fc)''^ < lmrfc(0, 

for all ^ sufficiently near Zk, where = G M" : Imrfc(^) = 0}. By assumptions of 
Theorem 12.151 we have the estimate 

codist(e,A^)' < lmrfc(0, 

for some cq > and s > max(si, . . . , sl), for ^ close to A4. We will need the following 
extension of Proposition 16.171 Its proof is similar to the proof of Proposition 16.171 if 
we consider the coordinate system associated to A^. As usual Ai" = G : 
dist(^, A^) < e}. 
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Proposition 7.3. Let U C M" be open and let (j) : U W be a continuous function. 
Suppose A4 C U is a set of codimension i such that 

codist{^,Mr<m, 

for some cq > 0, and all C, G for sufficiently small e > 0. Then, for any 
function a{^) that is bounded and compactly supported in U, and for all t > 0, 
f e C^{W), and r eR, we have 



e-<^«)Mist(e,-M)1a(0||7(0Me < Cil + ty 



[£+r)/s\ 



and 



distil, MY aiOfiOl 



< c{i + ty/' 



L2 



The proof of this proposition is similar to the proof of Proposition 16.171 and is 
omitted. Theorem 12. 151 states that we must have the estimate fl2.12p . which is 



did: 



k=l 



< Cil + ty 



'VP 9 / 



LP 



By Lemma 17.21 and Proposition 17.31 to estimate the sum in the amphtude, for all 
t > 0, we have 



Or ^ 



< c 



k=l 



<c [ (i+t)^'^e-*--''-^---i--^«)|x(oii/(e)|rfe 



<C(l + t)^-i~(^/^)||/llz.i. 
Also, using the Plancherel's theorem, we have 



L2(RS) 



did: ( J^^ e^-« ( g e*-^(«)*4(t, o) xiom d^ 

L 



L2(R5) 



k=l 



< (J^l _|_ ^^-f^-l llg-* minfc=i....,L Imrfc(5) 



xiowmi" 



< C{l+t) 



L-1 



L2 
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Therefore, interpolation Theorem 16.41 yields, for all t > 0, 

L 

< 0(1 + ty 



^ k=l ' 



'VP q I 



LP 



where - + - = 1, l<p<2; this, together with (16.221) proves Theorem 12 . 1 51 for roots 
meeting the axis with finite order. 



7.4 Phase function on the real axis for bounded frequencies 

Recall that in the division of the integral in Section 16. 2^ we have 

m 

e-e(^e^-^«)*4(t,0)/(0^e, 
fc=i 

which we then subdivide around and away from multiplicities. The cases where the 
root or roots are either separated from the real axis or meet it with finite order have 
already been discussed; here we shall complete the analysis by proving estimates for 
the situation where a root or roots lie on the real axis. These results can be also 
applied to the case of multiple roots. 

We note that in the case of nonhomogeneous symbols this analysis is essential 
since time genuinely interacts with frequencies. Unlike in the case of homogeneous 
symbols in Section II. 2[ where one could eliminate time completely from estimates by 
rescaling, here it is present in phases and amplitude and causes them to blow up even 
for low frequencies. Thus, we must carry out a detailed investigation of the structure 
of solutions for low frequencies, and it will be done in this section. 

A number of estimates can be already obtained using our results on multiple 
roots from Section 17. 1[ To have any possibility of obtaining better estimates, we 
must impose additional conditions on the characteristic roots at low frequencies — for 
large |^|, these properties were obtained by using perturbation results, but naturally 
such results are no longer valid for |^| < M. Also, we can impose the convexity 
condition on the roots to obtain a better result than the general case. We will give 
different formulation of possible results in this section. 

Again, throughout we assume that either > for all C, or t{C,) < for all C,- 
The key point is to use a carefully chosen cut-off function to isolate the multiplicities 
and then use Theorem 14.81 or Theorem 15.31 to estimate the integrals where there are 
no multiplicities (and hence the coefficients Aj{t,C,) are independent of t) and use 
suitable adjustments around the singularities. For these purposes, let us first assume 
that the only multiplicity is at a point G -B2m(0) and ti(^°) = T2(^°) are the only 
coinciding roots, and let x be a cut-off function around Then, we must consider 
the sum of the first two roots, where we have a multiplicity at C,^, 

1= / e^-«(^e^-^«)*4(t,0)x(0/(0rfe, (7.7) 
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and terms involving the remaining roots, which are all distinct, 

k=3 >'-S2i\f(0) 

7.4.1 Case of no multiplicities: Theorem 12.81 

For the second of these integrals II, we wish to apply Theorem 14.81 if Tfc(^) satisfies 
the convexity condition, and Theorem 15.31 otherwise. 

In order to ensure the hypotheses of these theorems are satisfied, however, we need 
to impose an additional regularity condition on the behaviour of the characteristic 
roots for the relevant frequencies (i.e. ^ G -B2A/(0)) to avoid pathological situations: 

Assume \d^Tk{Xuj)\ > Co for all uj G §"^\ 2M > A > 0. (7.8) 

Since this is satisfied for large |.^| (see Proposition 13. 8p and always satisfied for roots 
of operators with homogeneous symbols, it is quite a natural extra assumption. 



The other hypotheses of these theorems hold: hypothesis (i) is satisfied because 



\d^Tk{0\ — for all ^ since the characteristic roots are smooth in M"; hypothesis (ii) 



only requires information about high frequencies; and hypotheses (iv) holds by the 
same argument as for large |^|, where only IPart III of Proposition 13.51 is needed, 
and that holds for all ^ G M". Also, the coefficients ^^(0 smooth away from 
multiplicities, so the symbolic behaviour (i.e. decay, or bounded for small frequencies) 
holds. 

Now — L°° and — estimates can be found as in the case for large |^|, and 
the interpolation theorem used to give the desired results. Thus, with condition (17.81) . 
we have proved the on axis, no multiplicities case of Theorem 12.81 

7.4.2 Multiplicities: shrinking neighborhoods 

Now we can turn to the other integral given by (17.71) . Here we will analyse what 
happens in certain shrinking neighborhoods of multiplicities. First we will assume 
that only two roots intersect at an isolated point, and then we will indicate what 
happens in the general situation. 

To continue the analysis of an isolated point of multiplicity as in (17. 7p . we intro- 
duce a cut-off function ip G C^([0, oo)), < ip{s) < 1, which is identically for s > 1 
and 1 for s < |; then (17.71) can be rewritten as the sum of two integrals I = Ii + I2, 
where 



h = {2rcr^ [ e^^-«^(tie-e°i)x(oE4(^'Oe*^^^^^*/(0^e, 

k=l 

/2 = (27r)-" / e^^<{l-m\i-e\)x{i)Y.A%i)e"-^^^'mdi. 

k=i 

We study — L°° estimates for Ji and — estimates for both Ji and I2 in this 
section. 
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— estimates: For this, we use tlie resolution of multiplicities technique of 
Section mi By Lemma [7. we have, in particular, 

2 

k=l 

for 1^ — < t~^. Now, we may estimate the integral using the compactness of the 
support of ip{s): for < t < 1, Ji is clearly bounded; for t > 1, we have 



h\<ct / 



This argument can be extended to the case when L roots meet on a set of codimension 
i. In the following proposition we will change the notation for the cut-off function to 
avoid any confusion with point multiplicities in the case above. 

Proposition 7.4. Suppose that L roots intersect in a set M. of codimension i. Let 
M' = {i^W : dist({, A^) < e}, and let 6 e C^{M') for sufficiently small e > 0. 
Then we have the estimate 



e''<e{t dist(e, M)) J2 Oe'^'^^^'m d^ 



k=l 



< C{l + t) 



(7.9) 



Proof. By using Lemma [7.21 in the (bounded) neighborhood Ai"^ of Ai, we obtain 



k=l 



The size of the support of 6{t dist(.^, Ai)) can be bounded by (1 + 1) ^, which implies 
estimate dH]). □ 



— estimates: Let us now analyse the L^-estimate. This analysis will apply 
not only in a shrinking, but in a fixed neighborhood of the set of multiplicities. We 
will discuss first the case of two roots intersecting at a point in more detail, thus 
analysing mainly integral / in (17. 7p . We can have several versions of L^-estimates 
dependent on conditions on multiplicities and on the Cauchy data that we can impose. 
For example, by Lemma [7. II and Plancherel's theorem we get 

||/|U2<C(l+t)||/|U2. (7.10) 

On the other hand we can improve the time behaviour of the L^-estimate (I7.10p 
if we make additional regularity assumptions for the data. For example, we can 
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eliminate time from estimate flT.lOp if we work in suitable Sobolev type spaces taking 
the singularity into account. Let us rewrite 



= (2vr) 
(2vr)-" 



fc=i 



x(0 



k=l 



Using the representation from Lemma 16.11 we see that the expression in the square 
brackets is bounded. Hence by the Plancherel's theorem we get that 

ii/iu^ < ii(n(o - Mor'xmmL^ = ik^p) - r,iD)r'xiD)fu.. (7.11) 

An example of this is the appearance of homogeneous Sobolev spaces for small fre- 
quencies in the analysis of the wave equations, or more general equations with homo- 
geneous symbols. For example, in the case of the wave equation we have ti(^) = |^| 
and T2{^) = — ICI, so that fl7.11l) means that we have the low frequency estimate for 
the solution of the form 

with the homogeneous Sobolev space H~^. 

In the case of several roots intersecting in a set A4, we have similarly: 

Proposition 7.5. Suppose that L roots intersect in a set Ai. Let Ai'^ = G : 
dist(^,7Vl) < e}, and let 9 G C^{Ai'^) for sufficiently small e > 0. Let J denote 
the part of solution corresponding to these roots microlocalised near the set M. of 
multiplicities: 



k=l 



J{t,x)= / e'^<e{0j2A^{t,0e'^''^^^'mdC 
Then we have the estimate 



\L-1 



(7.12) 



Moreover, let us assume without loss of generality that intersecting L roots are labeled 
by Ti, - ■ ■ ,tl. Then we also have 



l<l<k<L 



< C\\f\\L2(Rn). 



(7.13) 



L2(R5) 



Estimate fl7.12p follows from Lemma 17.21 and Plancherel's theorem. Estimate 
(17.131) follows from Plancherel's theorem and formula (16.41) . 

Interpolating between Propositions 17.91 and 17.121 we can obtain different versions 
of the dispersive estimate in a region shrinking around Ai, depending on whether we 
use frmil or fl7:T3D . 
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7.4.3 Multiplicities: fixed neighborhoods 

Here, for simplicity, we will concentrate on the case of two roots ri and T2 intersecting 
at an isolated point We will discuss both — L°° and — estimates under 
additional assumptions on the roots ri and T2. 

£i _ j^oo estimates: For I2 we are away from the singularity, so we can use that 
2 

k=l 

Now, we would like to apply Theorem 14.81 (for the case where the root satisfies the 
convexity condition) and Theorem 15.31 (for the general case), as in the case of simple 
roots; however, the proximity of the multiplicity brings the additional cut-off function, 
(1 — ip){t\^ — ^°|), into play, and this depends on t. Therefore, the aforementioned 
results cannot be used directly. However, a similar result does hold, provided we 
impose some additional conditions, producing analogues of Theorems 14.81 and 15.31 in 
this case. 

Proposition 7.6. Let x ^ C^(]R"). Suppose Tfc(^), k = 1,2, satisfy the following 
assumptions on suppx-' 

(i) for each multi-index a there exists a constant > such that, for some 5 > 0, 

i,rkW + S7])]\ < C„(l + |r/|)-l"l, for small s > Q and \r]\ > 6 ; 

(ii) there exists a constant Cq > such that \d^Tk{C,^+Xu!)\ > C > for alluj G S""""^ 
and X > 0] in particular, each of the level sets 

AS';, = S, = {r/GM":rfc(e° + r/) = A} 

is non- degenerate; 

(iii) there exists a constant Ri > such that, for all X > 0, 

Furthermore, assume that A'^{^) satisfies the following condition: for each multi- 
index a there exists a constant Ca > such that 

(iv) we have the estimate 

|(9^[Aj(e° + sr])]\ < CaS-^{l + |r/|)"^"l"l , for small s>0 and \r]\ > 6 . 
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Finally, assume that ip G C^{{—6,S)) is such thatip{a) = 1 for \a\ < 6/2. Then, the 
following estimate holds for a// x G M", t > 0: 

I ^ r 

K2 e^(-S+-'=«W4(e)(l-V)(t|e-e°IM0rfe <C{l + ty-\ (7.14) 
k=i -^^^ 

forj > n— where'-f := supx>ol{^x{'^k)) , if^kiO satisfies the convexity condition; 
and for j >n-^, where 70 := sup;,>o 7o(SA(rfc)), if it does not. 



Remark 7.7. Conditions (i) , (ii) and (iv) appear and are satisfied naturally when 
roots Tfc(^) are homogeneous functions of order one — for example, the wave equation, 
or for homogeneous equations. 



Remark 7.8. Assumption (iv) is needed because Aj{^) has a singularity at so we 
must ensure we are away from that — this is the role of the cut-off function {l—ilj){\r]\) 
in this proposition; 

Remark 7.9. ^45 usual, for example in the convex case, taking j = n — we get 
the time decay estimate 

I Left hand side of ^^TJ^ \ < C{1 + 1)" 



n-l 



Proof. As before, cut-off near the wave front: let k G C^(M") be a cut-off function 
supported in B{0,r). Then, consider 

^ r 

h{t,x) := ^ / e^(-«+-'=(5)*)4(0(l - - eMO<t''x + VniO) d^, 

k=i -^^^ 

and 

hit,x):=f2 [ e^(-«+-'=(«)*)4(0(l-^)(t|e-e°l)x(0 
k=i 

{l-n){t-'x + VTk{0)dC 

Away from the wave front set: First, we estimate l2{t,x); we claim that 

\l2{t,x)\ <Cr{l+ty-'' for alH > 0,x G M''. (7.15) 
In order to show this, we consider each term of the sum separately. 



I^{t,x) 



and imitate the proof of Lemma [4. 101 (in which the corresponding term was estimated 
in TheoremSSD, but noting that in place of gni^) G C(^(M") we have (l-?/')(t(^-^°)), 
which depends also on t; in particular, this means that care must be taken when 
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carrying out the integration by parts when derivatives fall on (1 — 'ip){t\C, — To 
take this into account, use the change of variables ^ = + t^^r]: 



Integrating by parts, with respect to rj gives 



where P* is the adjoint operator to P = . | ^_ ^ j|^o - 1 1 2 ■ V^; this integration by 

parts is vahd as \t~^x + (Vgrfc)(^° + t^^f])] > r > 0, in the support of (1 — k) (t~^x + 
Vrfc(^° + t~^r])^ . For suitable functions / = f{ri; x, t), and ^ = + t^^r], we have 



t-^x+{v^T,m f 



V,-(V5r,)(0 , , t-^x + (Vgr,)(0 



t\t-^x + (V5rfc)(0r ^1^-'^ + (V5^fe)(OI 

2(t-ix + (Vgrfe)(0) ■ [V,[(Vgr,)(0] ■ {f-'x + (Vgrfc)(0)] 

z|t-ia;+(Ver,0(OI' 



Comparing this to fl4.16l) . observe that the first and third terms have one power of t 
fewer in the denominator due to the transformation; this is critical in this case where 
we are approaching a singularity in Aj{$,^ + t~^ri) when t — > 00. By hypothesis 
for T) in the support of the integrand of x), we get 

V,.|(V,rOK« + <c-.(l + |„|)-'; 



\t-^x + {V;n)(,e + t-'q)\ 



2 



thus, we have 

|P7l<a[(l + |r/|)'Vl + |V,/|]. 

In Lemma 14. lU^ we carried out this integration by parts repeatedly in order to 
estimate the integral. Here, however, note that differentiating (1 — ■?/')(|?7|) once is 
sufficient: by definition of ip{s), 

d,Mi - mv\)] = -^^{dsmvD 

is supported in | < |?7| < 1, so 

|5.J(l-^)(l^l)]l<Cli>M>3/4(r/), 
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where li>|r?|>3/4('7) denotes the characteristic function of {77 G M" : 1 > \rj\ > 3/4}; 
hence, by hypothesis 



IV 



for large t we have 



A'{e+t-'n)\m{i-mri 



\t-^x + {v^nW + t-^v)? 

\x{e + t-'r^) 1 1 (1 - [r'x + Vr.ie + t-^)) dr^ 



< Cr 

< Crt^ 



f<l^l<i 

1 

<|r,|<l (1 + 1^1)^' 



(7.16) 



which is the desired estimate f l7.15p . 

On the other hand, if, when integrating by parts, the derivative does not fall 
on iIj{\t]\), we use a similar argument to that in the earlier proof; let us look at the 
effect of differentiating each of the other terms: in the support of -i/^d?]!), for each 
multi-index a and t > 0, 



IV 



• |a^[A^'(^° + < CaPi^ + 1^71)"^'"'°' by hypothesis 

• l'9^[x(^° + i"^^)]| < Ca{l + |r/|)"l°l: for a = 0, take C„ = 1; for |a| > 1, note that 

d^ixie + t-'v)] = ri"i(9^"x)(e° + t-'v) , 

and that {d^x){^'^+t~^v) is supported in < |^°+t"^?7| < 2N, so r ^ < CN^^o\r]\-^] 

• |(9°[(1 - K){t-^x + (V5rfc)(^o + t-^r]))]\ < + 1^1)"'"': obvious for a = 0; for 
|a| > 1, note 

= -{d^K){t-'x + VMO)d^[i'^^n){e + r^r/)] , 

which yields the desired estimate by hypothesis ( 
Summarising, this means 

1(1 - mv\)d^ [A'jie + t-'v)xie + t-'v)a - {t-'x + iv^n)ie + t~'v))] \ 

<C,(l + |77|)-^'"l"ltn|^l>3(r/). 

So, repeatedly integrating by parts we find that either a derivative falls on {l — ip){\ri\) 
(in which case a similar argument to that in fl7.16p above works) or we eventually get 
the integrable function CP{1 + |?7|)^"~^l|,)|>3/4(77) as an upper bound; in either case, 
we have fl7.15p . 



103 



On the wave front set: Next, we look at the term supported around the wave 
front set, x). As in the case away from the wave front, set ^ = C,^+t'~^ri: consider, 
for A; = 1,2, 

I^{t,x) := e'"-«° [ e'(*"'"-''+"^(«''+*"''')*)4(e° + r^r/)(l -7/;)(|r/|) 

As in the proof of Theorems 14.81 and 15.31 let {"^ i{'r])}e=i be a conic partition of unity, 
where the support of "^eir]) is a cone K^, and each cone can be mapped by rotation 
onto Ki, which contains e^i = (0, . . . , 0, 1). Then, it suffices to estimate 

^Mxie + t-'vHt-'x + iV^rk)ie + t~'v)) dv , 

for k = l,2. 

Let us parameterise the cone Ki: it follows from hypothesis (ii) that each of the 
level sets 

is non-degenerate; so, for some U C M"~^, and smooth function hk{t, X,-) : U ^ M, 

K, = {i\y,\hkit,\,y)):\>0,yeU} . 

If Tk{^) satisfies the convexity condition, then is also a concave function in y. Now, 
we change variables r] i— > {\y, Xh^it, A, y)) and will often omit t from the notation of 
hk since the dependence on t will be uniform. We obtain: 



Ju 



(1 - ij){X\{y, hk{\, y)miiKy, hk{\, y)))x{e + r'X{y, h^iX, y))) 

K{t-'x + (V5r,)(e° + t~'Xiy, hiX, y))))^^dXdy, (7.17) 

where we have used rfc(^° + t~~^{Xy, A/ifc(A,|/))) = t~^A. As in the earlier proofs, we 



ensure Xn is away from zero in the cone — this requires hypotheses (i) and (iii)). So 



in the general case, we can write this as, with x = t ^x, X = Xxn = Xt 

poo p 

r" / / e*^^"(*"'^"'^'-^+*^''^^(^'J^)+^"')A^^(^° + r^A(?/,/ifc(A,?/))) 



Ju 



;i - ^)(A|(2/, h,{X, y)mi{X{y, h,{X, y)))x{e + t~'X{y, h,{X, y))) 

K{t-'x + (V^r,) ie + t-'Xiy, h,iX, y)))) dXdy . 
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If the convexity condition holds, then, as in the proof of Theorem I4.8[ we have 
the Gauss map 



Vc[rfc(e° + r^C)] _ (Ver,)(e° + r^C) 



and, as before, can define Zk{X) G t/ so that 

n^,{zkiX), hk{\, z{\))) = -x/\x\ . 

Then, 

— = -Vyhk{X,z{X)). 



So, in this case, (17.171) becomes: 



in 



'0 Ju 



+ t-'Xiy, h,iX, y)))K{x + iW^r,)ie + r h,{X, y)))) 



dr] 



d{X,y) 



dXdy, 



Let us estimate this integral in the case where the convexity condition holds. We 
have: 



The same argument as in the earlier proof (which uses hypothesis (ii)), shows 

dr] 



< CX 



n-l 



d{X,y) 

The constant C here is independent of t; 

Now, with Aiiu) = Ai{iy)x{iy)K{x + {V^Tk){iy))'^i{X{y,hk{X,y))), where u = + 
t~^X{y,hk{X,y)), we have 



JXx,4-{y-z{X))-\7yh^{X,z{^))+h^{X,y)+hk{X,z{^))] 



U 



t-^X^A'ie + t-'X{y, hk{X, y)))(l - ^){X\{y, /i,(A, y))|) dy 



A"-^-^' dX . 



Now, applying Theorem 14.11 — this may be used due to the properties of ^4*^(0 
Tk{^) stated in hypotheses (iv) and[(i)] — we find that 



,^Xx^[^{y^z{X))■\7yhkiX,zW)+hk{X,y)+h^:{X,z{m 



t-^X^A^ie + t~'X{y, hk{X, y)))(l - ^)(A|(y, h,{X, 2/))|) dy < CX^-^xiX) , 

where x(A) is a compactly supported smooth function that is zero in a neighbour- 
hood of the origin. 
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• Hence, 

/■oo 

Jo 

Finally, the general case without convexity can be estimated in a similar way, with 
the necessary changes used in the proof of Theorem 15.31 to account for the change in 
the phase function — in particular, the use of the Van der Corput Lemma, Lemma [5?T1 
in place of Theorem 14. 1[ This completes the proof of (17.141) . □ 



Using Proposition I7.6[ it is clear that 

2 

e"-«(l - ^P) (tie - e\)x{0 Yl Oe'-'^^^'m 



k=l 



< cil + ty 



if the roots satisfy the convexity condition, and 

2 



k=l 



L°°(M5) 

< c{i + ty 



otherwise. In comparison to (I6.16p . here we have L^-norms on the right hand sides, 
since x is a cut-off function to bounded frequencies. 

Finally, we must consider the case where L roots intersect; the above proof can 
easily be adapted for such a case, giving corresponding results. 

— L^ estimates: For the L^-estimates on the support of (1 — 'V')(t|e — e°|)x(0 



only need assumption (iv) of Proposition 17.61 with a = for the amplitude, namely 
that 



A^(e° + sri)\ < Cas-^{1 + \r]\)-^ , for small s > and \v\ > 5 . 



(7.18) 



Then, for the left hand side of (17.141) . we have 

j2 [ e'(-«+-'=(«)*)4(o(i-v^)(tie-e°i)x(o/(0^^e 
k=i 

2 

Ee^-'=(«)*4(0(i - - e\)x{om 



L2(RS) 



k=l 



L2(]R") 



<l|i^'(l + l^l)"V(^° + r^)lk.(M.), 



where we used Plancherel's theorem, (17.181) . and the notation s = t ^, = C,^ + t ^t], 
so that 7] = t{C, — e°). Then we can easily estimate 
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where D — Dq is & Fourier multiplier with symbol ^ — So, we finally obtain the 
estimate 



In the case of equations with homogeneous symbols (like for the wave equation), when 
roots are homogeneous, we have C,^ = 0, so that the right hand side becomes just the 
norm in the corresponding homogeneous Sobolev space. 

Due to the earlier bound near the multiplicity, we can combine the results with 
the interpolation Theorem 16.41 

8 Examples and extensions 

Theorem 12.181 gives estimates for operators provided the characteristic roots satisfy 
certain hypotheses. However, in order to test the validity of such an estimate for an 
arbitrary linear, constant coefficient m}^ order strictly hyperbolic operator with lower 
order terms, it is desirable to find conditions on the structure of the lower order terms 
under which certain conditions for the characteristic roots hold. For the case m = 2, 
a complete characterisation can be given, and some extension of this is discussed in 
Section ISTTl However, for large m, it is difficult to do such an analysis, as no explicit 
formulae are available in general; nevertheless, certain conditions can be found that 
do make the task of checking the conditions of the characteristic roots, and these are 
discussed in Section 18. 2^ where a method is also given that can be used to find many 
examples. Finally, in Section 18. 5^ we give a few applications of these results. 

8.1 Wave equation with mass and dissipation 

As an example of how to use Theorem I2.18[ here we will show that we can still 
have time decay of solutions if we allow the negative mass but exclude certain low 
frequencies for Cauchy data. This is given in (18.11) below. In the case of the negative 
mass and positive dissipation, there is an interplay between them with frequencies 
that we are going to exhibit. The usual non- negative and also time dependent mass 
and dissipation with oscillations have been considered before, even with oscillations. 
See, for example, [HR03j and references therein. 




L2(IR5) 



<C|||D-DorVI|L2(RS)- 
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Let us consider second order equations of the following form 



d^u - c^Au + 6dtu + fiu = 0, 
u{0,x) = 0, ut{0,x) = g{x) . 

Here S is the dissipation and fi is the mass. For simplicity, the first Cauchy data 
is taken to be zero. The general case when both Cauchy data are non-zero can be 
treated in a similar way. Let us now apply Theorem 12.181 to the analysis of this 
equation. The associated characteristic polynomial is 

- C^\^\^ - t6T - fi = , 

and it has roots 

r±(0 = |± Vc2|eP + /i-5V4. 

Now, we have the following well-known cases, which also correspond to different cases 
of Theorem \2AM 

• 6 = fi = 0. This is the wave equation. 

• 6 = 0, fi > 0. This is the Klein-Gordon equation. 

• fi = 0, 6 > 0. This is the dissipative wave equation. 

• 5 < 0. In this case, Imr-_(^) < | < for all ^, hence we cannot expect any decay 
in general. 

• 6 > 0, fi > 0. In this case the discriminant is always strictly greater than —5^/4, 
and thus the roots always lie in the upper half plane and are separated from the 
real axis. So we have exponential decay. 

Here is the main case for us, where we can show an interesting interplay between 
negative mass fi < and how it is compensated by positive dissipation 6 > for 
different frequencies: 

• dissipation 6 > 0, mass yU < 0. In this case, note that Imr_(,^) > if and only if 
c^KP + A* > 0, i.e. Imr_(^) = for |,^| = y/—fi/c. Therefore, the answer depends 
on the Cauchy data g. In particular, if supp^ is contained in {c^l^P + A* > 0}, 
then we may get decay of some type. More precisely, let B{0,r) denote the open 
ball with radius r centred at the origin. Then we have: 

— if (7 is such that supp^fl -8(0, ^^^) 7^ 0, then we have no decay; 

— if there is some e > such that supp^ C ]R"\i?(0, ^^^^ + e), then the roots are 
either separated from the real axis (if 5 > 0), and we get exponential decay, 
or lie on the real axis (if 6 = 0), and we get Klein-Gordon type behaviour 
(since the Hessian of r is nonsingular) . 
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— if, for all g, supp^ C M" \ 5(0, -^^-^) = ||^| > ^^^Tj ' then again we must 
consider 5 = and 5 > separately. 

If (5 = 0, then the roots lie completely on the real axis, and they meet on the 
sphere |^| = /i/c. It follows from fl2.17p (which is justified in Proposition 
17. 4p with L = 2 and i = 1 that, although the representation of solution as 
a sum of Fourier integrals breaks down at the sphere, the solution is still 
bounded in a (l/t)-neighbourhood of the sphere. In its complement we can 
get the decay. 

If 5 > 0, then the root r_ comes to the real axis at |,^| = in which case 
we get the decay 

||n(^,■)IU.<C(l + ^)"(^i)||^7||iP. (8.1) 

Indeed, in this case the order of the root r_ at the axis is one, i.e. estimate 
(EED holds with s = 1. Here 1/p + 1/g = 1 and 1 < p < 2. Note also that 
compared to the case of no mass when i = n, now the codimension of the 

sphere |^ e : |.^| = -^^-^j is £ = 1. We can apply the last case of Part II 

of Theorem 12.181 with L = 1 and s = i = 1 which gives estimate (18.11) . 



8.2 Higher order equations 

Let us now derive a simple consequence of the stability condition of Imrfc(.^) > 0, 
for all /c = 1, . . . , m and ^ G W^, for the coefficient of the D^~^u term in (11. ip . In 
fact, this coefficient plays an important role for higher order equations and can be 
compared with the dissipation term in the dissipative wave equation. 

Let L = L{Dt, Dx) be an m}^ order constant coefficient, linear strictly hyperbolic 
operator such that Imrfc(^) > for all = 1, . . . ,m and for all ^ G M". Recall that 
the characteristic polynomial corresponding to the principal part of L is of the form 

m 
k=l 

where the Pk{0 homogeneous polynomials of order k. Then, by the strict hyper- 
bolicity of L, Lm has real roots ipi{^) < ip2{0 — ' ' ' — V^m(0 (where the inequalities 
are strict when ^ 0). By the Vieta formulae, observe that 

m 
k=l 

On the other hand, the characteristic poljTiomial of the full operator is 

m m—1 

^r,0=r'" + J]P,(Or'"-'= + $^ J] c^./^V = . (8.3) 

k=l j=0 \a\+l=j 
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In particular, the coefficient of the r"* ^ term is 

m 

Pi(0 + co,m-i = - J]r,(0, (8.4) 

fc=i 

where the rfc(^), k = 1, . . . ,m are the roots of fl8.3l) . Comparing (18.21) and fl8.4l) . we 
see that Im ( ^^^^ rfc(^)) = — Imco,m-i- Therefore, since Imrfc(^) > for all k = 
1, . . . ,m and ^ G M*^, it follows that Imco,m-i < 0, or, equivalently, Reico^m-i > 0. 
Furthermore, if Imco,m-i = then it must be the case that Imrfc(.^) = for all ^ G 
and k = 1, ... ,171 since the characteristic roots are continuous. Hence we have shown 
the following: 

Proposition 8.1. Let L = L{Dt,Dx) be an m*^ order linear constant coefficient 
strictly hyperbolic operator such that all the characteristic roots Tk{C,), k = 1, . . . ,m, 
satisfy Imrfc(^) > for all ^ E W^. Then the imaginary part of the coefficient 
of D^~^u is non-positive. Furthermore, if in addition the {imaginary part of the) 
coefficient of D^~^u is zero then each of the characteristic roots lie completely on the 
real axis. 

If we transform our operator back to the form L{dt,dx), this result tells us that 
in order for the characteristic polynomial to be stable, that is for lmTk{$,) > for 
all A; = 1, . . . , m, ,^ G M.^, it is necessary for the coefficient of <9™~^n to be non- 
negative; this is the case for the dissipative wave equation. In some sense this may 
be interpreted as a higher order dissipation, since it is necessary for the characteristic 
roots to behave geometrically like those of the wave equation with a dissipative term, 
where they lie in the half-plane Im^; > and lie away from Im2; = for large |^|. 

In the next section, we look at the case where characteristic roots must lie com- 
pletely on the real axis. First, though, let us consider one case where a root lies 
completely on the real axis but the coefficient Co,m-i is nonzero, Co,m-i 7^ 0. 

Consider a constant coefficient strictly hyperbolic operator of the form 

Lm{dt, <9x) + Lm-i{dt, dx) + Lm.^2{dt, d^) = 0, (8.5) 

where = Lr{dt,dx) denotes a homogeneous operator of degree r with real coef- 
ficients. This is an example of a hyperbolic triple, which will be discussed in more 
generality in Section [8731 Furthermore, assume that Lm~i is not identically zero. Let 
r(^) G M be a characteristic root of (18. 5p which lies completely on the real axis. So, 
denoting as usual D^. = —idx^, Dt = —idt, we have that t{^) is a root of 

Lm^r, - iL^-i{,r, C) - L^_2(r, = 0. 

This means that Lm_i{^, r(^)) = 0, and so t{^) is homogeneous of order 1, and thus 
for such roots Theorem 12.181 applies to yield results similar to those described in 
Section 11.21 
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8.3 Hyperbolic triples 

We now turn to the case when all the characteristic roots lie completely on the 
real axis. This section is devoted to showing some more examples of appearances 
of real valued non-homogeneous roots and some sufficient conditions for this. In 
order to study this case we first recall some results of Volevich-Radkevich |VR03j 
on hyperbolic pairs and triples. Throughout this section only, Lr(r, ^) denotes a 
homogeneous polynomial in r and ^ = {^i, . . . , ^n) of order r such that Lr{T, i^) has 
real coefficients. 

Definition 8.2. Suppose Lm = Lm{T,^) and Lm.^i = Lm.^i{r,^) are homogeneous 
polynomials as above. Furthermore, assume that the roots of Lm, Ti{C,) , . . . , Tm{0 > 
and those of Lm-i, <yi{C)i ■ ■ ■ i <^m-i{C)! ^'^^ real-valued {in which case we say Lm and 
Lm-i are hyperbolic polynomials). Then, {Lm,Lm-i) is called a hyperbolic pair if 
{possibly after reordering) 

ri(0 < ^i(e) < 7-2(0 < ■ ■ ■ < rm^iiO < ^m-i(0 < r^(0- (8.6) 

If, in addition, the roots of Lm, Lm-i are pairwise distinct for ^ ^ {in which case 
they are called strictly hyperbolic polynomials) and the inequalities in \8.6^ are all 
strict, then we say {Lm, Lm-i) is a strictly hyperbolic pair. 

Definition 8.3. Let 

Lm = Lm{T,^) , Lm-1 = Lm-l{T, ^) , Lm-2 = -^m-2(T, 

be {homogeneous) hyperbolic polynomials. If {Lm,Lm-i) and {Lm-i, Lm-2) both 
hyperbolic pairs then we say that {Lm, Lm-i, Lm-2) is a hyperbohc triple. //, in 
addition, all the polynomials and all the pairs are strictly hyperbolic {in the sense of 
Definition 18.21) then {Lm, Lm-i, Lm-2) is called a strictly hyperbolic triple. 

Theorem 8.4 ( |VR03j ). Suppose that {Lm, Lm-i, Lm-2) is a strictly hyperbolic triple. 
Then Lm{T,^) + Lm-i{T,C,) + Lm-2{T,C) ^ /(^"i^ '^^^ lmr<0. Furthermore, any two 
of the polynomials Lm, Lm-i, Lm-2 have no common purely imaginary zeros. 

We also recall a theorem of Hermite (see, for example, [NisOO] ): 

Theorem 8.5. Suppose pm{z), pm-i{z) are real polynomials of degree m,m — 1, 
respectively, and that all the zeros of p{z) = Pm{z) — ipm-i{z) He in the upper half- 
plane {that is, if p{z) = then Imz > 0). Then all the zeros of Pm{z) and Pm-i{z) 
are real and distinct. 

Now we will give some rather constructive examples of how non-homogeneous real 
roots may arise, and some sufficient conditions for this. 

Assume that L is of the form Lm{Dt,Dx) + Lm-2{Dt, Dx), where the L^ are as 
in Definition 18.31 and neither is identically zero. Suppose that there exists a homoge- 
neous operator of order m — 1, Lm-i{Dt, Dx), such that the characteristic polynomials 
Lm{T,^), Lm-i{T,C) and Lm-2{T,i) form a strictly hyperbolic triple. Then, by The- 
orem [8]U we have 

Lm{T,i)+Lm-l{T,i)+Lm-2{T,i)^^i0I Im T < . 
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Thus, by Theorem 18.51 all the zeros of Lmij,^) + Lm_2(r, ^) are real, but clearly 
non-homogeneous if Lm-2 ^ 0. So, using this construction, we can obtain examples 
of operators for which all the characteristic roots lie completely on the imaginary 
axis (so that «r(^) are real, which would be the notation for the rest of this paper), 
but for which we cannot automatically expect the standard decay for homogeneous 
symbols to hold. 

8.4 Strictly hyperbolic systems 

Our results can also be used to establish — decay rates for strictly hyperbolic 
systems. Let us briefly sketch the reduction of systems to the situation covered by 
results of this paper. Let 

iUt = A{D)U, f/(0) = f/o, 

be an m X m first order strictly hyperbolic system of partial differential equations. 
That is, the associated system of polynomials may be written as A{S,) = Ai{^)+Aq{C,), 
with Ai being positively homogeneous of order one in ^ and Ao(0 ^ 'S'°q(M"). If 
A{C,) is a matrix of first order polynomials, then Ao is constant. It is known that 
A[D) is hyperbolic if and only if det A{D) is hyperbolic (see e.g. Atiyah, Bott and 
Carding jABG] ). Moreover, if detAi{D) is strictly hyperbolic, then A{D) is strongly 
hyperbolic. 

Now, the strict hyperbolicity of A{D) means that the roots ^Pi{^), ■ ■ ■ ,(pm{0 '^^ 
equation det{(pl — Ai{$^)) = are all real and distinct away from the origin. Denote 
the roots of the equation det (rJ — A{$^)) = (which is an m}^ order polynomial in 
r with smooth coefficients) by ri(^), . . . , rm(0- Now, by analogy to the case of the 
^th QTi-(^Qj- scalar equation, we can, via perturbation methods, show that for large |^| 
the rfc(^) behave similarly to the (fikiO^ that they are distinct, analytic and belong 
to S'j'^Q(M"). For bounded |^| we will need similar regularity assumptions on the 
characteristic roots Tk{C,) as for the scalar equations. Furthermore, we assume that 
there exists Q e SIq{W) such that |det(5(OI > C > and such that 

g-iAg = diag(ri(0,...,r„(0)=:T. 

The existence of such Q is a very interesting question itself, especially in the presence 
of variable multiplicities, but we will not go into such details here. Now, we use the 
transformation U = Q{D)V, so that 

Ut = QVt =^ iQVt = AiD)QV =^ iVt = TV ; U{0) = QViO). 

This systems decouples into m independent scalar equations: 

dtVk = n{D)Vk, = 1, . . . , m, \4(0) = {Q~'U{0))k 

each of which is solved by 

Vk{t,x)= [ e*(^-«+^'=(«W\4(0,e)rfe- 
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Now, Q G S'°(M"), so it is a bounded map L'^ L'^, 1 < g < oo, and we can get our 
estimates for Vk as in the case of m**^ order scalar equations; thus, we can conclude 
that 

||f/||L. = WQVWl. < C\\V\\l. < CK{t)\\V\\Lr> = CK{t)\\Q-'U\\Lv < CK{t)\\U\\L. , 

where K{t) is as in Theorem 12.181 

8.5 Application to Fokker— Planck equation 

The classical Boltzmann equation for the particle distribution function / = f{t, x, c), 
where x, c G M", n = 1,2, 3, is 

{dt + c-V,)f = S{f), 

where S{f) is the so-called integral of collisions. The important special case of this 
equation is the Fokker-Planck equation for the distribution function of particles in 
Brownian motion, when the integral of collisions is linear and is given by 

n 

S{f) = Ve ■ (C + Ve)/ = ^^.(C*^ + 9c,) f. 

k=l 

In this case the kinetic Fokker-Planck equations takes the form 

(n \ n 

dt + J2 ^''^^^ I ^' = XI ^^fe ^^'^ + 
k=l / k=l 

The Hermite-Grad method of dealing with Fokker-Planck equation consists in de- 
composing f{t,x, ■) in the Hermite basis, i.e. writing 

f{t,x,c) = Y] ^m^{t,x)4j"'{c), 

\a\>0 

where ip°'{c) = (27r)~"/^(— c?c)" exp(— ■^) are Hermite functions. They are derivatives 
of the Maxwell distribution ip'^ which annihilates the integral of collisions and form 
a complete orthonormal basis in the weighted Hilbert space Z/^(M") with weight 
w = This decomposition yields the infinite system 

x) + dx,Tnpj^ek{t,x) + \(5\mj3{t, x) — 0. 

The Galerkin approximation of the solution / is 

f{t,x,c)= V -m^{t,x)r{c), 
^ — ^ cd 

0<\a\<N 

^Thus, the convergence of the series of such decomposition is understood as a convergence of the 
decomposition with respect to a basis in a Hilbert space. 
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with m{t^x) = {mp(t,x) : < < A^} being the unknown function of coefficients. 
For m(t, x) one obtains the following system of equations 

Dtm(t, x) + AjDx_.m{t, x) — iBm{t, x) = 0, 



where i? is a diagonal matrix, -Bq,/? = |«|5q:,/3, and the only non-zero elements of the 

a—ej,a a- 

I J = aj, a J 



matrix Aj are a" = aj, a^~^^^'°' = 1. Hence, the dispersion equation for the 



system is 

P(r, = det(r/ + ^ A,^, - iB) = 0, (8.7) 

i 

which we will call the A^*'* Fokker-Planck polynomial, and we have, in particular. 
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P{t, 0) = det(rJ - iB) = t JJ(r - jip , (8.8) 

j=i 

for some powers 7j > 0. Properties of this polynomial P{t,^) have been extensively 
studied by Volevich and Radkevich in |VR04] . who gave conditions and examples 
of situations when Imrj(^) > 0, for all ^ 7^ 0. They also described more general 
(necessary) conditions in terms of coefficients of P. See also |VR03t IZR04j . In our 
situation here we have to take additional care of possible multiple roots, as is done 
in Theorem 12. 16[ 

From formula fl8.8p it follows in particular that there is a single characteristic root 
at the origin. Let M = Y[f=i ■ 

Let us examine the structure of the operator P{t,C,)- It is a polynomial of order 
m which can be written in the form 

m 
j=0 

with Pj being a homogeneous polynomial of order j. Moreover, we have 
Po = 0, Pi = Mt, P2 = M^^r^-M|e|^ 

k=2 

The case n = 1 was considered in |VR03] . where one has M = N\ 

Let P(r(^), ^) = 0, where r(0) = is the simple root at the origin. Differentiation 
with respect to r yields ^(0) = 0. Differentiating again we get 



as 



(0) = 22/„. 



So, for small frequencies we obtain the decomposition 

Imr(0 = 2|eP + ... + c(logm)||eir + ..., 
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where 

m = 1 + 7i + . . . «i CnN, 

and 1 1^1 1^ denotes a fourth order polynomial in ^. We also easily have a rough estimate 
for M of the form 

<M < (iV!)^, in > 2). 
It follows then that for small frequencies we get the estimate 

|m(t, x) I < C(l + + Ce-^(^)*, 

where, in general, it may be that s{N) — > as ^ oo. For medium frequencies 
we get exponential decay in view of the result of Theorem 12.1^ also in the case 
when there are multiple characteristics, where we can use Theorem 12. 2[ Here, there 
is an additional polynomial growth with respect to time caused by the resolution 
procedure of Section 17.11 but this is compensated by the exponential decay given by 
characteristics with strictly positive imaginary part (see Theorem 12. 2p . 

Let us discuss the situation with large frequencies. For operators of general form, 
away from points where roots coincide, the roots are analytic. For large |^|, per- 
turbation arguments of Section [3] give properties of roots Tk{^) related to v^fc(O) the 
characteristics of the principal part. Here Tk{C) ^"^^ (pkiO defined as roots of 
equations P(r, = and its principal part Pm{^p,0 = O5 respectively. Let K be the 
maximal order of lower order terms. Then we can summarise the following properties 
of P established in Section [31 

• there are no multiple roots for large ^; 

• |9^"rfc(0l<C(i + |e|)'-i"i,i.e. neS\ 

• the exits ifk such that |(9"rfc(0 - d^MOl < C{1 + for all ^ 
and all multi- indices a; 

• Since (f>k are real- valued, we get IniTfc G 5*^^+1-™. In particular, Imr^ G S^. 

The statements above are obtained by perturbation arguments and rely on the 
strict hyperbolicity of the principal part. However, this does not have to be the 
case for polynomials P that we obtain in the Galerkin approximation. Moreover, in 
general, it might happen that Imrfc(^) ^ as |^| — 00, the case which is discussed 
in Section 16.81 To avoid these problems we impose the condition of strong stability. 
First, we will say that P(r, ^) is a stable polynomial if its roots r(^) satisfy Im r(^) > 
for all ^ e M", and if Imr(^) = implies ^ = 0. Then we will say that P(r, ^) is 
strongly stable if, moreover, Imr(^) = implies ^ = and Rer(^) = 0, and if its 
roots r(^) satisfy liminf|^|^oo Imr(^) > 0. Thus, the condition of strong stability 
means that the roots t{^) may become real only at the origin of the complex plane 
at ,^ = 0, and that they do not approach the real axis asymptotically for large ^. 

In Section [873| as well as in [VR03[ I VR04j . there are several sufficient conditions 
for the stability of hyperbolic polynomials. In this case we have a consequence of 
Theorem 12.161 and Remark 12.171 in the form of estimate fl2.15p : 
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Corollary 8.6. Let P be a strongly stable polynomial with characteristic roots with 
non-negative imaginary parts. Let 1 < p <2 and 2 < g < oo be such that ^ + ^ = 1. 
Then the solution to Cauchy problem \2.1\) satisfies dispersive estimate Ii2.15\) . i.e. 
we have 



Till 11 ■'''^1 ^ > 



i=o 



3 ll^yj^P + l^l+'-^J 



with Np > — ^) for 1 < p < 2 and Ni > n for p = 1. 

From this, we can conclude the following estimates for solution to the Galerkin 
approximations of Fokker-Planck equation: 

Theorem 8.7. If the N^^ Fokker-Planck polynomial P in (8.1) is strongly stable, we 
have the estimate 



||/Ar(t, X,c)||ioo(]jn)2,2(H 

with w = exp(— |cp/2) and e{N) > 0. 



<C(l + t)-"/' + C;ve 



-e(N)t 



Here the constant C is independent of A^, but, in general, we may have asymp- 
totically that e{N) — > as — oo. The validity of the assumption of Theorem 18.71 
for all N is an open problem. 
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